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ABSTRACT 



The object of this thesis is to study the spin-spin correlations for top quark productions and 
their decay at polarized e + e~ linear colliders. For the analyses of the spin-spin correlations 
in top quark production, we introduce a generic spin basis and show which spin basis is 
the most optimal. We discuss the spin-angular correlations for the top quark decay, and 
demonstrate that the directions of electro-weak decay products of polarized top quarks are 
strongly correlated to the top quark spin. We also explain how to measure the top quark 
spin using decay products. The spin correlations using the generic spin basis are extensively 
studied at the next to leading order of QCD for top quark productions at e + e _ colliders. The 
radiative corrections induce an anomalous 7/Z magnetic moment for the top quarks and allow 
for real gluon emissions. Therefore they might change the leading order analyses. We show 
that higher order effects on the top quark spin orientation are very small. The final results 
are that the top (or anti-top) quarks are produced in an essentially unique spin configuration 
for a particular spin basis even after including the Q(a s ) QCD corrections. Our results imply 
that measuring the spin correlation for top quarks is very promising to reveal new aspects of 
the Standard Model, as well as new physics beyond the Standard Model. 



"Truth emerges more readily from error than confusion. 

-Francis Bacon- 
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Chapter 1 



Introduction 



High-energy e + e~ colliders have been essential instruments to search for the fundamental 
constituents of matter and their interactions. Merged with the experimental observations at 
hadron colliders, the structures and properties of the "fundamental" particles are adequately 
described by the Standard Model. The matter particles, which are leptons and quarks, can 
be classified into three generations with identical symmetries. It has been known that the 
electro-weak and strong forces are described by gauge field theories. The strong force is de- 
scribed by the SU(3) C gauge theory, which is called Quantum Chromo Dynamics (QCD) (]]]. 
The vector bosons of this model are grouped into eight colored massless gluons which mediate 
the strong interactions. The QCD has only one arbitrary or free parameter, Aqcd- Neverthe- 
less, it explains, all hadronic dynamics including confinement, asymptotic freedom, nucleon 
structure, and baryon and meson spectroscopy. On the other hand, the electro-weak theory 
has many arbitrary parameters and is described with the gauge group of SU(2) £ ® U(l)y 0, 
where SU(2) L is the weak isospin group, acting on left-handed fermions, and U(l)y is the hy- 
percharge group. At "low" energy (at the scale ~ 250 GeV), the SU(2)x®U(l)y symmetry is 
spontaneously broken. This mechanism is referred to as Higgs mechanism [|J, which induces 
mass of gauge bosons, W ± ,Z, and other fermions. The heavy gauge bosons, W ± ,Z, were 
discovered at CERN in the middle of 1980's. Thanks to these impressive successes, the gauge 
theory with the group SU(3) C <g> SU(2) L <g> U(l) y has been called the "Standard Model". 

The Standard Model has been tested and checked through precision experiments with 
remarkable accuracy during the past few decades. With those numerous data, the Standard 
Model has been successful in describing the properties of particles and the structure of their 
basic interactions. Despite this success, however, the Standard Model has not been considered 
as the final theory of nature. As a matter of fact, the Higgs mechanism is still hypothetical and 
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has not been established experimentally. The properties of the latest member of the matter, 
top quark, have not yet been fully revealed. To make a further progress and to step on to a new 
stage, we must investigate both the "strengths" and the "weaknesses" of the Standard Model 
and must understand it in detail. There are two complementary and sometimes mutually 
related strategies to explore the new aspects of the Standard Model and to reveal signals 
of new physics beyond it. The first is to look for "new" phenomena at high energy collider 
facilities, the e + e~ colliders LEP II, the ep colliders HERA, pp colliders TEVATRON, the 
pp colliders LHC and e + e~ linear colliders JLC The high energy colliders bring us a good 
chance to have a clue to new physics. The second is to utilize good probes, which might be 
Higgs particle and top quarks. Precision studies of these heavy particles also may reveal new 
aspects of the Standard Model and clues to the physics beyond the Standard Model, since 
these particles are considered to have important information on the mechanism of breaking 
of the electro-weak symmetry. In this thesis, we will focus on the top quark as a fascinating 
probe. A high energy e + e~ linear collider provides a very impressive tool to carry out detailed 
studies of top quarks. 

The first stage of top quark physics was opened with its discovery at FNAL [|], |5| . The top 
quark with a very large mass 173.8 ± 5.2 GeV makes us expect that the top quark plays a 
special role in particle physics. It might contain a hint of new physics beyond the Standard 
Model. The presence of Higgs boson and its Yukawa coupling to fermions play a fundamental 
role in the Standard Model, and the Yukawa coupling is proportion to the fermion mass. 
Therefore the top quark will likely offer a possibility for its direct measurement, which has not 
yet been established experimentally. Among many other possibilities, for instance, T or CP 
violation in the top quark sector and the possibility of the presence of anomalous couplings 
have been discussed in many papers |7[]. The very heavy top quark will provide us with a 
unique opportunity to understand nature deeply. 

Now, the following fact is typical to heavy top quarks that the top quark rapidly decays 
through electro-weak interactions before hadronization effects || |J come in. Therefore there 
are significant angular correlation between the decay products of the top quark and the spin 



of the top quark [|X0| | . This implies that if the produced top and anti-top quarks correlate their 
spins, there will be a sizable angular correlations between all the particles, both incoming and 
outgoing, in these events. There are many papers on the angular correlations for top quark 



events [[IT], |T2[ produced at both e + e colliders [13] and hadron colliders [13]. In most of these 
works, the top quark spin is decomposed in the helicity basis. For ultra-relativistic particles, 



this approach is appropriate. However, in general, the helicity basis is not a unique choice for 
massive particles. It may depend on the energy (speed) of the produced quark which spin basis 
is the most optimal to study the spin correlations. Mahlon and Parke [IS| have considered 



various decompositions of the top quark spin which result in a large asymmetry at hadron 
colliders. Parke and Shadmi [|16J extended this study to the e + e~ annihilation process at the 
leading order in the perturbation theory and found that the "off-diagonal" basis is the most 
efficient decomposition of the top (anti-top) quark spin. In this spin basis, the top quarks 
are produced in an essentially unique spin configuration. This result is of great interest and 
importance from both theoretical and phenomenological viewpoints. However, it is of crucial 
importance to estimate the radiative corrections to this process which are dominated by QCD. 
The radiative corrections, in general, add two effects to the leading order analysis: the first 
is that a new vertex structure (anomalous j/Z magnetic moment) is induced by the loop 
corrections to the tree level vertex, the second is that a (hard) real gluon emission from 
the final quarks can flip the spin and change the momentum of parent quarks. Therefore, 
compared to the radiative corrections to physical quantities which are spin independent, it is 
possible that spin- dependent quantities may be particularly sensitive to the effects of QCD 
radiative corrections. The analytical study of QCD radiative corrections to heavy quark 
production was pioneered in Ref. [|17] (see e.g. Ref. |18) for a recent article). Polarized heavy 
quark productions, in helicity basis, have been calculated by many authors^, |(J. Recently, 
the polarized cross-sections for not only longitudinally but also transversely polarized heavy 
quarks have been investigated |H|] at the next to leading order of QCD corrections. Meanwhile, 
polarization phenomena in top quark productions near threshold have also been analyzed [ [22 



The main part of this thesis is a discussion on the polarized top quark production at the 
QCD one-loop level. We focus on the issue of what is the optimal decomposition of the top 
quark spin for e + e~ colliders. (The physics of top quark production at muon colliders and 
e + e _ colliders is identical provided the energy is not tuned to the Higgs boson resonance.) 
We calculate the cross-section in a "generic" spin basis which includes the helicity basis as a 
special case. 

The thesis is organized as follows. In Section 2, we review some techniques for calculating 
helicity amplitudes, since these techniques are useful and effective to investigate, in particular, 
the spin dependent amplitudes. One can see that these calculational methods are indispensable 
techniques to analyze the spin correlations. In Section 3, we discuss the spin-spin correlations 
in top quark pair productions and the spin-angular correlations in the decay process of top 
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quarks at the leading order in perturbation theory. The measurement of top quark spin is 
also discussed. We demonstrate that top quark pairs are produced in an essentially unique 
spin configuration in polarized e + e~ linear colliders, and that the angular distribution of top 
quark decay products is strongly correlated to the direction of top quark spin axis. In section 
4, we examine the QCD corrections to the polarized top (anti-top) quark production, and 
present our analytic calculations of one-loop corrections to the polarized top quark production 
in a generic spin basis. We also give the numerical results both in the helicity, beam-line 
and the off-diagonal bases. We compare the full one-loop results with those in the leading 
order and soft gluon approximations. The final chapter contains a summary. In Appendices, 
the notations, conventions and some useful formulae for the helicity amplitude we use in this 
thesis are included. The phase space integrals which are needed in Section 4 are summarized. 
The unpolarized total cross-section for top pair production, using our results, is also given as 
a cross check. 



Chapter 2 

Spinor Helicity Method 



With the progress of precision in experimental data, we must calculate higher order Feynman 
diagrams in the perturbation theory in order to refine the theoretical predictions. For scat- 
tering process with more than two particles in the final state, the spin projection operator 
methods quickly become unwieldy This difficulty is overcome by employing the spinor helicity 
method, which is a powerful technique for computing helicity amplitudes for multi particle 
process involving spin 1/2 and spin 1 particles. 

The technique of helicity amplitudes in the high-energy limit (which corresponds to mass- 
less limit) was pioneered in the paper by Bjorken and Chen p3| , and later was extended to 
the massless gauge theories by the CALKUL collaboration [p4]. Xu, Zhang and Chang p5 



had proposed a new definition of the polarization vector of the gauge boson based on the 
helicity spinor basis and the decomposition of amplitudes into gauge-invariant subsets. They 
generalized the helicity amplitude techniques developed by the CALKUL collaboration to the 
non-abelian case. Meanwhile, a technique |H| , which utilizes the super symmetry connecting 
vector particle scattering with scalar particle scattering, was suggested, and the combined use 
of these techniques was shown to be more advantageous by Parke and Taylor. The helicity 
amplitude method, which is useful in the calculation of cross-sections for processes with many 
final state particles, is reviewed in the article [S7| by Parke and Mangano. After these develop- 
ments, these techniques for massless particles were extended to the case of massive particles by 
Kleiss and Stirling and other authors f28fl . This useful extension makes us possible to analyze 
and calculate the processes involving massive particles more easily. 

These techniques enable us to obtain easily analytically compact expressions for various 
amplitudes, so they are extremely useful for the analysis of the process we are interested in. 
We use these techniques for top quark pair productions and their decay. In this chapter, we 
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present a brief introduction to the spinor helicity method following Ref. |27| . 

2.1 Massless Spinor States 

In this section, we derive and summarize expressions for the spinor products [^, |27|, [29] 
for massless fermions. The spinor helicity method utilizing the spinor products is extremely 
efficient for QCD with only massless quarks, and the massless case is a basis for the extension 
to the massive spinor case. We demonstrate an example using this technique in the next 
section. The extension of this technique to massive fermions is discussed in Sec. [2.3| . 
We will begin with massless helicity spinors and denote their chiral projection by, 



(l± 76 )b=F> = 0, 

(p±|(l± 75 ) = 0. 



3. Normalization and Completeness: 



4. Massless Spinor Relations: 



5. Chisholm Identity: 



(2.1) 



\p±) = u±(p) = v T (p) = ^u(p), 
(p±\ = u±(p) = v T (p). 

These massless spinors have the following properties. 

1. Dirac Equation: 

p\p±) = (p±\p = 0,(p = fp M ). (2.2) 

2. Chirality Conditions: 



(2.3) 



(p±| 7 "b±) = 2p», (2.4) 

\p+){p+\ + \P-){P-\ = P- (2-5) 



\p±)(p±\ = ^pp, (2.6) 

(p+\q+) = (p-\q-) = 0, (2.7) 

(P~\q+) = -(q-\p+), (2-8) 

(p-\p+) = (p+|p-) = 0. (2.9) 



(p±\ ^ \q±) [7m1 = 2[|g±> (P±\ + \PT) (g=F|]- (2.10) 



2.2. CIRCULAR POLARIZATIONS OF MASSLESS BOSONS 7 

Equations ( |2.7| )^( |2T9| ) show that most of spinor products vanish, and only a few spinor prod- 
ucts become non-zero. Therefore, it is convenient to introduce the following notation for spinor 
products: 

(pq) = (p-\q+) = ~(qp), ^ ^ 

[pq] = (p+ \q-) = ~[qp). 

These spinor products satisfy 

(pq)* = -[pq], (2.12) 

where all spinors are assumed to have positive energy. The squared spinor products become 

\(pq)\ 2 = \[pq}\ 2 = 2p-q. (2.13) 

The spinor product (p q) plays a fundamental role in this paper. Other properties of the spinor 
products are summarized in the Appendix D. 

2.2 Circular Polarizations of Massless Bosons 

We shall explain how to construct the polarization vector of massless gauge bosons by em- 
ploying the spinor techniques |^, |27], |2^]. The polarization vectors, £±(k), corresponding to 
states of definite helicity, satisfy, 

e±(k)-k = 0, (2.14) 

e£(*) = [4(*)]*, (2-15) 

e x {k) ■ e* x ,{k) = -6 xy , (2.16) 

where fc M is the momentum of massless gauge boson, and the subscripts A, A' denote the 
helicity states (±). 

In the four dimensions, the physical Hilbert space of massless vector is isomorphic to the 
physical Hilbert space of massless spinor, since they both lie in the one-dimensional represen- 
tations of SO (2). This isomorphism is realized through a linear transformation which relates 
the vector states to like-helicity fermion states: 

et{k) = N{k+\r\p+), r917 x 

4(fc) = N*(k-\-f\p-), { ' 

where s^(k) is the polarization vector of an incoming massless gauge boson, the momentum p 
is a light-like arbitrary reference momentum and assumed to be not parallel to the momentum 
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k. The factor N is a normalization constant. The factor N is determined from Eqn. fl2TT6|) : 

e-(k)-e + (k) = -2\N\ 2 p-k = -l. (2.18) 

Therefore, the polarization vectors of massless gauge bosons are given by 

*<*> = V^w 1 ^' (219) 

-i(t>(k,p) IU_\ ryfl |„_\ 

£ + W = - V2 W ■ (2 - 20) 

where e i ^ k,p ' is a phase factor which depends on the momentum k, and the reference momen- 
tum p. In this thesis, we set this phase factor unity and we obtain the polarization vector: 

t[) ~ T V2 ( P ±\kT) ' { } 



In Eqn. (^.21[), the reference momentum p can be chosen freely in a convenient way because 



of gauge invariance of the theory. The proof proceeds in the following way. By changing the 
reference momentum p in Eqn. ([2.21 ) to another reference momentum q, we have 



4(M = e^ + (3 ± (p,q,k)k^ (2.22) 

*-<***) - J^Ty (2-23) 

P + (p,q,k) = (f3 + y. (2.24) 

When the polarization vector £±(k) is replaced by its momentum k, the scattering amplitude 
KA should vanish: 

M = AV± -> MJP = 0. (2.25) 

Therefore the second term in the Eqn.( |2.22]) does not contribute to the helicity amplitude. 
This implies that we can choose a different reference momentum p for each polarization vector 
in the process, and furthermore different reference momentum for each gauge invariant subset 
of the full amplitude, without worrying about the relative phase. 

Now we show an application of these techniques to a real process [ 2"7| , fZPH . The process we 
consider as an example is the photon-photon scattering process: 

7(fci,Ai) +7(^2,A 2 ) -> g(pi,As) + g(p2,A 4 ), 

where k±, k 2 are the momenta of real photons, pi,p 2 are the momenta of final massless quarks, 
and Aj, (i = 1,2,3,4) denotes the helicity state (±). Two diagrams contribute to this process 
at the leading order, t channel (Fig. |2.1| (a)) and u channel fermion (Fig. [2.1| (b)) exchanges. 
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k u X { 



viM 



P2.^4 k 2 ,A 




2.^2 



Pl.^3 



P2M 



(a) (b) 

Figure 2.1: The tree level contributions to the 77 — » qq process 



The corresponding amplitudes are 

A*(A 3 ,A 4 ;Ai,A 2 ) = M(A 3 , A 4 ; Ai, A 2 ) +M U (X 3 , A 4 ; Ai, A 2 ) , (2.26) 

-M*(A 3 ,A 4 ;Ai,A 2 ) = -e 2 Q 2 q u(p 1 ,\ 3 )\e 1 (ki,\ 1 )- — — ^ e 2 (k 2 , A 2 )j^(p 2 , A 4 ), 

A1 u (A 3 ,A 4 ;Ai,A 2 ) = -e 2 (3;?u(pi, A 3 ) £ 2 (/c 2 ,A 2 ) r ^ — -2-ei(fci, Ai) v(p 2 , A 4 ), 

where Q q is the quark charge in the unit of electron charge e. It is easy to show that the 
helicity amplitudes with Ai = A 2 , and with A 3 = A 4 vanish. Therefore, only four helicity 
amplitudes corresponding to, 



M(A 3 ,A 4 ;Ai,A 2 ) = M(±,t;+,-),M(±,t;-,+) 



(2.27) 



should be considered. Here we consider the amplitude «M(+, —;—,+), and use the photon 
polarization vectors, 

1 



i(h, -1) 



V2(p 2 k 1 ) 
V2 



(h + \Y\p 2 +) i» 



(P2 h) 



e(k 2 ,+l) 



V2\pi h 

V2 



(\h-) (p 2 -\ + \p 2 +) (k 1+ \) 
(k 2 - |7 M bi-)7^ 



(2.28) 



bi ^2 



(|*2+><Pl + | + |Pl-><*2-|) 



(2.29) 



where the reference momentum of the polarization vector for the photon with momentum 
k\jk 2 has been chosen to be the momentum p 2 jp\. We obtain, 



M t (+, -;-,+) = e 2 Q 



1 ^2 (gi fe i) (P2 fci) [ fe i Pi] (fc2 P2) 
9 (Pi • &i) (P2 &i) bife] 
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2^2 [h Pi] (fa P2> 



(fa Pi)* (pi fc 2 ) 

-2e z Q 2 e 



2 n 2J<P(k 1 ,p 1 )(fa L P2}_ 



(Pi fa) ' 
, _ In-, n-i \ /rir, 

M u (+, -;-,+) = e 2 g 



2^2 (Pi Pi) (P2 fel) (&2~| (Pi - fc 2 ) |pi-) (p 2 P 2 ) 



9 (Pi • fo) (P2 fa) [Plfa] 

0, 



where e l ^ kl ' Pl > is the phase factor and given by e l ^ kl ' Pl > = [fa pi] / (fa pi)* ■ The amplitude 
Ai u gives zero, because (P2P2) = (piPi) = 0. Thus we immediately arrive at the final result, 



i^fa-P2 

/q fa -P2 



|^(+,-;-,+)| 2 = 4e 4 Q 4 -^ 2 , (2.30) 



By using the invariance under the CP transformation, the remaining non-zero helicity ampli- 
tudes are, 

\M(+,-;-,+)\ 2 = |-M(- +;+,-)| 2 = 4e 4 Q 4 ^, (2.31) 

H fa ■ Pi 

|M(+,-;+,-)| 2 = |^(-,+;-,+)| 2 = 4e 4 Q 4 ^ 2 -. (2.32) 



This example shows the efficiency of the spinor helicity method to calculate the cross-section. 
The squared amplitude for 77 — > qq process, averaged over initial spins and summed over final 
and colors (N c = 3) reads, 



,t 2 +M 2 



\M\' = 2N c e^ q (—- r ). (2.33) 

Here we use the Mandelstam valuables, t — (fa — pi) 2 and u — (fa — P2) 2 ■ 
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2.3 Massive Spinor States 

We consider, in this section, massive spinor states and show that the massive spinor states 



also can be described in terms of massless spinor states [p3i [27 |. At first, we summarize the 
properties of the massive spinors u(p, s), v(p, s) with momentum p and polarization s. 

1. Dirac Equation: 

(ft — m)u(p, s) = 0, u(p,s)(p — m) = 0, (2.34) 

(p + m)v(p, s) = 0, v(p, s)(p + m) = 0. (2.35) 



2. Normalization, Orthogonality and Completeness: 

u(p, si) u(p, s2) = 2m5 sltS2 , 
v(p, si) v(p, s2) = -2mS sltS2 , 

u(p, si) v(p, s2) = 0, 
v(p, si) u(p, s2) = 0, 



(2.36) 
(2.37) 



E s u(p,s)u{p,s) = p + m, 
J2s v (Pi s ) v(p, s) = p — m. 

The spin vector s M of massive fermion has the form in the fermion's rest frame: 

a" = (0, s), (2.39) 

with s ■ s = 1. The spin vector, when the fermion is moving with momentum p, become 

y m m(E + m) J 



where p = (E, p) and spin vector satisfies p-s = 0,s-s = — 1. The projection operators PU 
are given by 

1 ' " " u(p), (2.41) 

v(p). (2.42) 



(2.43) 
(2.44) 
(2.45) 
(2.46) 



u(p, s)u(p, s) 


- (p + m) 


v(p,s)v(p,s) 


= (p m) 2 


Therefore, the polarized fermion states are projected as, 


«(p, T) 


l + 7sS , . 
2 u\P), 


u(p, i) 


1-75S , s 

2 w(PJ, 


v(p, t) 


1+75S , V 

2 u ^' 


v(p, i) 


1-755 , X 

"IP). 
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where spin states, f/l refers to fermions with spin in +s/—s direction in their rest frame. 
Now, we rewrite Eqns. ( fr%3| ) ~ ( 2~4l]| ) in terms of massless helicity spinor states and show how 



to describe the massive spinor states in terms of massless helicity states. Let us consider, for 
example, the state u(p, |) and rewrite it using properties satisfied by gamma matrices, 

/ t x 1 + 75 1 + 75S 1 - 75 1 + 75S , . 

u(pA) = — - — u{p) + — — u(p), 

I + 75I + S I-75I-S 

= — g y~ u (p) + — 2 y~ u ^' ^ ' 

In Eqn. (|2.47|) , the quantities [(1 ± s)/2] u(p) can be identified with massless spinors having 



momentum (p ± ms)/2 due to the following reasons [31 



1. Dirac Equation: At first, note that the momentum p^ (i = 1,2) is light-like, 

Pi = P -±P, (Pi = 0), (2.48) 

p — ms , . . , . 

P2 = P -^— , {vl = 0), (2.49) 

and simple relations between p if p and s, 

P = P1+P2, (2.50) 

ms = P\ — Pi- (2.51) 



Next it is easy to show, 



Pi^—u{p) = P2^—u{p) = 0. 



Therefore, the quantities [(l±s)/2] u(p) are proportional to the massless spinors Q(pi), 
Q{p 2 ) respectively. 

Q(p x ) = A { -^-u(p), (2.52) 

Q(p 2 ) = B { -^-u(p), (2.53) 

(2.54) 

where factors A, B are normalization constant. 

2. Projection Property: The normalization factors A and B will be fixed by constructing 
the projection operators, which read 

Q(pi)Q(pi) = \A\ 2 p u 
Q(P2)Q(P2) = \B\ 2 p 2 , 
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Therefore, we conclude that [(1 ± s)/2] u(p) are the massless spinors Q(pi) except for phase 
factors. Above consideration tells us 



u(p, T) 



L +75^„ , , ^(p 1)P2 ) 1 - 75 



Q(pi) + e l 



Q(P2 



= \pi+) + e 



if t>(Pl,P2) 



|P2-> 



(2.55) 



Here, we have chosen the overall phase factor such that there is no phase in front of \pi+). 
Therefore, 4 > {pi,P2) represents the relative phase between the states \pi+) and \p%— ), This 
phase factor can be determined by the normalization condition for u(p, |) in Eqn.( 2~36l ). It is 
not difficult to obtain, 

U, _ (P2 Pi) 



III 



where we have used the relation: 



(P2 Pi) I 2 = 2pi -p 2 = m 2 . 



(2.56) 



(2.57) 



In this way, all massive spinor states can be written in terms of massless chiral spinors as (for 
the anti-fermions, we use the superscript "bar" for their momenta), 



u(p,t 
u(p,i 

«(p,t 

u(p,i 

u(p,t 

v(p,i 

t?(p,t 

v(p,l 



\P1+) 


+ lp 2 ) {P2 P1) 

m 


\P1-) 


+ b+> hPl1 , 

771 


(Pl+I 


+ ¥l P3] (■» 1 

771 


(Pl-I 

IP1"> 


+ (Pl ^ (P 2 + | 

m 

-ift + ) tePl1 , 


IP1+) 


- la ) (fc ft > 


(Pl-I 




(Pl+I 


[P1P2},- | 

(P2 1 



III 



The light like momenta pi, p%, p\ and pi are defined by, 



Pi 

Vi 



p + ms p 

2 ' 

p + msp 



P2 



P2 



p — ms p 



p — msp 



(2.58) 
(2.59) 
(2.60) 
(2.61) 
(2.62) 
(2.63) 
(2.64) 
(2.65) 

(2.66) 
(2.67) 
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For the convenience of later calculations, let us write down pi explicitly in the rest and 
moving frame of the massive fermions. In the rest frame, the four momentum of the massive 
fermion and the spin vector take 



P 



So, the light like momenta pi, P2 are 

m 



(m, 0), 

(0, «), (M 



I)- 



Pi 



1, w 



m 



Pi 



u> 



(2.68) 
(2.69) 



(2.70) 



2 v 7 /7 " 2 

These equations show that the direction of the particle's spin is in the same direction of the 
spatial part of p\ or in the opposite direction of the spatial part of pi- For the moving fermion, 
p = (E, Vn), (\n\ 2 = 1), pi are obtained by boosting Eqn. (|2.69|) ~ (|2.70| ). The spin vector and 
decomposed light-like momenta pi, P2 become 

V 2 (n-u)) n 



V- -,w4 

m 



m(m + E) 



Pi 
Pi 



-[E + V (n-u),V n + m u + —± ^~ 

2 V (m + E) 

2 V v ;> (m + E) 



(2.71) 
(2.72) 
(2.73) 



Here we consider the familiar helicity state, that is, the spin of massive particle is projected 
along its direction of motion, 

(2.74) 
(2.75) 



n - 


= w, 




Pi = 


= \(e + v, 


L v2 1 

V + m+- — 

(m + E) 




- U E -^ 


L v 2 1 


Pi 


' ''*• "" / 7-1 \ 




2 V 


(m + E) 



n 



n 



In the large momentum limit (which is equivalent to the massless limit), p\ 
Then, the spinors become the chiral eigenstates: 



P, Pi 



u(p, R) - 


-» |Pi+), 


u(p,L) - 


- IPi") 


v(p,R) - 


- \Pi-), 


u(p,L) - 


* |pi+) 


u(p, R) - 


■* (Pi~l> 


u(p,L) - 


-> (pi+| 


v(p,R) - 


■* (Pi+I, 


v(p,L) - 


* (Pi~l 



(2.76) 
= 0. 

(2.77) 
(2.78) 
(2.79) 
(2.80) 
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where we have labeled helicity states by R (L) instead of j (J,). This demonstrates the relation 
between the right-handed (left-handed) helicity state, u(p,R) (u(p,L)), in the massless limit 
and the right-handed (left-handed) chiral state. 
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2.4 Circular Polarizations of Massive Bosons 

The polarization vectors for massive gauge boson |28|, ^, |32] can also be described in terms 
of massless spinors. Transversely polarized states satisfy Eqns. ( p7T4[) ~ ( |2.16|) in Sec. |2.2| . For 



massive gauge boson, we must also take into account the longitudinally polarized state. Let 
e± denote transversely polarized states and e$ is the longitudinally one. The properties of 
polarization vectors for a massive gauge boson with momentum k and mass m are 

e x (k)-k = 0, (2.81) 

e x (k) ■ e* x ,(k) = -5 xy , (2.82) 

_ k^k u 

E^^rw = -<r+^-, (2.83) 

where A and A' denote +, — and 0. 

Let s be the spin vector associated with the massive vector boson, which satisfies, 

k-s = 0, (k 2 = m 2 ). (2.84) 

Here three polarization vectors e x are defined such that in the vector boson rest frame, the 
boson has the spin projection A = (+, 0, — ) with respect to the spatial part of the spin vector 
s. Three polarization vectors e x (k) are given in terms of two light-like vectors, 

, ,. k^ — ms^ 

K = — 2 — ' (2 ' 85) 

K = k(1+ 2 mS ". (2-86) 

Then, we get polarization vectors for the incoming gauge boson: 

(k 1 ±\Y\k 2 ±) 
el = ^ » ( 2 - 87 ) 

(h+\r\h+)-(k2+\r\k2+) k - 14 r9SS , 

e - - . {2.XK) 

2m m 

In this expressions, the helicity basis is given by choosing the spatial part of the spin vector s 
to be in the same direction as the spatial part of the momentum vector k. 



Chapter 3 

Spin Correlations at Leading Order 



In 1994, the top quark was discovered by the CDF and DO collaborations |4|, |!|. The measured 
top quark mass is approximately 175 GeV, which is nearly twice the mass scale of electro- 
weak symmetry breaking. Thus the top quark with large mass brings a good opportunity 
to understand electro-weak symmetry breaking and to reveal physics beyond the Standard 
Model. The top quark decays electroweakly before hadronizing because its width is much 
greater than the hadronization timescale set by Aqcd [0- Therefore, there are significant 
angular correlations between the decay products of the top quark and the spin of the top 
quark. These angular correlations depend sensitively on the top quark couplings to the Z° 
boson and photon, and to the W boson and b quark. These angular distributions are considered 
as good information to constrain the couplings of top quarks in the Standard Model. Most 



works Jy], [12] are performed by using top quark spin decomposed in the helicity basis. In 
the case of ultra-relativistic particles, this decomposition is appropriate. However this basis is 
not necessarily optimal to investigate the spin correlations since the produced top quarks are 
non-relativistic at realistic linear colliders: e.g. /3 = v/c ~ 0.5 for a 400 GeV collider. As a 



matter of fact, Mahlon and Parke [15] have shown that this decomposition in helicity basis is 
far from the optimal decomposition for top quark productions at hadron colliders. Parke and 
Shadmi |TE| have extended these analyses to lepton colliders. 

In this chapter, we discuss the spin-spin correlations for the top quark pair production at 
polarized e + e~ linear colliders at the leading order in perturbation theory. We introduce the 
"generic spin basis" and investigate which decomposition of top quark spin is optimal to study 
spin correlations for moderate energies. The spin-angular correlations for the top (anti-top) 
quark decay are also discussed. 
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3.1 Spin-Spin Correlations for Top Pair production 

In this section, we derive the leading order spin dependent differential cross-sections for top 
quark pair productions in a generic spin basis. We find the most optimal spin decomposition 
for top and anti-top quark to study spin correlations. The one-loop analyses will be given in 
the next chapter. 




Figure 3.1: The tree level contributions to the e e + — > tt process 

The Feynman graph for the e + e _ — > ti process at the leading order in perturbation theory 
is given in Fig.[TT| We calculate the amplitude using the spinor helicity method. The momenta 
of electron, positron, top quark and anti-top quark are labeled by q, q, t and t respectively, 
and s t (sf) represents the spin vector of top (anti-top) quark. As explained in the previous 
chapter, we decompose the top (anti-top) quark momentum t (t) into a sum of two massless 
momenta t 1 ,t 2 (ii,i 2 ), 

t = h + t 2 , (i = h + t 2 ), (3.1) 

ms t = h-t 2 , (msf = ti — t 2 ). (3.2) 

The amplitudes for e^e^ (e^ej) scattering processes are obtained to be, 

M(eie%-> t st t 8l ) 
= (q-\ 1l !<?-) u(t, stfrf, \JLLlh + IlrIr] v(t, Si) , (3.3) 

M{e R et -* t St t s -) 

Ano/. - 

= {q+\ 1r \q+) u(t, sAjn [JrlIl + IrrIr] v(t, Sf) , (3.4) 

where e^/e^ (e^/ej) means the electron/positron with left-handed (right-handed) helicity 
state. We have neglected the electron (positron) mass. The quantity a is the fine structure 
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constant a = e 2 /(47r), and s is the total energy, 7^, L = 7^7_r/l, 1r/l = (1 ± 7s)/2. The 
quantities /jj's are the sum of photon and Z° boson couplings to fermions (electron and top 
quark) including the propagator and are given by 

1 s 



fu = -Qt + QiQ 



eVt ,2 



sin 2 9 W s- Mf + iM z T z 

Here Mz is the Z boson mass, Y z is the width of Z boson, 9w is the Weinberg angle and 
I, J G (L, R). The parameter Q t = 2/3 is the electric charge of the top quark in the unit of 
electron charge e. The electron couplings to the Z boson are given by 



Qi 



2 sin 2 9w — 1 ^p sin 2 9 W 



Q 



2 COS 9w ' COS 9\y 

The top quark couplings to the Z boson are given by 

R 2 sin 2 9w 



Qt 



3-4 sin 2 9 W 



6 cos 9w 
The squared amplitude for e^e^ scattering reads 



Qi 



3 cos 9 



w 



l>t(e £ e5->t T t T ; 



41 T) 2 









\M(e L e+ -* t^; 

. 4 (^Y 



|/ LL | 2 (2g-t 1 )(2g-t 2 ) + |/ Li? | 2 (2g-t 2 )(2g-t 1 
1 



+ 



77?/ 



;{fLLft R T 1 c[y L qt 1 t2qt2ti} + c.c} 



\f LL \ 2 (2q ■ h)(2q ■ t x ) + \f LR \ 2 (2q ■ t 2 ){2q ■ i 2 

+ — J {/LL/ifl1t[7L9*l*29*l*2] + CC} 
1 1 L 



\f LL \ 2 (2q ■ t 2 )(2q ■ t 2 ) + \f LR \ 2 {2q ■ t 1 )(2 9 • t^ 
+ — {/LL/^Tr[7 L g£ 2 tig£ 2 ti] + c.c] 

1 1 L 



\f LL \ 2 (2q ■ t 2 )(2q ■ t x ) + \f LR \ 2 (2q ■ t 1 )(2q ■ t 2 
+ —{fLLflRTrilLqhtiqtxh] + c.c} 

1 1 L 



(3.5) 



(3.6) 



(3.7) 



(3.* 



(3.9) 



(3.10) 



where all momenta, p, under the "Tr" operator are understood to be p. The squared amplitude 
for e^ej scattering are given by interchanging L, R as well as f, J, in Eqn.( |3.7p ^( |3.10| ). 
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Before discussing the spin basis, let us fix the coordinate system. In the center of mass 
(CM) frame, we take the production plane to be x — z plane and set the direction of the motion 
of the top quark to be +z direction (see Fig. 




+x 



Figure 3.2: e e + — > tt process in the center of mass frame. 



Now, the generic spin basis we consider here is based on the following facts: First, CP is 
conserved in the Standard Model. Hence the spins of top and ant-top quark can be defined 
to be back-to-back. Second, there is no transverse polarization of the top (anti-top) quark 
at the tree level in the Standard Model. Here the "transverse" means the direction normal 



to the production plane [Q. Hence, the spin vectors are on the production plane and they 
are parametrized by one parameter. We define the generic basis in the following way; For 
top quark, we got to its rest frame and its spin is decomposed along the direction s t in the 
rest frame of the top quark which makes an angle £ with the anti-top quark momentum in 
the clockwise direction. Similarly, the anti-top quark spin states are defined in the anti-top 
rest frame along the direction sj having the same angle £ from the direction of the top quark 
momentum (see Fig. |3.3j ). 
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Figure 3.3: The generic spin basis for the top (anti-top) quark in its rest frame. s t (sj) is the 
top (anti-top) spin axis. 



In the top (anti-top) quark rest frame, the top (anti-top) quark momentum t RF , (t RF ) and 
the spin vector become s^ F (sf F ) as, 



e RF = m(l, 0,0,0), F RF = m(l, 0,0,0), 

[sf F Y = (0,-sin£,0,-cosO, {s? F y = (0, sin£, 0, cosO- 



(3.11) 



Now let us decompose the top (anti-top) quark momentum t RF (t RF ) in terms of massless 
momenta tf F and tf F (if F and tf F ), 



t RF — ti + 1 2 5 Tns-f- — ti — £ 2 , 



'RF 



j-RF _j_ ±RF qRF 4.RF ±RF 



(3.12) 
(3.13) 



From these relations, the four massless momenta are expressed as 



m 



(tf F r = -(l,-sin£,0,-cosO, 



m 



(t* F Y = f(l, sine, 0,oobO, 

(tf F r 



m 

m 



l,sin£,0,cos£), 

1, -sin£,0, -cos£)- 



(3.14) 
(3.15) 
(3.16) 
(3.17) 



To get back to the CM frame, we boost back all quantities in the +z (—z) direction with the 



top (anti-top) quark speed (3 = J\ — Am 2 / s. The momenta of particles in the CM frame are 
given by, 



t" = m 7 (l,0,0,/3), i" = m 7 (l, 0,0,-/3), 

q^ = 7717(1, sin 9, 0, cos 9) , (f = 7717(1, — sin 9, 0, — cos 1 



(3.18) 
(3.19) 
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m 



-77(7(1 - /3cos£), - sin£, 0,7(/3 - cos£)) 



/?? 



-(7(l + /?cosO,-sine,0, 7 (/5 + cosO) 



/» 



— (7(1-/3 cos £),- sin ^,0, -7(/3- cos £)) , 



/?? 



-(7(1 + /?cos£),- sine, 0,- 7 (/? + cosO) 



(3.20) 
(3.21) 
(3.22) 
(3.23) 



where 7 = l/\/l — (3 2 and 6 1 is the scattering angle of the top quark with respect to the electron 
in the CM frame [see Fig. |3.2[ |. 

Let us write down the squared amplitudes using kinematical variables in the CM frame. 
Since each term in Eqn.(|3.7]) reads as 



{2q-t 1 ){2q-t 2 ) 



= m 4 7 4 [(l - /3cos^) 2 - {(cos 6 - (3) cos£ -\]l- /3 2 sinflsin£} 2 ] 

{2q-t 1 ){2q-t 2 ) 

= m 4 7 4 [(l + [3cos6) 2 -{((3 + cos 9) cos«£ + y/l - /? 2 sin#sin£} 2 ] 

Tr[7 Lqhhqhh] = Tr [7^9*1*2^2*2] 



= m 6 ^! - [3 2 - {\/l - (3 2 cos6 cos£ + sin6 sin£} 2 } , 

the squared amplitude becomes 

\M(e£e% -v t T t T )| 2 = l67r 2 a 2 \A LR cos£-B LR sm£\ 2 , 

where the quantities A LR and B LR are defined by 

Alr = [(f L L + fLR)^/l-(3 2 sme}/2, 
B L r = [f L L(cos8 + [3) + f LR (cos6-[3)}/2. 

Similarly, other squared amplitudes are 



\M(e T e 



L^R 



tit 



IHj 



\M(e L e+^t ] t i ort ; t T 



with 



D 



LR 



= 1 6-K 2 a 2 \A LR cos £ - B LR sin £| 2 , 
)| 2 = leTrVlA^sine + E^cos^iD^I 2 

/ LL (1 + (3 cos 6) + / M (l - $ cos 0)1/2 . 



The differential cross-section is given by 

da 1 

dcosO 327r s 



E\M? 



3.24) 
3.25) 
3.26) 

3.27) 

3.28) 
3.29) 

3.30) 
3.31) 

(3.32) 



N c 
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where N c is the color degree of freedom, N c = 3 and s is the total energy. Thus, we arrive at 
the final expression of the polarized cross-sections for the top quark production, 



da 



da 



dcos9 



(e L ek 



y e L e R 



Ut 



da 



1 L V 



dcos9 

3ira 2 



^ e L e R 



td 



iH> 



hh or *i*t) 



2s 

37T« 2 



/3\A LR cost, - B LR sin £\ , 
(3\A LR sm£ + B LR cos£ ± D LR \ 



(3.33) 



Now, consider numerical behavior of the parameter fjj. The production threshold for 
top quarks is far above the Z boson mass. The contribution of the Z boson width to fjj is 
suppressed by a factor 

(, - M^% Iz T zf ~ (°' ° 15 » ' 
when y/s = 400 GeV, Mz = 91 GeV, Y z = 2.49 GeV. Therefore, in our discussions bellow, 
we neglect the Z boson width. We use the values for the parameters of the standard model 



in Table |3T1 . 



m Mz a sin 2 9w 


175 [GeV] 91.187 [GeV] 1/128 0.2315 



Table 3.1: Input parameters of the Standard Model. 



We show the j3 dependence of fi/s in Fig. 3.4, 
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Figure 3.4: Top quark speed, /3, dependence of f LL , f LRl f RL and f RR . 



For the top quark production, fjj have a weak f3 dependence and are given by 

- 1.23 < f LL < -1.19 , -0.882 < f LR < -0.867 , 
-0.433 < f RR < -0.417 , -0.217 < f RL < -0.184. 

When fermions do not couple to Z boson and couple only to the photon, all parameters fi/s 
go to —2/3. The difference, fu — (—2/3), shows the contribution from the Z boson, and we 
can investigate their effects in all energy region. 

Since we are interested in maximizing the spin correlations of top quark pairs, we vary the 
spin angle, £, to find an appropriate spin basis. Following the articles |L5], [I6| , we introduce 
"helicity" , "beamline" and "off-diagonal" bases and derive the polarized cross-sections in these 
bases from Eqn.( |3.33|) . 

1. Helicity basis 

The helicity basis, the most familiar spin basis, is given by 

cos£ = ±1, (3.34) 
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for which the top quark spin is defined along its direction of motion. It is noted that the 
helicity basis can not be used in the threshold energy region, because the top (anti-top) 
quark momentum is undefined in this region. Substituting Eqn.( |3"mD to the general 
polarized cross-section in Eqn. (|3.33|) , we can obtain the well known polarized cross- 
section: 



da 



da 



dcosO 



dcosO 



v e L e R 



■£lZr 



trt 



da 



L^L) 



e,e, 



dcos9^ L R 
3ira 2 



tpi 



R^R) 



t R t L or t L t R ) 



8s 
37ra 2 



8s 



(3\f LL + f LR \ 2 (l-f3 2 )sm 2 6, 
P\f LL (lT (3) + f LR (l± f3)\ 2 (lTcosi 



(3.35) 



Here ti (t R ) state corresponds to t-j- (tj_) state with cos£ 



2. Beamline basis 

In the rest frame of top quark, there are three natural choices for the direction of the 
top quark spin; direction of the electron, the positron or the anti-top quark momentum. 
In the threshold region, the anti-top quark momentum direction, which is undefined, is 
excluded leaving only the electron or positron momentum directions. The situation is the 
same for the anti-top quark spin. The natural choice (see the discussion at the end of this 
section)is that the top (anti-top) quark spin is in the direction of the positron (electron) 
momentum in its rest frame. This spin basis is called "beamline basis " |T5|, [16], Q , The 
spin angle £ is determined by 

cos 6 + p 



C ° Se 1 + 0COS0" 

The polarized cross-sections in the beamline basis are 



da 



dcos9 



K e L e R 



hh, 



da 



^ e L e R 



hh) 

32/1 o2\„-„2 



dcos9 

(3na 2 \ (3 2 (l-(3 2 )sm 2 9 

-P \JLR\ 



da 



[e T e 



dcos6 K " L R 

do- i _ + 

~, n\ e L e R 

dcosv 



t\t 



i L h 



hh) 



2s 
2s 
2s 



(3)\fLR 

p 



(l + (3cos6) 2 
/? 4 sin 4 fl 



v l + /3cos#) 2 ' 
f LL {l + 13 cos 6) + f LR - 



p? 



1 + P COS I 



(3.36) 



(3.37) 



3. Off-diagonal basis 
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There exists the third possibility called "off-diagonal" basis which makes the contri- 
butions from the like-spin configuration vanish ||16|| . The spin angle, £, is determined 

by 

** B LR f LL (cos9 + (3) + f LR (cose-P)- { - ' 

In this basis, we obtain the polarized cross-sections, 

da , _ , - do , 

eZ4 ~» *t*l or *l*t) 

^](^L + 5LT%)' (3.39) 

/L£(l + /3cos0) + / Lfl (l-/3cos0) 



dcos9 y n n ' dcosfl 

da 



dcos^ 

/37ra 2 \ / rr^ — — „ \2 



2s 
'3?ra 2 



S.s 



/? 



T V[/ll(1 + ^os6) + f LR (l - (3cos6)Y - 4f LL f LR (3 2 sm 2 9 
It is to be noted that we can consider another "off-diagonal " basis, 

. t_A R L (fRR + fRL)VT=W Sing 

^^ B RL f RR (cos9 + (3) + f RL (cos9 - (3) ' l ' ' 

which makes the contributions from the like-spin configuration vanish in the e R e~£ scatter- 
ing process. So, we can consider two "off-diagonal " bases for top quark pair production 
because f LR /f LL ^ Irl/Irr- However, the difference between f LR / f LL and Jrl/Jrr 
is not large numerically (see Fig. p.4p . Therefore we call the basis defined by Eqn.(|3"l 
as "off-diagonal" basis in the following discussions. 



Let us consider the difference between above three bases. In Fig.|3.5|, we plot the dependence 



of the spin angle £ on the scattering angle 9 for the helicity, beamline and off-diagonal bases 
at ft = (threshold), (3 = 0.5 (y/s ~ 400 GeV) and (3 = 1 (ultra high energy). 
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(3=0 



P=0.5 



P=1 




-1.0 -0.5 0.0 0.5 -1.0 -0.5 0.0 0.5 -1.0 -0.5 0.0 0.5 1.0 
cose cose cos9 

Figure 3.5: The dependence of the spin angle, £ on the scattering angle 9 in helicity, beamline 
and off-diagonal bases at (3 = 0,0.5 and f. It is noted that the off-diagonal basis at (3 — 1 
changes from cos£ = -1 to cos£ = +1 at cos6» = -(f LL - f LR ) /(f LL + Ilk)- 



When (3 is zero, the helicity basis can not be defined. However, the beamline and off-diagonal 
bases are defined also at threshold and they are given by the same equation, £ = 9. This 
fact says that it is correct that the top (anti-top) quark spin decomposes into the direction of 
positron (electron) momentum in the threshold region. When (3 = 0.5, we can see apparently 
the difference among these three bases. In particular, for the central region in which cos 9 ~ 0, 
there are remarkable differences between the three bases. However, the difference between 
the beamline basis and off-diagonal basis are small for all scattering angle compared to the 
helicity basis. The beamline and off-diagonal bases become close to the helicity basis with 
cos£ = ±1 near cos 9 — ±1. When (3 = 1, the beamline basis coincides with the helicity basis 
with cos£ = +1. While, the off-diagonal basis is equal to the helicity basis with cos£ = —1 
from cos 9 = — I to cos 9 = — (j'll — /l_r)/(/ll + Ilr), and equal to the helicity basis with 
cos£ = +1 form cos 9 = —(fLL — Ilk)/ '(III + Ilk) to cos 9 = +1. The off-diagonal basis 
corresponds to the different helicity bases at cos# ~ 1 and cos# ~ — I. 

Now we will show the polarized cross-sections for the helicity, the beamline and the off- 
diagonal bases in Fig.^76] at y's = 400 GeV. 
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Figure 3.6: The polarized differential cross-sections in the helicity, beamline and off-diagonal 
bases at y/s = 400 GeV for the e~ L e + (e~ R e + ) — > ti process. We plot the tiXh (LL), t R in (RR), 
tail (RL), tnii (RL) productions for the helicity basis and t^t^ (UD), ijjj (DU), t-ft^ (UU) 
and i^fj (DD) productions for the beamline and the off-diagonal bases. Note that the sub- 
dominant and sub-sub-dominant configurations for both e~[e + and e^e + process are amplified 
by the factor of 10 in the beamline basis and by the factor of 100 in the off-diagonal basis. 



e L e+ 


Helicity 


Beamline 


Off-Diagonal 


Dominant Frac. 


0.5305 (LR) 


0.97846 (UD) 


0.99876 (UD) 


Sub-Dominant Frac. 


0.2033 (RL) 


0.00964 (UU,DD) 


0.00124 (DU) 


e «e + 


Helicity 


Beamline 


Off-Diagonal 


Dominant Frac. 


0.5804 (RL) 


0.98979 (DU) 


0.99744 (DU) 


Sub-Dominant Frac. 


0.1608 (LR) 


0.00457 (UU,DD) 


0.00089 (UU,DD) 



Table 3.2: The fraction of e L , R e + cross-sections for the dominant and the sub-dominant spin 
configurations in the helicity, beamline and off-diagonal bases at yfs = 400 GeV. 
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The dominant component for the helicity basis is tijtn (LR) for e~ L ~e + scattering and tj$L (RL) 
for e~ R e + scattering at y/s = 400 GeV. This dominant component occupies 53% of the total 
cross-section for e~ L ~e + scattering and 58% of the total cross-section for e^e + as shown in Table 



3721 Other components make up more than 40% of the total cross-section. Thus, the helicity 
basis does not show strong spin correlations for the polarized cross-sections at this energy. 
In contrast, in the beamline and the off-diagonal bases, only one component dominates the 
total cross-section. The dominant component for the beamline and the off-diagonal bases is 
t|t| (UD) and makes up more than 97% of the total cross-section for e^e + scattering. For 
e] i e + scattering, the dominant contribution in the beamline and off-diagonal bases is t±ii (DU) 
component and make up more than 98% of the total cross- sect ion. Other components give 
almost zero contribution. Apparently the top (anti-top) quark spin is strongly correlated with 
positron (electron) spin in the beamline and off-diagonal bases. 

Here we focus on the beamline and off-diagonal bases and try to explain (l)why the top 
(anti-top) quark spin is associated with the positron (electron) spin and (2) why the polarized 
cross-section in the beamline and off-diagonal bases show similar behaviors. To answer these 
questions, recall the expression Eqn.(|3.3|) for ej j e + —> t St t Sr 

M{ele^^t st t Sl ) 

Afcct - 

= (q-\ r rL\q-)u(t,S t )7n[fLLjL + fLRjR}v(t,Si). (3.41) 

Using Fierz transformations in Eqn. (|3.41| ), we obtain, 



M{e L e + R ^t St t Si ) = f LL 



u{t,s t ) \q+) (q+\v(t,St) 



+-j£ u(t, s t ) \q-) (q+\ v(t, Si) 

ILL 



(3.42) 
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Figure 3.7: Top quark speed, (3, dependence of ratios, JlrIJll and JrlIJrr- 



Figure ( |3.7| ) shows the top quark speed (3 dependence of the ratio //j's. For the top quark 
production, the //j's have a weak energy dependence; Numerically the parameters fu fall in 
the range, 



0.34 <i^- < 0.37, 



i 



LL 



0.21 < ^- < 0.25. 



./' 



RR 



(3.43) 
(3.44) 



Because the ratio |/l.r//ll| is less than 1, the first term of Eqn. (|3.42j ) is dominant. This 
implies that the top (ant-top) quark spin is associated with positron (electron) spin. 

We now come to the question (2). Now remember the definition of the spin angle £ in 
the beamline and the off-diagonal bases in Eqns.([T36) and Q3.38| ). The spin angle £ for the 
beamline basis is obtained from Eqn.( |3.38|) with Jlr = 0. Then we focus on the factor 
/lr to answer the question (2). Since the factor |/l_r| is smaller than |/z,l|, let us take an 
approximation of fm = in Eqns. ( pT37| ) and ( p.39|) . The polarized cross-sections in the 
beamline and off-diagonal bases become the same expressions as 



da 



dcos9 

da 



(e L e 



R 



dcos6 [eLeR 



t T t T/ 



hh) 



da 



dcos9 

0. 



y e L e R 



tit 



J, iy 



(3.45) 
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do , , - s /3na 2 ' x 



2/1 , /3„„„/3\2 



eltR^hh) = \-z-P\\fLL\\l + Pcm 



dcos6 K ^ n ' +y V 2s 

Thus the difference between the polarized cross-sections in the beamline and the off-diagonal 
bases is C(|/l_r//ll|)- The subdominant contributions, (UU, DD, DU), come from the 
0(\fLR./ Ill\) corrections. This is the reason why polarized cross-sections in the beamline 
basis behaves like that in the off-diagonal basis. 

Finally we discuss the top quark speed j3 dependence of the polarized cross-sections in 
the helicity, beamline and off-diagonal bases. We define the fraction a(ej i e + — > t St t Sl )/<jT of 
top pair productions for each spin configurations, where <jt is the total cross-section for the 
process, e~ L e + — > ti and the suffix s t (s?) represents the spin of top (anti-top) quark. We 
present the fraction in the helicity, beamline and off-diagonal bases in Fig. |3.8| . 
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Figure 3.8: The fraction of top quark pairs in the helicity, beamline and off-diagonal bases, as 
a function of the top quark speed, (3. 



In the helicity basis, as the CM energy get larger, the spin correlation becomes stronger for 
e~[e + and e^e + scattering. The polarized cross-sections in the beamline basis coincide with 
those in the helicity basis in the limit /3 to 1. This is in agreement with the fact that the helicity 
and beamline bases become identical in the limit of (3 to 1 (see the discussion on Fig. |3.5|) . In 
the beamline and off-diagonal bases, the UD (t^i) channel dominates for e~[e + process for all 
values (3 as shown in Fig. |3.8| . (The DU (tjtj-) channel for e^e + gives the dominant fraction 
in the beamline and off-diagonal bases.) However, there is a slight difference between in the 
beamline basis and in the off-diagonal basis near (3=1. This is due to the fact that the 
off-diagonal basis corresponds to the helicity basis with either cos£ = +1 or with cos£ = — 1, 
depending on the scattering angle, while the beamline basis coincides the helicity basis with 
cos£ = 1 in the limit (3 to 1. 
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3.2 Top Quark Decays at Leading Order 

In the previous section, we have shown that the top quark pairs are produced in an unique 
spin configurations at polarized e + e~ linear colliders in a appropriate basis. We discuss in this 
section the spin-angular correlations for the top quark decay and measurement of the top quark 
spin. Because of the heavy mass of top quark, the width of top quark amounts to nearly 1.57 



GeV ||10| . In consequence, the top quark decays rapidly through the electro- weak interaction 
before the hadronization (which is governed by the scale Kqcd) take place. Therefore, we 
can obtain information of the top quark spin by measuring its decay products without being 



suffered from the complicated hadronization effects ||35|| . We show the electro- weak decay 
products of polarized top quark are strongly correlated to its spin axis at the leading order in 
perturbation theory of the Standard Model. 




v, u 




V, U 



Figure 3.9: The hadronic and semi-leptonic decay process of top (anti-top) quark 



The top (anti-top) quarks decay dominantly, 



t -> W + b, t 



W~b. 



(3.46) 



For the W + (W ) decay, there are hadronic decay and leptonic decay modes (see Fig. |3T9j ); 

(3.47) 



W 
W 



ud or cs, W~ 



du or sc, 



Now let us calculate the decay distribution for the top quark decay. (The anti-top quark 
decay can be calculated similarly.) Momenta of the particles are assigned by its symbol, and 
we decompose the top momentum t into a sum of two massless momenta t% and £2, 



t = U+U. 
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The spin vector of top quark is related to the momenta t\, £2 by the following relation, 

St = 7^(tl-h). 

The amplitudes for t^/t^ — > blui are given by, 

Mto^Hk) = -V, Z ,h{ \^ U i F ^. (3-48) 

V r ; y (t - bf - M& + iM w Y w m ' v ; 

*(»,-*,) = - 9 \ t Z* (b £i% v ■ < 3 - 49 > 

The mass and the width of W boson are M w and T w , and SU(2) coupling g is related to 
the electro-magnetic coupling through g = e/ sin9 w . V t b is the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix element. For the hadronic decay of the top quark, one should replace I with d 
and v with u in Eqns. ( fT%£| ) and ( 3.49| ). The squared matrix elements for the top (anti-top) 
quark decay, t — > blu h are given by 

|A ^ - 6 ^' 2 = ett^iStmp^ ' " )(r 4l) • (3 - 51) 

To calculate the decay distributions, we take the top (anti-top) quark rest frame in which the 
direction of top quark spin axis coincides with the z— axis. The kinematical variables for top 
quark decay process are, 

st = (0,s t ) = (0,0,0,1), 

P = m(l, 0,0,0) , 

h = (m/2)(l, 0,0,1), (3.52) 

h = (m/2)(l, 0,0,-1), 

/ M = Ef(l, smOj cos 4>, sin. 9i sin. (fi, cos Of) . 

Here the lepton mass is neglected, and the energy of lepton is denoted by Ej. The variables 
9j, cf) are angles to specify the orientation of top quark spin as in Fig. [3.10[ 
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Figure 3.10: The momentum configurations for decay products, lepton and neutrino, in the 
top quark rest frame. The direction of top quark spin is defined by the direction of +z axis. 



It is convenient to introduce the scaled variables, w, y, y, xj and x v 



w 



(l + u) 2 , y 



w 



M 2 W 



9 1 y •) i 

m m z 



7 



w 



M, 



w 



XI 



2Ej 



m 



x,, 



2E V 



m 



The three particle phase space integral can be parametrized in terms of scaled variables xj, y, 
and the angles 9j, 0. 



(PS), 



d 3 b d 3 l d 3 u 4 4 - 

(2tt)3(260) (2vr) 3 (2ro) (2tt)3(2z/°) ( } ( " " " ^' 



m 



2(47r) 4 7o l Jo 



dxj / dy dcosO 



(3.53) 



The differential distribution for the lepton, I, in the decay of the polarized top quark can be 
written as 



dT(s t ) 



it -► bin) 



mg 



'• /IT/ I 2 " 
/ \ V tb\ 



Foist) 



dxidydcos6j y ' 16(47r) 3 \ y 2 J (1 — y/y) 2 + 7 2 ' 

where s t denotes top quark spin states (| or |), and the function F (st) is given by 

F (T) = Xi(l-Xi)(l + cos6 T ), 
Fo(i) = Xi(l-Xi)(l-cose T ). 



(3.54) 



(3.55) 
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After integrating over variable y in Eqn.([T54]), we obtain the double differential (xj — 9j) 
distribution. 



dT{s t ) 
dxfd cos 9\ 



— TflQ / L \ 

(t — > blui) = — - — — — Fn(s t ) arctan 
V ; I6(4tt) 3 \yjj UV ; 



jx T 



[l + -f 2 )y-x T 



(3.56) 



This result shows that the double differential distribution in written in a factrized form, the 
energy distribution and the angular distribution. 

The above formula Eqn. |3.56| will be reduced to a simple form if we take into account that 
the W boson propagator in Eqn. (|3.54 ) can be approximated as for the process t — > blv. 



I 



(l-2/A/) 2 + 7 2 7 



''K 



0( 7 ) 



In fact, the value of 7 is small, ~ 0.02 and we can neglect the contribution from terms of order 
0(7) (the narrow width approximation). In the narrow width approximation for W boson, 
7^0, Eqn. (|3.56|) is reduced to, 



dT 



Ft dxf d cos 0j 



-F{xi) [1 + aj-cos 9j] , 

3xi(l — xj) 
(l-yy(l + 2y) ' 



(3.57) 
(3.58) 



with the region of kinematical valuable xj, 



y < xj < 1. 



(3.59) 



Here IV is the total decay width, and the parameter aj {aj = 1 for the j quark spin state and 
Oij = ~ 1 f° r the j state.) denotes the strength of correlation between the momentum direction 
of decay particle I and the direction of top quark spin axis. 

In the same way, we can calculate the differential distribution for any decay product i 
(i = b, I, v, u, d, W). The differential angular distribution, in general, is written as after 
integrating over x t |36|, |37|, [38| . 

1 dT l + a { cos 9i 
T^l F = o ' ^ 3 - 60 ) 

1 t d COS V{ I 

where 9i denotes the angle between the top quark spin axis s t and the direction of motion of 
decay product i (i = b, I, u, u, d, W). The correlation coefficient c^ may be computed from 
the squared matrix elements (|3.50| )^( [3.51|) . (see Refs.|3^, |3^| for the details.) 
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i 


OLi 


OLi 


I or d 


1 


1.00 


v or u 


(g-l)(2^+lly-l)-12^1ny 


-0.318 


(25+l)(5-l) 2 


W + 


y-'i 

5+2 


0.406 


b 


y-i 

5+2 


-0.406 



Table 3.3: Correlation coefficients a» for the polarized top quark decay in the narrow width 
approximation for W boson. The b quark mass is neglected, and the values m = 175.0 GeV 
and M w = 80.41 GeV are used. 



For the spin-up top quark, the correlation coefficients «j are given in Table |T3| . For the 
spin-down top quark, the correlation coefficients have the opposite sign. 

Using the values in Table. |3~3| , we plot the single particle spin-angular correlation Eqn.( 3~60| ) 



in Fig.3.11 




Figure 3.11: The angular distributions of decay products i, (J 1 d 1 u,i , i ) b and W + ), in the top 
quark decay process. The 9i is the angle between the direction of top quark spin and the 
moving direction of decay product i. 



The lepton and d-type quark are maximally correlated with the top quark spin axis compared 
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with other particles, u,ui,b and W + . It is surprising that the charged lepton or ci-type quark 
come from the decay of the W boson which is less correlated with the top quark spin. 

Here it might be interesting and instructive to investigate spin component of W boson 
contributes to the decay process, (t — > blv). To answer this question, we consider the top 
quark cascade decay t —>■ bW + —>■ blv [|T6| , |34| , |39| . The differential decay distribution of 



a polarized top quark depends on three angles. First is the angle, x\vi which specifies the 
moving direction of the W boson with respect to the top quark spin in the top quark rest 
frame. Second is the angle, ii — xf , between the direction of motion of b quark and lepton in 
the W boson rest frame. Finally the azimuthal angle, $, between the lepton and top quark 
spin around the direction of W boson. The differential decay distribution is 

_1 d 3 T 

T T d cos x l w d cos xf d® 

U 2 (l + cosxW) sin 2 %r + M^(l - cosxW)(l - cosxD 2 



167r(m 2 + 2M 2 ,) - 

+ 2mM w (l — cosxf) sin x^, sin xj^ cos $ . (3.61) 



In Eqn.( |3T6l]) , the first and the second terms are contributions from the longitudinal and 



transverse W boson, and the third term represents the interference between the longitudinal 
and transverse W boson. The contribution of the transverse W boson is suppressed by the 
factor M^r/m 2 , and the interference contribution, by the factor 2M^//m compared to the 
contribution of the longitudinal W. Note that the interference term does not contribute to 
the total decay width because of the factor cos $. 
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Figure 3.12: The 3- dimensional plot of the top quark differential angular distribution, 
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Figure 3.13: The 3-dimensional and contour plot of the top quark differential angular dis- 



tribution, 
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— . Each contour-line shows the angular distribution with the values 
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We plot the 3- dimensional differential angular distribution in cosxf^ verses cosxw i n 
Fig. ^.l2[ Fig. |3.12| (a) presents the contribution of longitudinal W boson, Fig. |3.12| (b) shows 
the transverse W boson. Fig. |3.13| shows 3-dimensional and contour plots of the differential 



angular distribution in cosxi verses cosxw- When the W boson momentum is parallel to 
the top quark spin cosxty = 1, only the longitudinal W boson contributes to the angular dis- 
tribution. When the W boson momentum is anti-parallel to the top quark spin cosxw = ~~ 1? 
the transverse W boson purely contributes. In the region, —1 < cosx^ < 1, both the longi- 
tudinal and the transverse W boson contribute and significantly. Although the contribution 
of transverse W boson to the total decay width is small, ~ 30%, the transverse W boson has 
a significant impact on the angular distribution of the top quark decay products. 

Now we move the to the discussion on the measurement |27|, f|(J of the top quark spin. 
As mentioned before, there are two top quark decay modes, hadronic decay and semi-leptonic 
decay. The rate of the hadronic decay mode is approximately 9 times larger than the rate 
of each leptonic decay mode. However, we must take into account the efficiency to identify 
the d-type quark from u-type (or vice versa) in order to utilize the hadronic decay products. 
Therefore, we consider the the semi-leptonic decay mode of the top (anti-top) quark to simplify 
our discussions. From Eqn.( J3l)0| ), we can obtain the probability P (| | in the cone) that the 
top quark has spin up when we pick up events which satisfy the condition cos 9i > y, 

p / + |. ., v (l + (s t ))(2 + ai(l+y)) 

PIT m the cone) = ; — ; - r , (3.62) 

U ' 4 + 2(^)^(1 + 2/) 

where the probability function P(A\B) denotes the conditional probability of the event A, and 
(s t ) denotes the average value of the top quark spin defined by 

(**)= E s t xP{s t ). (3.63) 

st=-l,l 

P(s t ) is the probability that the top quarks with spin s t are produced. 

It is straightforward to extend this analysis to the top pair production process. The double 

decay distribution for the decay product V from the top quark and the decay product "i" 

from the anti-top quark is given by 

d 2 T 1 

- {1 + oti (st) cos 6i + ai (st) cos Oj 



d cos 9 id cos 9\ 4 

+aj«i (s t St) cos 9i cos 9j} , (3.64) 

where 9i {9j) is the angle between the direction of the top (anti-top) quark spin and the 
momentum of decay product % (i) in the top (anti-top) quark rest frame. As can be seen 
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in the above expression, it is obvious how the spin correlation between the top and anti-top 
quarks is related to the distribution of their decay products. 

When we collect events for which the charged lepton i from the top quark lies in the cone 
defined by cos#j > y, the effective a in the expression Eqn. (|3.60| ) for anti-top quark becomes, 

eff {2(s 1 } + a l (s t st)(l + y)} 

a-/ = ; — 7—. : X «7, M.o5) 

2 + a i {s t ){l + y) " l ; 

which determine the i distribution under the condition mentioned above. Then the distribution 
of the decay product i from the anti-top quark is given by 

1 dT l + a- e// cos# ? 

; 7T = l -■ (3-66) 

Tt d cos &i 2 

By observing this distribution, we can obtain the averaged value of the spin (s t ) and the spin 
correlation (s t St) between the top and anti-top quarks. 
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Chapter 4 

QCD Corrections to Spin Correlations 



A basic consequence of the large top quark mass is its short lifetime. The top quark lifetime 
is t ~ 10~ 24 second and top quark decays very rapidly to the h quark and the W boson 
before QCD long distance effects come in. Therefore, top quarks are produced essentially as 
unbounded fermion. However they still feel the strong interactions and will radiate gluons. 
Especially, there are two "new" effects, which are caused by the QCD radiative correction. 
First, the QCD corrections induce the new structure to 7/Z — t — i vertex. Second, the 
emission of a real gluon with spin 1 leads to the spin-flip effects. Thus it is important and 
urgent to study how the tree level picture will be affected by QCD correction. The QCD 0(a a ) 
corrections to the top quark productions have been studied in several papers |0| |2(| , including 
analyses of the effects on the production angle distributions and their polarization. However, 
in these works, the top quark spin is decomposed in the helicity basis. In previous chapter, 
we studied that the top quark spin decomposition in the helicity basis is not appropriate at 
moderate energies, and the "off-diagonal" basis is optimal spin decomposition in this energy 
region at the leading order. 

In this chapter, we present the differential cross-sections for top quark productions in the 
"generic" spin basis at QCD one-loop level, and investigate which is the optimal decomposition 
of top quark spin at QCD one-loop level at the energy far above the threshold for top quarks. 
We quantitatively analyze spin-spin correlations for top quarks, and examine the spin flip 
effects caused by the QCD radiative corrections. 

4.1 Vertex Corrections and Soft Gluon Case 

In this section, we derive the QCD correction at 0(a s ) to the spin dependent differential 
cross-section for top quark pair production in the soft gluon approximation (SGA) pi] , I2|. In 
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this approximation, only the QCD vertex corrections modify the spin correlations of the top 
quarks. Therefore, in order to investigate a new structure induced by the QCD corrections, 
it is instructive to consider the soft gluon approximation. We give analytic expressions for 
the polarized cross-sections at QCD one- loop level in the soft-gluon approximation using the 
generic spin basis. 

4.1.1 QCD Corrections to Vertex 7/Z — t — t 

The QCD corrections to the cross-sections are given by the interference between the tree and 
the one-loop vertex diagram in Fig. [4.1[ 





Figure 4.1: The tree level and the QCD one-loop contributions to the e e + — > tt process. 



There are both infrared and ultra-violet singularities in the one-loop integrals of the virtual 
corrections. We treat the ultra-violet singularities within the frame work of dimensional regu- 
ralization in D = 4 — 2e, and introduce an infinitesimal mass for the gluon to avoid infrared 
singularities. 

At the one-loop level, the 7 — t — t and Z — t — i vertex functions can be written in terms 
of three form factors Aq, Bq, Cq as follows: 
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(4.2) 



where Q t , 9\y, fn and t^ (t M ) have been defined in the previous chapter. The form factor A 
represents the enhancement effect to the tree vertex structure. So, this factor does not change 
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the tree-level results except for the multiplicative enhancement. Meanwhile, the factors Bq, Cq 
are the coefficients of new vertex structures induced by the QCD corrections. In particular, 
the factor Cq represents the contribution from only the QCD corrected Z — t — i vertex. 
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s + it 

where e = 2 — -D/2, D is the space-time dimension, 7^ ~ 0.57721 is the Euler constant, 
/3 is the speed of the produced top (anti-top) quark and the strong coupling constant is 
a s = ^a s = ^-g 2 with C 2 (R) = 4/3 for SU(3) color. 

We have introduced an infinitesimal mass A for the gluon to avoid infrared singularities. 
The wave function renormalization factor z 2 i n the one-shell renormalization scheme is calcu- 
lated to be, 

* 2 A 2 \ 



Z-2 



1 + a s [ -C uv + In — - 4 - 2 In 
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After multiplying the wave function renormalization factor, the "renormalized" form factors 
A, B and C read 
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(4.9) 



46 



CHAPTER 4. QCD CORRECTIONS TO SPIN CORRELATIONS 

(4.10) 
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/? 1-/3 

where suffices R, I stand for the real part and the imaginally part of the form factors A, B, C 
respectively. In the above expressions, we take only the real part of the form factors. The 
reasons are (1) The contribution of Z boson width is negligible in the region of CM energy 
yfs far above the production threshold for top quarks. (2) We do not consider the transverse 
top quark polarization normal to the scattering plane, since we are interested in how QCD 
corrections modify the tree level spin correlations and which spin basis is the most effective 
for spin correlation studies (Although it is known |12] , |2l| that the transverse top quark 
polarization becomes non-zero when the higher order QCD corrections are included and this 
transverse polarization is very important and related to the phenomena of CP violation, we 
do not consider transverse polarization of the top quarks). 

We get the amplitudes for the e2/ R e R , L — ► t St i s - process using above renormalized vertex: 
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(4.13) 



where the electron (positron) mass is neglected. Note that the contribution from C R vanish 
because it is proportion to the electron mass in the amplitudes. 

Since we do not consider the transverse polarization of top quark, we can define the spin 
of top and anti-top quark in the production plane. We use the generic spin basis defined in 
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Chapter 3 and again, we decompose the top (anti-top) quark momentum t (i) into the two 
massless momenta t±, t 2 (ti, t 2 ). 

Using spinor helicity method, we can easily obtain the squared amplitudes at one-loop 
level, 
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The squared amplitudes for e R e\ scattering process are obtained by interchanging suffices 
L, R as well as |, I in Eqns.(^T4D~(T4T7l). 

The differential cross-section at the one-loop level in the CM frame as shown in Fig. |3.2| is 
given by 
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Here, the angle 9 is the scattering angle of the top quark with respect to the incident electron, 
a is the QED fine structure constant and 7 = l/\/l — (3 2 . The quantities Alr, Bl R , Blr, Dlr 
and Dl R are defined by 
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(4.20) 



The cross-sections Eqns.( [47L8"D and ( |4.19|) contain an infrared singularities (in the form factor 
A R ) that will be canceled by the contributions from the real gluon emission. 
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4.1.2 Spin Correlations in the Soft Gluon Approximation 

Now we consider the cross-sections for the real gluon emission. We isolate the soft gluon 
singularity by splitting the tig phase space into a soft and a hard region. Since the hard gluon 
with high momentum make the momenta of the top and anti-top quarks change, the top and 
the anti-top quarks are not produced back to back. Hence, in the hard gluon region, it is 
difficult to investigate the j/Z — t — i vertex. Then, we concentrate on only the soft gluon 
region to study the vertex structure induced by the QCD corrections ||1], fi^| . 

In the soft gluon approximation, it is very easy to calculate the real gluon contribution. 
As is well known, the amplitude for the soft gluon emissions can be written in the factorized 
form proportional to the tree amplitude. This means that the soft gluon emission does not 
change the spin configurations or momentum of the produced heavy quark pairs from the tree 
level values. Therefore, the QCD radiative corrections enter mainly through the modifications 
of the vertex parts Eqns.( |Q| ) and (O). The cross-section for the soft gluon emissions can be 



rfc°=^max d 3 k f tfj, tfj, V 



written as 

A = J » A • < 4 - 21 > 

a cos a a cos 6 

where the subscript denotes the tree level cross-section. The soft gluon contribution Jir is 
defined by 

J m = -4*C 2 (R)a.J (27r)32fc0 yt k - kJ , 

where uj mSuX is the cut-off of the soft gluon energy, and 

"max < V"s -2m , 

where the maximum value of co> max is constrained by the momentum conservation. At large 
u> m3X , the soft-gluon approximation is spoiled, because it is the leading order of the amplitude 
in the expansion of co> max . Therefore u; ma x should be 

Wmax = x min ( v / s -2m) , (4.22) 

where x m j n is a sufficiently small number. The integral of Jir can be easily performed and we 
obtain 
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By adding the one- loop contributions Eqns. ( fOED and ( |4.19| ) and the soft gluon ones 
Eqn.( P~2"T|) , one can see that the infrared singularities, In A, are canceled out and the finite 
results are obtained by replacing 2A R by 
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inEqns. (^T8|) and ( gig ). 

Thus we obtain the 0(a s ) polarized cross-sections in the soft-gluon approximation using 
the generic spin basis: 
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Since we are interested in maximizing the spin correlations of the top quark pairs, we 
vary the spin angle, £, to find the appropriate spin basis. As shown in the previous chapter, 
we analyze the polarized cross-section in "helicity", "beamline" and "off-diagonal" bases at 
one-loop level in the soft-gluon approximation. 
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1. Helicity basis: 

The helicity basis is given by 



cos£ = ±1. 



Above condition leads to the following differential cross-sections; 
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dcosQ 
3vra 2 



e L e fl ~> *L*L 



{e L e%^t R t R ) 



8s 



da 



dcos9 
3na 2 



8s 



dcosO 
(III + hnf sin 2 9 [(1 - /3 2 )(1 + Sj) - 2S H 

e L e R -* *«*l or tifn) 

1=Fcos0)[/ ll (1=f/?) + /l*(1±/?)] 



(4.27) 



(4.28) 



x [{/ LL (1 =F P) + f LR (l ± /?)}(1 + &) - 2{/xl(1 ± P) + f LR (l T P)}S n ] • 

2. Beamline basis: 

In the beamline basis, the top (anti-top) quark spin is defined in the positron (electron) 

direction in its rest frame. The spin angle £ is obtained by 

cos 9 + (3 



cos£ 



(4.29) 



1 + (3 cos 9 
The differential polarized cross-sections in the beamline basis are obtained by Eqns. (f4.23|) . 

flop and qogp . 
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T l T> 



dcos9 



L^R-* ttfl 



,e r e 



(3 



(3 2 sin 2 



/I*(i-/? 2 )(i + S/) 



2s "/ (l + /3cosfl) 2 
- "{2 - /?(/? - cos9)}f LL f LR + (3((3 + cos9)f 2 LR 



S 



ii 



da 



dcos9 



x e L e R ~~ * Mt/ 



^a 2 \ (3 4 sm 4 9 



2s "7 (l + /3cos#) 2 



/lr(1 + Si + Sn) — fnfLRSn 



(4.30) 



(4.31) 
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l^r ~* *T*4 



,e r e 



37Ta 5 

~27 



/9 



f LL (l + Pcos9) + f LR - 



a-pf 



X 



+ 



/ LL (l + /3cos^) + / LR 

1 



(l-0) s 



l+/3cosfl) 
-1(1H-^ Z ) 



1 + /? cos 



1 + /3 cos 6 
{(/m - f L R){{3 2 {l + cos 2 0) - 2) - Af LR (l + (3cos9)}Sn 



(4.32) 
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3. Off-diagonal basis: 

At the tree level, we have shown that the "off- diagonal" basis makes the contributions 



from the like-spin configuration vanish | lq , [34|. At order 0(a s ), we can find two off- 
diagonal bases for the e^e^ scattering. The first basis is given by 

A L r cos £ - B LR sin £ = . 

This equation leads to the relation between the spin angle £ and the scattering angle 9, 

A LR 



tan£ 



(/ll + /la) VT=F sin 



(4.33) 



B LR f LL (cos 9 + f3) + f LR (cos 9 - (3) ' 

This basis is equivalent to Eqn. flOED in the leading order analysis, and not modified by 
the QCD corrections. The first order QCD corrected cross-sections in this basis are 

da 



dcos9 

da 
dcos9 



(e L e% -► t T t T and ifa) = 



(4.34) 



e L e R -* *T*J or *i*T 



37ra" 

~27 



(3 (y/Al R + B\ R T D LR ) (4.35) 



( V^L^L T D LR ) (1 + Sl ) - 2 f 7 ^* + ^f* T D L . ) g 
The second basis is defined by, 



^l_r + B LR 



(A LR cos£- B LR sin£){l + Si) -2{ 1 2 A LR cos£- B LR sing)S n = 

Therefore the spin angle, £, satisfies the following relation, 

(1 + Si) - 2 7 2 S 



tan£ 



Ul 



A 



LH 



(4.36) 



(4.37) 



BlrO- + Si) — 2B LR S n 
Since Sj and Sjj are C(a s ), there is only difference of 0(a s ) between the Eqn. (|4.33| ) 



and (|4.37|) . Hence, we use the off-diagonal basis defined by Eqn. ([4.33|) . A similar result 
holds for e R e^ scattering. 

In the soft-gluon approximation, the polarized cross-sections depend on the unknown value 
of Wmax fc> r the soft gluon. Before discussing the spin correlations in the soft-gluon approxi- 
mation, we will examine the the uj max dependence of the cross section. Since we are interested 
in the cross-section mainly in the off-diagonal basis, we present the w max dependence of the 
differential cross-section in the off-diagonal basis. Fig.|4~2] shows the ^ max dependence of the 
differential cross-section, where the numerical values of parameters are taken from Table . pTT| 
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in Chapter 3. The running strong coupling constant is obtained by evolving a s (M§) = 0.118 
with 5 flavor to a s (s) with 6 flavor. Here we use the renormalization equation for the running 
coupling constant up to QCD two-loop level. We plot the cross-section at values for the a> max ; 



UJ„ 






i(\/i - 2m) 






1 1 1 
50' 10' 5 



The cross-sections behaves quite uniformly as the value of co> max is changed. This tells us that 
the behavior are qualitatively the same for any reasonable value of CG> max . 
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Figure 4.2: The cj max dependence of the cross-sections in the off-diagonal basis at y/s = 
400 GeV. The DU (UD) component for the e^e + (e] i e + ) process at \/i = 400 GeV is multiplied 
by 100. 
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Figure 4.3: The cross-sections in the helicity, beamline and off-diagonal bases at \/s = 
400 GeV, c<J max = lOGeV for the e~e + — > ti process: tiXn (LR), t^L (R-L), tjJ,L (LL) and 
tfiiji (RR) in the helicity, t^ti (UD), ijjf (DU), t^ (UU) and t^i (DD) in the beamline and 
off-diagonal bases. The suffices "Tree" and "SGA" mean the differential cross-section at the 
tree level and at the one-loop level in the soft gluon approximation. It should be noted that 
DU (UD,UU,DD) component for the ej / e + (e^e + ) process in off-diagonal basis is multiplied by 
100, and that DU, UU, DD (UD,UU,DD) components for the e^e + (e i? e + ) process in beamline 
basis is multiplied by 10 
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e L e+ 


Helicity 


Beamline 


Off-Diagonal 


Dominant Frac. 


0.53610 (LR) 


0.97964 (UD) 


0.99871 (UD) 


Sub-Dominant Frac. 


0.19992 (RL) 


0.00905 (UU,DD) 


0.00129 (DU) 


e R e+ 


Helicity 


Beamline 


Off-Diagonal 


Dominant Frac. 


0.58709 (RL) 


0.99084 (DU) 


0.99699 (DU) 


Sub-Dominant Frac. 


0.15662 (LR) 


0.00406 (UU,DD) 


0.00110 (UU,DD) 



Table 4.1: The fraction of eJ j , R e + cross-sections for the dominant and the sub-dominant spin 
configurations in the helicity, beamline and off-diagonal bases at 0(a s ) in soft-gluon approxi- 
mation. Where y/s = 400 GeV and tu max = 10 GeV. 



In Fig.O, we show the differential cross-sections in three bases, at y/s = 400 GeV and 
w max = 10 GeV. Table |4.1| shows the fraction of the e^m — » ti cross-section for the dominant 
and the sub-dominant spin configurations. The dominant components in the helicity basis 
are txjtn (LR) for e~£e + scattering and tpfi, (RL) for e~ R e + scattering at y/s = 400 GeV. 
This component occupies 53% of total cross-section for e" L e + and 58% of total cross-section 
for e^e + scattering. In the beamline basis, only UD (tftj.) component is dominant for e^e\ 
cross-section, which makes up more than 97%. On the other hand, DU (£;?]-) component is 
dominant for e^ej cross-section, which makes up more than 99%. For e^e\ cross-section in 
the off-diagonal basis, the UU (t|tf) and the DD (£[?j.) components are identically zero. The 
total cross-section is dominated by the UD (£-|-£j_) component which amount to more than 
99% and DU it\t{) component contributes less than 1 %. For e^e\ scattering the UU (t^) 
and the DD (£j,£j.) components are non-zero because we have used the off-diagonal basis for 
e^e^ scattering. However the DU (£j.£t) component is still more than 99% of the total cross- 
section. Although there exist a magnetic moment modification to the j/Z — t — i vertex 
from QCD corrections, this does not change the behavior of the spin dependent cross-sections 
in the helicity, beamline and off-diagonal bases. The QCD corrections, however, make the 
differential cross-sections larger by ~ 30% compared to the tree level ones at this y/s. Thus 
the off-diagonal and the beamline bases continue to display very strong spin correlations for 
the top quark pairs even after taking the QCD corrections into account, at least in the soft 
gluon approximation. 
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4.2 Single Spin Correlations in e + e~ Process 

The result in the soft gluon approximation imply that the QCD corrections do not change the 
behavior of the spin dependent cross-sections in any basis. However, the following two points 
are shortcomings of this approximation. First, soft gluon emission cannot change the spin of 
the heavy quarks. Second, the heavy quark pairs are produced always back to back. These 
points are invalid for the hard gluon, and it is possible that the full 0(a s ) QCD corrections 
might completely change the conclusions of the previous section. In this section, we investigate 
the full 0(a s ) QCD corrections p4] , |4"I|| . Since, in the presence of a hard gluon, the top and 
anti-top quarks are not generally produced back-to-back, it is instructive to consider the 
single heavy quark spin correlations. Namely, we investigate the inclusive cross-section for 
the production of the top (or anti-top) quark in a particular spin configuration. We have 
organized this section as follows. After defining the kinematics and our conventions, we give 
the polarized cross-section for the top quark using a generic spin basis closely related to the 
spin basis of the previous section. The numerical analysis will be relegated to the next section. 



4.2.1 Amplitudes and Kinematics 

The principle result of this section will be the inclusive cross-section for the polarized top quark 
production, e + e~ —>■ t^ or t[ + X, where X = i or ig (The cross-section for the anti-top quark 
inclusive production can be easily obtained from the results in this section.). 





Figure 4.4: The real gluon emission contributions to top quark pair production. 



Since the vertex corrections are the same as in the previous section, all that is required is the 
full real gluon emission contributions. The real gluon emission diagrams to leading order in 
a s are given in Fig.fO. Figure JO also defines the momenta of particles. We also use the 
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spinor helicity method for massive fermions to calculate the squares of these amplitudes for a 
polarized top quark (see Chapter 2). 

The top quark momentum t is decomposed into a sum of two massless momenta t\ , t 2 such 
that in the rest frame of the top quark the spatial momentum of t\ defines the spin axis for 
the top quark. 

t = ti + t 2 , ms t = ti-t 2 , 



where St is the spin four vector of the top quark. The amplitude for Fig. 4.4 is given by 



M (e L CR -► t at ts t g) = (q~ \{il)M-) 



x ufa s t ) 



I (O-LL n , &LR 



2t-k\ 2 

1 



1l + ^-7« )(-t-k + m) lu 



(4.38) 



2t-k 



-f u (i +k + m) [~y1l + -7^1 r 



T a v(i, Si )e:(k) 



where e is the polarization vector of the gluon. T a is the color matrix, and p = p M 7 M . The 
coupling constants au are defined as follows: 

a L i _ e 2 g 

The expressions for the squares of the amplitudes given below have been summed over the 
spins of the unobserved particles (the anti-top quark and gluon) as well as the colors of the 
final state particles. Let us write the square of the amplitude for the top quark with spin "up" 

as 



|M(eZ4-M"(7)| : 

= N C C 2 (R) 

After some calculation, we find 



(i-ky 



Ti + jt 



1 



T 2 + 



1 



(t-k)(t-k) (t-k) 



T, 



(4.39) 



To 



= 4 \a LL \ 2 [{i- k){q ■ k) - m 2 q ■ (t + k)\ (t 2 ■ q) 
+ 4:\a LR \ 2 [(i-k)(q-k)-m 2 q-(i+k)] fa ■ q) 
a LLa* LR ((i ■ k) +m 2 )Tr('j R t 1 t 2 qq) + c.c.\ , 
= \a LL \ 2 [4(t -i)(i- q)fa ■ q) + fa ■ q)Tr( lR qktt) + (t ■ q)Ti{ lL t 2 qktj\ 
+ \a LR \ 2 [4(t -t)(t ■ q)fa ■ q) + fa ■ g)Tr( 7L q k ti) + (t ■ q)Tr( lR t iq ki)} 



+ a LL a 



LR 



{t-t)Tr(^ R t 1 t 2 qq) - (k ■ g)TDr(7 fl ti t 2 q k) 

+-Ti(-f R t 2 t 1 qktq) + -Tr(^ R tit 2 qqkt) 
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+ 0*LL a LR 



+ c.c. 



n 



2|a LL | 



+ 2 \a LR \ 



(t ■ t)Tr( lL t 2 hqq)-(k- g)Tr( 7L t 2 h q k) 

+-Tx(j L t 1 t 2 qktq) + -Ti(^ L t 2 t 1 qqki) 



-m 2 (k ■ q) - 2(m 2 + (t ■ k))(t 2 ■ q) + 2((t + k) ■ q)(t 2 ■ k)] (t ■ q) 
-m 2 (k ■ q) - 2(m 2 + (t ■ k)){tt ■ q) + 2((t + k) ■ q)(t x ■ k)] (t ■ q) 



m 



[a L La* LR {2Tr(j R h t 2 qq) + Ti(^ R kt 2 qq) + Ti(^ R txkqq)} + c.c] 



where 7_r/l = —^ an d that all momentum, p, under the "Tr" operator are understood to be 
p. By interchanging the t± and t 2 vectors in the above expressions, we can get the amplitude 
square for the top quark with spin "down". Since we neglect the Z width, all the coupling 
constants a LI are real. 

To define the spin basis for the top quark, we naturally extend the spin definition of the 
previous section to the present case. The top quark spin is decomposed along the direction s t 
in the rest frame of the top quark which makes an angle £ with the sum of the anti-top quark 
and the gluon momenta in the clockwise direction, see Fig.O. 




Figure 4.5: The spin basis for the top quark in the process e e + — > ttg. 



(We consider the case in which there is no transverse polarization.) To calculate the cross- 
section from Eqn.( |4T39| ), we take the CM frame in which the e + e~ beam line coincides with 
the z-axis, 



v* (1 

2 



g = ^1,0, 0,1 



g=^(l,0,0,-l) 



We define the variables x, y and z which are related to CM energies of the gluon, top and 
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anti-top quarks by 



x = 



1 2k-(q + q) 2t-(q + q) = 2i ■ (q + q) 



s s s 

The momenta of the final state particles, in terms of these variables, are 

k = y^(l-x,(l-x)k),t=y^(l-y,a(y)t), t = ^(1 - z, a(z) t) 

where k, t and i mean the unit space vectors respectively and 



a (v) = ^Ji 1 -y) 2 -a , a(z) = y/(l -z) 2 -a , 



with a = Am 2 /s. Fig. 46 defines the orientation of the top and anti-top momenta and by 
energy-momentum conservation the momentum of the gluon is also determined. 




Figure 4.6: The momentum (unit vectors) configuration of the top and anti-top quarks in the 
CM frame. The momentum of e~ (e + ) is in the +z (—z) direction. 



One can easily obtain the spin four vector of the top quark in the CM frame by boosting 
the spin vector characterized by £ in the top quark rest frame in the direction of top quark 
momentum by j3(y) (the speed of the top quark in the CM frame). The explicit form for 
t\ (t 2 = t — t\) is given by 



m 



t°i = 2"[7(v)(1-/%)cob£): 



;;/ 



t\ = — [l(y) (fi(y) — cos£) sin^cos^ + sin £ cos # cos <^ 



(4.40) 



m 



t 1 = — [^(y) ((3 (y) — cos ^) sin 6 sin ip + sin ^ cos 9 simp] 



2 
m 



*i = 77 [7(y) (/%)-cos£)cos6»-sin£sin6»] 
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where 

y/aj(y) = l-y , y/aj(y) (3(y) = a(y) . 

If we eliminate the gluon momentum k using the energy momentum conservation and use 



the angular variables x > 4> m Fig. 4.6 to specify the orientation of the anti-top quark (if one 
eliminates the anti-top momentum, one can proceed in the similar way by introducing other 
angular variables), the square of the amplitude Eqn.( |4.39|) can be written as 



\M[e- L e\ -> t^ig)\ 2 = - N C C 2 (R) a 2 LL M 1 + a\ R M 2 + a LL a LR M : 



(4.41) 



where Mj are the functions of y , z, angles defined above and the spin orientation £. 
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(t ■ k) 2 
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m 2 (2t ■q-s)(t-q)+s(z + -)(t ■ q)(St ■ q) + (2t ■ q - s)(t ■ q)(St ■ k) 



(t-k)(t-k) 



{£(1-1$ -(!-£(* ■*)}(«■ 9) 
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'l-p-q)}(t-q) 



M 2 
M 3 



-s{^l-y-^)-(l-^Kt-q)}(t.q) 
+ s{~(l-z-~)-2(t.q)}(t.q)(5t.q) 

+ 2(i- q)(t-q)(8t ■q) + {2t-q-s)(t- q)(St ■ t) 
Mi(q++ q St -»• -St) , 

j^ S -(y + ^(-rn 2 s-2(t.q)(St.q) + 2(t.q)(St.q)) 



+ 



+ 



m 



(t ■ k) 2 
1 



~(z + ~) - (s - 21- q){St -q) + {s- 2f • q){St ■ q) 



(t-k)(t-k) 



m 



'■{i(2- a-y - z) - 4(fc- q){k- q)} 



+ 



where St is defined as 



~Z + s{2-a + z)(t ■ q) - s(z + |)(fc ■ q)\ {St ■ q) 

s 2 a \ 

--z + s(2-a + z)(t-q)- s(z + -)(k-q)\ (St ■ q) 



+ (s(l - y){s -(*•?)}- s(3 + z)(t- qj) (St ■ k) 



St = ~(t t - t 2 ) 

(7 
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We transform the scalar product of the momenta to the kinematical variables in the CM 
frame. The functions Mi's are rewritten as 

Mi = — \(l - y) 2 + a 2 (y) cos 2 9 + (I -z) 2 + a 2 (z) cos 2 9 

yz L ■ 

a I — + — J (1 - y 2 + a 2 (y) cos 2 9) - a I — + — ) (1 - z 2 + a 2 (z) cos 2 



+ 2 



jjz y 

2 — 2y — a a 



yz z* 



yz 



" \ / \ n I 2 — 2z — a a\ , , ^ 
— - a(y) cos 9 — 2 a(z) cos 9 

y \ yz z A 



l 

yz 1 
1 



[1 — z — a(z) cos( 



'1 — z — a(z) cos( 



a (I 1 



y(l — a(z) cos#) — z(l — a(y) cos#) (5t ■ t) 
y — z + (y + z)(z — a(z) cos 9) (St ■ (q + q)) 



y(l — a(z) cos 9) + z(l + a(y) cos 9) (St ■ q) 



yz \y z / 

2 r -. -i 

(1 — y + a(y) cos (9) + (1 — y — z)(l — z — a(z) cos 9) (St ■ q) 



M 2 
M 3 



yz 

Mi (cos 9 — > — cos# , cos 9 — > — cos# , St — > —St) 
( 4 4 4 (y + z) 2 (y + z) 2 ^ 



2a { — 

[yz y z yz 



+ 



cos 2 9 h \ - 2a 2 I - + - 



yz 



M z , 



+ — — I — (yaiz) cos9 — za(y) cos9)(St ■ (q + q)) 

yz \y zj 

2 r -i 

— — a(y) cos 9 + (1 — y — z)a(z) cos 9 (St ■ (q + q)) 
yz L J 

/l l\ 2 2 

+ - + - (<Mg-g)) + — [z(y + z) -2(1 -y-z)](St-(q-q)) 
\y z yz 



+ 



yz 



(1 — y)a(z) cos 6* + (3 + z)a(y) cos # (^ • t) ■ 



In the above equations, 9k is the angle between the z axis and the gluon momentum 



(y + z) cos 9k = —a(y) cos 9 — a(z) cos 6* 



Using Eqn.(OU) we find that the products of St with momenta q , q and t are 



^•? = {(1 — y) cos9 — a(y)} cos£ + i/asin^sin^ , 
<5£ • q = { — (1 — y) cos 9 — a(y)} cost; — y/a sin 9 sin £ , 
St -I = { — (1 — z)a(y) + (1 — y)a(z) cos x} cost; 

+ \faa(z) sin £ sin x cos ( 
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The unpolarized top quark production process is given by dropping the spin dependent parts 
(terms proportional to St) in Eqn.( |4.4l| ). As a check we have reproduced the results of Ref. ]T7) 
by putting a LL = a LR = -^3.Q q . 
The cross-section is given by 



da(e~ L e + R - t } tg) = ^|M(eZe+ - t^ig)\ 2 (PS) 3 



(4.42) 



where (PS) 3 is the three particle phase space. 

dH dH d 3 k 



(PSh 



(2TtyS\t + t + k-q-q) . 



(2n) 3 2t° (2ir) 3 2t° (27r) 3 2k° 

We introduce a small mass A for the gluon to regularize the infrared singularities. It is easy 
to rewrite the above phase space integral as, 



(PS), 



(An) 



fv+ fz+(y) 
dy I dz 



z -(y) 



x / dVtd cos xd<t> S cos x 



y + z + yz + a — 1 — 2a x 



a (y)a(z) J 

where dVL = d cos 9dip is the solid angle for the top quark and a\ = X 2 /s. The integration 
regions over y and z are determined by the condition | cosxl < 1, 



y + = 1 - Va , y- = </aa\ + a x , 

z±(y) = y(l-y- - + a x ) +a x ±a(y)^(y-a x ) 2 -aa x 

The integration over the angle is not difficult if one uses the relation 



cos 9 = cos 9 cos x + sin 9 sin x cos < 



The integrals we need are the followings: 



cos 6* 

cos 2 # 

cos 2 9k 

cos 9 cos (ft 
cos 2 6* cos 



2n cos # cos x 



2n 
2n 



cos 2 + -(1 — 3 cos 2 9) sin 2 % 



1 " 2 
cos 2 9 + - (1 — 3 cos 2 i/ / 

2 (2/ + 2) 



a 2 (z) 2 

2 sin x 



7r sin 6 1 sin x , 

2n sin 6 1 cos 9 sin x cos x • 
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Due to the S function in the phase space integral, the angle x is a function of y and z: 

y + z + yz + a — 1 
cos x = - 



a(y)a(z 

4yz( 
sin" x - 



2 4yz(l - z-y)-a(y + zf 



a 2 (y)a 2 (z) 

where we have put a\ to be zero since the A — ^ limit does not produce any singularities 
in the (squared) amplitude. The remaining integrals to get the cross-section are over the 
variables y and z. According to the type of integrand, we group the phase space integrals 
(after the integrations over the angular variables) into four distinct classes { Jj} , {Aj} , {Li} 
and {Ki} [EDI. The individual integrals of these classes are summarized in Appendix. 



4.2.2 Cross-Section in the Generic Spin Basis 

We write the inclusive cross-section for the top quark in the following form. 

<Sa 'e- L e+ - t T X) = ^- ]T (D klmn + a s C klmn ) cos fc 9 sin' 9 cos m £ sin" £ , (4.43) 



dcosO^ «« ' ' 4s klmn 

where D k \ mn are the contributions from the tree and the one-loop diagrams and Ckimn are 
from the real emission diagrams. Let us first write down the Dkimn, 

Ax>oo = /3[fh + fl R + 2af LL f LR ]{l+a a V J ) 

-PWll + hnf - P 2 Ull - hR?]&sVn , 
£2000 = P\fl L + fl R )(X + asV I )+ (3\f LL - f LR ) 2 a s V n , 
Dxooo = 2(3 2 (f 2 LL -f 2 LR )(l + a s V I ), 

^0010 = P 2 (fL-flR)(l+&sV I ), 

#2010 = Pifh-fLdd+Wl), 

D 10W = f3[(fLL + fLR) 2 + P 2 (fLL-fLR) 2 }(l + &sV I ) 

-2(3[(f LL + f LR f - (3 2 (f LL - f LR ) 2 }a s Vu , 
Anoi = -^(III + f L R) 2 [a{l + a s Vj) - (1 + a)a s V n ] , 



'a 
D 1101 = ^(/u-/L)kl + «^)-(l-«)W 



with 



/5 = /?(0) = v^a , a 8 V I = 2A + 2B R , a s V n = B R . 
A R , B R are defined in Eqns.(fL9j) and ( |4.10| ). 
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For the C, 



klmn 



c\ 



0000 



a 



2000 



c 



1000 



Cooio 

C2010 
C1010 



c, 



0101 



Cl 



101 



we find, 
2 (/L + fin) 

2 (fh + f LR ) 
+ 4a /ll/lr 
2 (/x,L — /l_r) 



Ji R - (4 - a) J 3 + (2 + a) J 2 + aJ ± + -R x 



JlR ~ 4^3 — J2 — J\ 



:R, 



1 - a)Jl K - (4 - a) J 3 + (2 - a)J 2 - a.J x - -i?i 



■h + J\ ~ -7R2 
[1 - a)J? R + aiVio - 2(4 - 3a)N 9 

- (4 - 5a)iV 8 - 2iV 7 + 2N 6 + 6N 3 + 2N 2 



C 



1000 



'./../, /l_r) 



AN. 



aN A - aN 3 - aN 2 + (4 - a)iVi + -i? 3 



)i?3, 



+ (// 

C0010 - 2(/ LL - f LR 

UIl + flu) [ 2(2 - a)Jl R + 2aJ 4 - 2(8 - 5a) J 3 

- 2o(l - a)L 8 - 2o(l - a)L 7 - 2(4 + 3a)L 6 + 2(4 - 3a)L 5 
+ 2aL 4 + o(10 - a)L 3 + (8 - 6a - a 2 )L 2 - 2(4 - 3a + a 2 )Li 

+ 2a f LL f LR \ 2 J/ R - 8 J 3 + 4L 6 + 4L 5 + aL 3 + aL 2 - 2(2 - a)L x 



\ UIl + flu) [ 4Jm - (8 + a) J 3 - (8 - a)(l - a)L 7 - (12 + a)L 6 

- 2aL 5 + 2aL 4 + (8 + a)L 3 + (4 - 5a)L 2 + 4(2 - a)L 1 
+ v 7 ^ /ll/lb [ 4JfR - (16 - a) J 3 - o(l - a)L 7 + (4 + a)L 6 

+ 2aL 5 + aL 3 - (4 - 3a)L 2 + 2aLi 

— C1000 + Va (f LL - f LR ) 

-aiVio - 9aN 9 + 2iV 7 - 2iV 6 - -(12 + a)N 3 + -N 2 - (12 + a)N x 

- a(l + a)K 8 + o(l - a)K 7 + 9a(l - a)K 6 - 3a(5 + a)if 5 
- (12 + a)K A + (4 - 50)^3 - 3(4 + 5a)K 2 + 3(4 + a)(l - a)K x 



where we have defined 



n 



J\ 



IK 



Ri 



(2-a)J 6 -aJ 5 , 

2(2 - a)N 13 - aA^i 2 - aN n 
-4a(l -a)L 7 -4(2-a)L 6 



1 - a)L^ 
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+ a 2 L 4 + a(2 + 3a)L 3 - o(2 - a)L 2 + (8 - 8a + 3a 2 )Li , 

R 2 = -4:L G - 4L 5 - aL 3 - aL 2 + 2(2 - a)Li , 

R 3 = a 2 N 10 + 3a(2 + a)N 9 + 4aN 3 - 2aN 2 + 8(1 + a) N x 

+ 2a 2 K 8 - a 2 (l - a)^ - 3o(l - o)(2 + a)if 6 + 6o(l + 2a)K 5 
+ 2(4 + 3a)K 4 + (a 2 + 6a - 8)^3 + 3a(8 + a)K 2 - 12a(l - a)K x . 



Note that the integrals Jl K , J 2 R namely J5 , Jg , Nu , A^ 12 and A^ 13 contain the infrared sin- 
gularity. This singularity is exactly canceled out in the sum Eqn. fl4~43"D by the contributions 
from Dki mn .The numerical value of the coefficients Cki mn and Dki mn are show in Appendix. 
We are now ready to discuss our numerical results. 
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4.3 Numerical Results 

We are now in the position to give the cross-section for polarized top quark production for any 
e + e~ collider. Since we did not specify the spin angle £ for the top quark, we can predict the 
polarized cross-section for any top quark spin. The spin configurations we will display are the 
helicity, the beamline and the off-diagonal bases for center of mass energies \fs = 400 GeV, 
500 GeV, 800 GeV, 1000 GeV, 1500 GeV and 4000 GeV. 



v^ 


400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


4000 GeV 


p 


0.4841 


0.7141 


0.8992 


0.9368 


0.9723 


0.9962 


a s 


0.09804 


0.09564 


0.09096 


0.08890 


0.08539 


0.07794 


<l [Pb] 


0.8707 


0.7728 


0.3531 


0.2312 


0.1047 


0.0149 


offi [pb] 


1.105 


0.8578 


0.3643 


0.2358 


0.1058 


0.0150 


*T,L [P b ] 


1.113 


0.8719 


0.3734 


0.2418 


0.1084 


0.0153 


a% [Pb] 


0.3827 


0.3562 


0.1710 


0.1132 


0.0519 


0.0074 


°%% A [Pb] 


0.4865 


0.3964 


0.1768 


0.1157 


0.0525 


0.0075 


°tm [Pb] 


0.4898 


0.4031 


0.1813 


0.1187 


0.0537 


0.0076 


KL 


0.2778 


0.1283 


0.0575 


0.0457 


0.0347 


0.0259 


Kr 


0.2799 


0.1316 


0.0602 


0.0461 


0.0360 


0.0262 



Table 4.2: The values of /3, a s , tree level cross-section, one-loop cross-section in soft-gluon 
approximation and full one-loop cross-section and Kl/r for eJ y , R e + scattering. The superscript 
0, 1 and SGA stand for tree level, full one-loop, soft-gluon approximation with cu max = (\/~s — 
2m)/ 5 = 10 GeV. 



Table fO| contains the values of the maximum center of mass speed /3, the running a s , and the 
tree level cross-section (J^ L , R , the one-loop total cross-section cr^PA Rl the full one-loop total 
cross-section Ctl/r ^ OT e L/R e+ scattering, Kl/r stands for the fractional 0(a s ) enhancement 
of the tree level cross-section for top quark pair production in eJ / , R e + scattering. 



kl/r 



a T,L/R a T,L/R 
U T,L/R 



At y/s = 400 GeV, the QCD corrections enhance the total cross-section by ~ 30% compared 
to the tree level results. At higher energies, 800 ~ 4000 GeV, the enhancements are at the 
2 ~ 6% level. 

Next we present the numerical values for the coefficients (Dkimn + &sCkimn) in Eqn. (|4.43|) . 
Since the dominant effect of the 0(a s ) corrections is a multiplicative enhancement of the tree 
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level result, we have chosen to write 



Dklmn + O^sCklmn = (1 + K L/R)D klmn 



S, 



klrnn- 



The (1 + ^L/i?)-D^ mn terms are the multiplicative enhancement of the tree level result whereas 
the Skimn give the a s deviations to the spin correlations. The coefficients D\ lmn are given by 
the Dkimn with limit, a s — > 0. The numerical value of these coefficients are given in Table |4.3| 



and 4.4. 





400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


4000GeV 


(1 + K L )D° 000 


1.511 


1.724 


1.722 


1.700 


1.671 


1.644 


{1 + K L )D° 2000 


0.2404 


0.6724 


1.241 


1.383 


1.526 


1.624 


(1 + k l )D» 000 


0.7809 


1.466 


2.126 


2.269 


2.406 


2.494 


(1 + k l )D° 0010 


0.3905 


0.7330 


1.063 


1.134 


1.203 


1.247 


(1 + K L )D» 010 


0.3905 


0.7330 


1.063 


1.134 


1.203 


1.247 


(1 + k l )D% 10 


1.751 


2.396 


2.963 


3.083 


3.197 


3.268 


(1 + K L )D» m 


1.452 


1.502 


1.100 


0.9048 


0.6186 


0.2353 


(1 + k l )D° U01 


0.3416 


0.5131 


0.4650 


0.3970 


0.2807 


0.1091 


SoOOO 


-0.002552 


-0.005000 


0.0005645 


0.004836 


0.01172 


0.02039 


5*2000 


0.007655 


0.01500 


-0.001682 


-0.01450 


-0.03516 


-0.06116 


S'lOOO 


0.02154 


0.03170 


0.01224 


-0.0008102 


-0.02085 


-0.04516 


^OOIO 


0.01094 


0.01785 


0.01586 


0.01304 


0.008357 


0.002274 


52010 


0.01059 


0.01356 


-0.006158 


-0.01791 


-0.03614 


-0.05987 


'S'ioio 


0.004433 


0.005069 


-0.02030 


-0.03590 


-0.06134 


-0.09742 


'S'oioi 


-0.006564 


-0.02193 


-0.04413 


-0.04852 


-0.04952 


-0.03639 


Siioi 


0.007920 


0.004308 


-0.01270 


-0.01734 


-0.02047 


-0.01642 



Table 4.3: The values of (1 + K L )_D°, /mn and Ski mn for e L e + scattering. 
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V~s 


400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


4000 GeV 


(1 + k r )D° 0000 


0.6613 


0.7901 


0.8313 


0.8307 


0.8266 


0.8204 


(1 + n R )D° 2000 


0.1174 


0.3318 


0.6180 


0.6899 


0.7622 


0.8112 


(1 + k r )D% 00 


0.4408 


0.8357 


1.222 


1.306 


1.386 


1.437 


(1 + K R )D» 010 


0.2204 


0.4178 


0.6109 


0.6529 


0.6931 


0.7187 


(1 + k r )D° 20W 


0.2204 


0.4178 


0.6109 


0.6529 


0.6931 


0.7187 


(1 + k r )D% 10 


0.7787 


1.122 


1.449 


1.521 


1.589 


1.632 


(1 + k r )D° 0101 


0.6216 


0.6547 


0.4875 


0.4024 


0.2759 


0.1052 


(1 + K R )D» W1 


0.1929 


0.2925 


0.2673 


0.2285 


0.1617 


0.06289 


5*0000 


-0.001287 


-0.002358 


0.0006297 


0.002742 


0.006079 


0.01024 


5*2000 


0.003862 


0.007074 


-0.001884 


-0.008222 


-0.01824 


-0.03072 


5*1000 


0.01143 


0.01556 


0.003835 


-0.003159 


-0.01367 


-0.02639 


5ooio 


0.005814 


0.008914 


0.007500 


0.006143 


0.003972 


0.001128 


52010 


0.005612 


0.006475 


-0.005121 


-0.01164 


-0.02164 


-0.03468 


5ioio 


0.002248 


0.002290 


-0.01078 


-0.01854 


-0.03107 


-0.04878 


5oioi 


-0.003764 


-0.01122 


-0.02028 


-0.02189 


-0.02208 


-0.01640 


5noi 


0.004154 


0.001587 


-0.007977 


-0.01043 


-0.01198 


-0.009478 



Table 4.4: The values of (1 + k r )D 



,. lmn and S k i mn for e R e 



scattering. 
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e L e + 


400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


4000 GeV 


Roooo 


-0.001689 


-0.002900 


0.0003278 


0.002845 


0.007014 


0.01240 


-R2000 


0.03184 


0.02231 


-0.001356 


-0.01048 


-0.02304 


-0.03767 


-R1000 


0.02758 


0.02163 


0.005757 


-0.0003571 


-0.008667 


-0.01811 


-R0010 


0.02803 


0.02435 


0.01492 


0.01150 


0.006948 


0.001823 


-R2010 


0.02711 


0.01849 


-0.005793 


-0.01579 


-0.03004 


-0.04802 


-R1010 


0.0025310 


0.002115 


-0.006849 


-0.01165 


-0.01919 


-0.02981 


-R0101 


-0.004520 


-0.01460 


-0.04013 


-0.05362 


-0.08005 


-0.15465 


-R1101 


0.02318 


0.008395 


-0.02731 


-0.04369 


-0.07294 


-0.15050 


e^e + 


400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


4000 GeV 


-Roooo 


-0.001947 


-0.002985 


0.0007576 


0.003301 


0.007354 


0.01248 


-R2000 


0.03289 


0.02132 


-0.003048 


-0.01192 


-0.02393 


-0.03787 


-RlOOO 


0.02593 


0.01862 


0.003139 


-0.002419 


-0.009864 


-0.01836 


-R0010 


0.02638 


0.02133 


0.01228 


0.009409 


0.005731 


0.001570 


-R2010 


0.02546 


0.01550 


-0.008382 


-0.01782 


-0.03122 


-0.04826 


-R1010 


0.002886 


0.002041 


-0.007438 


-0.01219 


-0.01955 


-0.02989 


-R0101 


-0.006056 


-0.01714 


-0.04161 


-0.05439 


-0.08003 


-0.1559 


-Rnoi 


0.02154 


0.005427 


-0.02984 


-0.04566 


-0.07406 


-0.1507 



Table 4.5: The ratios, R k i mn = S Hmn /(l + K)Dl lmn , for e L/R e + scattering. 



Table 4.5 shows the ratios, 



Rklmn — Sklmn/v- + K )Dklmn- 

The ratios are never larger than 10% and are typically of order a few percent. Hence the 
0(a s ) corrections make only small changes to the spin orientation of the top quark. 

To illustrate the cross-sections in different spin bases, we present the top quark production 
cross-section in the three different spin bases discussed in the previous Section. One is the 
usual helicity basis which corresponds to cos£ = +1. The second is the beamline basis, in 
which the top quark spin is aligned with the positron momentum in the top quark rest frame. 
In this basis, £ is obtained by Eqn. (|4.29| ). The third corresponds to the off-diagonal basis 
which has been defined in Eqn.( |4.33|) . Note that as (3 — > 1, the helicity and off-diagonal 
bases coincide. Therefore, at an extremely high energy collider, there will be no significant 
difference between these bases. But a remarkable differences among these bases exists at 
moderate energies. Now let us show the polarized differential cross-sections at y/s = 400 and 
500 GeV. 
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1.2 
1.0 r 

S °- 8 

cd 0.6 ■ 
w 

8 0.4 

"B 0.2 - 
■a 

0.0 
1.2 
1.0 

S °- 8 

cd 0.6 
w 

8 0.4 
■a 

"6 0.2 - 
■a 

0.0 



Helicity Basis Beamline Basis Off-Diagonal Basis 




- R (Tree)' 
R (SGA) 

- R (Full) 
L(Tree) 

L (SGA) 
L (Full) 



e e 



+ 



+ 



-- L(Tree) 
L (SGA) 
-- L (Full) 

R(Tree) 

- R (SGA) 

R (Full) 



— — ux1 6 (Tree) 

Dx10(SGA) 

-- Dx10(Full) 




-+- 



+ 



Ux1 (Tree) 
Ux10(SGA) 
Ux10(Full) 
D (Tree) 
D (SGA) 
D (Full) 




-- D (Tree) x 100 

- D(SGA)x100 

D(Full)x100 




+ 



-+- 



+ 



U (Tree) x 100 
U(SGA)x100 
U(Full)x100 
D (Tree) 
D (SGA) 
D (Full) 



-1.0-0.5 0.0 0.5 -1.0-0.5 0.0 0.5 -1.0-0.5 0.0 0.5 1.0 
cosG cosG cos6 



Figure 4.7: The cross-sections in the helicity, beamline and off-diagonal bases at y'i = 
400 GeV. Here we use a "beamline basis" , in which the top quark axis is the positron direction 
in the top rest frame, for each ej j e + and e] t e + scattering. For the soft gluon approximation 
(SGA) cu max = {y/s- 2m)/ 5 = 10 GeV. 
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1.0 
0.8 



Helicity Basis Beamline Basis Off-Diagonal Basis 




1.0 
_ 0.8 

.Q 

S 0.6 

8 0.4 - 

o 
■a 



0.2 
0.0 



R (Tree)' 
R (SGA) 
R (Full) 
L (Tree) 
L (SGA) 
L (Full) 



+ 



+ 



L (Tree) 
L (SGA) 
L (Full) 
R (Tree) 
R (SGA) 
R (Full) 



— — 6 (Tree)' x1 6 
D(SGA)x10 

— D(Full)x10 

U (Tree) 

U (SGA) 

U (Full) 



+ 



+ 



U (Tree) x1 
U(SGA)x10 
U(Full)x10 
D (Tree) 
D (SGA) 
D (Full) 




D (Tree)x10 
D (SGA) 



+ 



-+- 



U (Tree)x10 
U(SGA)x10 
U(Full)x10 
D (Tree) 
D (SGA) 
D (Full) 



-1.0-0.5 0.0 0.5 -1.0-0.5 0.0 0.5 -1.0-0.5 0.0 0.5 1.0 
cosG cosG cos6 

Figure 4.8: The cross-sections in the helicity, beamline and off-diagonal bases at y/s = 
500 GeV. Here we use a "beamline basis" , in which the top quark axis is the positron direction 
in the top rest frame, for each ej j e + and e] t e + scattering. For the soft gluon approximation 
(SGA) cu max = (y/s - 2m) I h = 30 GeV. 



In Fig. [4.7| and Fig. [4.8| , we give the results for y/s = 400 and 500 GeV for both e^e + and e] t e + 
scattering using the helicity, beamline and off-diagonal spin bases. These figures show the tree 
level, the SGA and the full QCD results for all three spin bases. Since the SGA results almost 
coincide with the full QCD results, the probability that hard gluon emission flips the spin of 
the top quark is very small. Clearly, the qualitative features of the cross-sections remain the 
same as those in the leading order analysis. That is the top quarks are produced with very 
high polarization in polarized e + e~ scattering. In Table JO)] and [4.7| , we give the fraction of 
the top quarks in the sub-dominant spin configuration for e2/ R e + scattering, 

a (ele + -> t T A(t or ig)j /a^ L 
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for the three bases. Similar results also hold for e R e + scattering. 



e L e+ 


Helicity 


Beamline 


Off-Diagonal 


Tree 


0.336 


0.0119 


0.00124 


SGA 


0.332 


0.0113 


0.00129 


0(a s ) 


0.332 


0.0115 


0.00150 


e R e+ 


Helicity 


Beamline 


Off-Diagonal 


Tree 


0.290 


0.00564 


0.00167 


SGA 


0.285 


0.00510 


0.00191 


0(a s ) 


0.285 


0.00530 


0.00214 



Table 4.6: The fraction of the e L , R e + cross-section in the sub-dominant spin at y/s = 400 
GeV for the helicity, beamline and off-diagonal bases. For the soft gluon approximation 
(SGA) u max = (y/s-2m)/5= 10 GeV. 



e L e+ 


Helicity 


Beamline 


Off-Diagonal 


Tree 


0.249 


0.0323 


0.00745 


SGA 


0.243 


0.0321 


0.00775 


0(a s ) 


0.243 


0.0311 


0.00907 


e^ e+ 


Helicity 


Beamline 


Off-Diagonal 


Tree 


0.191 


0.0157 


0.00629 


SGA 


0.185 


0.0145 


0.00682 


0(a s ) 


0.185 


0.0156 


0.00816 



Table 4.7: The fraction of the e L , R e + cross-section in the sub-dominant spin at yfs = 500 
GeV for the helicity, beamline and off-diagonal bases. For the soft gluon approximation 
(SGA) cu max = (Vs - 2m)/ 5 = 30 GeV. 



In Fig.|Q|~ |4.12| , we have plotted the similar results for 800,1000,1500 and 4000 GeV colliders. 
The fraction of top quarks in the sub-dominant spin component for ej j e + is given in Tables 



478] and |4.11j Our numerical studies demonstrate that the QCD corrections have a small effect 
on the spin configuration of the produced top (or anti-top) quark for any spin basis. The 
off-diagonal and beamline bases are clearly more sensitive to the radiative corrections than 
the helicity basis. However, in the off-diagonal basis, the top (and/or anti-top) quarks are 
produced in an essentially unique spin configuration even after including the lowest order QCD 
corrections. 
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Figure 4.9: The cross-sections in the helicity, beamline and off-diagonal bases at y'i = 
800 GeV. Here we use a "beamline basis" , in which the top quark axis is the positron direction 
in the top rest frame, for each ej j e + and e] t e + scattering. For the soft gluon approximation 
(SGA) cu max = {y/s- 2m)/ 5 = 90 GeV. 
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Figure 4.10: The cross-sections in the helicity, beamline and off-diagonal bases at ^/s = 
1000 GeV. Here we use a "beamline basis" , in which the top quark axis is the positron direction 
in the top rest frame, for each e~[e + and e^e + scattering. For the soft gluon approximation 
(SGA) u max = (y/s- 2m) /h = 130 GeV. 
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Figure 4.11: The cross-sections in the helicity, beamline and off-diagonal bases at ^/s = 
1500 GeV. Here we use a "beamline basis" , in which the top quark axis is the positron direction 
in the top rest frame, for each e~[e + and e^e + scattering. For the soft gluon approximation 
(SGA) u max = (y/s- 2m) /h = 230 GeV. 
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Figure 4.12: The cross-sections in the helicity, beamline and off-diagonal bases at yfs = 
4000 GeV. Here we use a "beamline basis" , in which the top quark axis is the positron direction 
in the top rest frame, for each ej j e + and e^e + scattering. For the soft gluon approximation 
(SGA) cu max = {y/s- 2m)/ 5 = 730 GeV. 
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Table 4.8: The fraction of the e L , R e + cross-section in the sub- dominant spin at ^/s = 800 
GeV for the helicity, beamline and off-diagonal bases. For the soft gluon approximation 
(SGA) uw = (y/s-2m)/h= 90 GeV. 
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Table 4.9: The fraction of the e L , R e + cross-section in the sub-dominant spin at ^/s = 1000 
GeV for the helicity, beamline and off-diagonal bases. For the soft gluon approximation (SGA) 
Wmax = (y/s-2m)/5= 130 GeV. 
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Table 4.10: The fraction of the e L , R e + cross section in the sub-dominant spin at ^/s = 1500 
GeV for the helicity, beamline and off-diagonal bases. For the soft gluon approximation (SGA) 
w max = (y/s-2m)/h= 230 GeV. 
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Table 4.11: The fraction of the e L , R e + cross-section in the sub-dominant spin at ^/s = 4000 
GeV for the helicity, beamline and off-diagonal bases. For the soft gluon approximation (SGA) 
Wmax = (y/s-2m)/5= 730 GeV. 
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4.4 f LR = f RL = Limit 

In the previous chapter, we investigated why the polarized cross-sections for the beamline and 
off-diagonal bases give similar behaviors. At the leading order, the polarized cross-sections 
under the approximation fin = Jrl = in the beamline and off-diagonal bases coincide, the 
top spin-down configuration for the e~ L e r scattering process is identically zero, alternatively, 
the top spin-up configuration for the e R e + scattering process is identically zero. The smallness 
of the parameters fiR ~ fni ~ is a key point to understand the relation between the 
beamline and the off-diagonal bases. 

Since these analyses were useful to get further understanding on the off-diagonal and 
the beamline bases, we examine the polarized cross-sections at the one-loop level by setting 
Ilr = Irl = 0. 



v^ 


400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


Tree 


1 


1 


1 


1 


1 


0{a s ) 


0.9998 


0.9988 


0.9969 


0.9962 


0.9952 



Table 4.12: The fraction of the e L , R e + — > ty±t cross-section with f LR = f RL = in the 
dominant spin at -y/i = 400, 500, 800, 1000 and 1500 GeV for the beamline and off-diagonal 
bases. Note that cross-section in the beamline basis coincides with that in the off-diagonal 
basis with f LR = 0. 



Table |4.12j shows the fraction of the e~ L , R e + cross-section with fi R = f R i = in the dominant 
spin at the energy region 400 ~ 1500 GeV. The QCD corrections to these processes slightly 
change the tree level result. Although the top spin- up configuration for the e~[e + scattering 
process is 100% at the tree level, its configuration at the QCD one-loop level is less than 
100%. However this change is of order 0.002 ~ 0.05% at energy region 400 ~ 1500 GeV 
and the QCD corrections are almost negligible. We plot the differential cross-sections with 



Ilr = Jrl = in the helicity, beamline and off-diagonal bases in Fig. ^4.13| ~ Fig.[07]. The 
characteristic features of the polarized cross-sections in the beamline and off-diagonal bases 
do not change even after including the lowest order QCD corrections. 
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Figure 4.13: The cross-sections for f LR 
bases at y^ = 400 GeV. 
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Figure 4.14: The cross-sections for /^ = /#£ = in the helicity, beamline and off-diagonal 
bases at y/s = 500 GeV. 
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Figure 4.15: The cross-sections for f LR 
bases at ^ = 800 GeV. 
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Figure 4.16: The cross-sections for f LR 
bases at v^ = 1000 GeV. 
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Figure 4.17: The cross-sections for f^x 
bases at ^~s = 1500 GeV. 



Irl — in the helicity, beamline and off-diagonal 
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4.5 Validity of Soft Gluon Approximation 

In the previous section, we have arrived at the following results. The spin-flip effects caused 
by the hard gluon are quit small, and the QCD corrections produce only the multiplicative 
enhancement to the leading order results. Moreover, this enhancement is almost attributed to 
the soft-gluon effects. In this Section, we consider the possibility that the full QCD corrections 
are reproduced by the soft-gluon approximation with effective cut-off cu max . We repeat the same 
procedure as in the previous analyses in the soft-gluon approximation. At first, let us rewrite 
the spin dependent cross-sections in the generic spin basis as follows, 
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Now we separate the multiplicative enhancement factor from the coefficients, Dl lmn + B\ 
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The deviations 5f^ from the multiplicative enhancement and the ratio Rlh^n are defined by 
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At first we determine the effective energy cut-off oo max of gluon such that the total cross-section 
in the soft-gluon approximation reproduces the full 0(oe s ) total cross-section. 

8=1,1 



/ dcosO J2 r T ' L/ n ( e L/R e+ -* f s x (t or tg)) 
J-i ~f, acosd ' 



s=T,J. 

Next using this cut-off co> max , we compare the above deviation Rf^, with Rumn in the full 
calculation. 

We show the effective u; max and the numerical results of B\ lmn1 Sf^ and -Rfj^ in Ta- 



ble. ^T3| and [4.14| . Compared to Table [4.3[ -^ [4.5| in the previous Section, we find that all values 



are in good agreement with the full calculation results at y^ = 400 and 500 GeV. As the 
total energy get larger, the difference of results between the soft-gluon approximation and in 
the full 0(a s ) calculation become more remarkable. However the spin-flip effects, S^^, are 
still small. These become smaller and negligible, as the total energy get larger. Therefore 
we conclude that the soft-gluon approximation can be an effective tool to represent the spin 
dependent cross-section at 0(a s ). We also present the polarized cross-sections in the helicity, 
beamline and off-diagonal bases. 
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-0.002615 


-0.0005169 


o9GA 
°2000 


0.007821 


0.01787 


0.01724 


0.01365 


0.007846 


0.001551 


o9GA 
°1000 


0.02134 


0.03293 


0.02514 


0.01910 


0.01057 


0.002038 


qUGA 
°0010 


0.01067 


0.01646 


0.01257 


0.009549 


0.005284 


0.001019 


qSGA 
°2010 


0.01067 


0.01646 


0.01257 


0.009549 


0.005284 


0.001019 


c9GA 
°1010 


0.005214 


0.01191 


0.01149 


0.009101 


0.005231 


0.001034 


qSGA 
°0101 


-0.005870 


-0.01676 


-0.02604 


-0.02585 


-0.02235 


-0.01181 


o9GA 
°1101 


0.007915 


0.005621 


-0.005706 


-0.008136 


-0.008968 


-0.005391 


dSGA 

-"-oooo 


-0.001725 


-0.003455 


-0.003337 


-0.002677 


-0.001566 


-0.0003144 


dSGA 
-"-2000 


0.03253 


0.02658 


0.01389 


0.009870 


0.005141 


0.0009552 


tdSGA 
-"-1000 


0.02733 


0.02246 


0.01183 


0.008418 


0.004393 


0.0008172 


tdSGA 

-"-0010 


0.02733 


0.02246 


0.01183 


0.008418 


0.004393 


0.0008172 


tdSGA 
-"-2010 


0.02733 


0.02246 


0.01183 


0.008418 


0.004393 


0.0008172 


pSGA 
-"-1010 


0.002977 


0.004972 


0.003879 


0.002952 


0.001636 


0.0003164 


tdSGA 
-"-0101 


-0.004042 


-0.01115 


-0.02368 


-0.02857 


-0.03614 


-0.05018 


tdSGA 

-n-noi 


0.02317 


0.01096 


-0.01227 


-0.02049 


-0.03195 


-0.04941 



Table 4.13: The values of B\ lmn1 Ski mn and R^j^ n for e L e + scattering 
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400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


4000 GeV 


^max 


13.55 


38.84 


107.2 


150.2 


255.2 


776.1 


R 1 

^0000 


0.1433 


0.1907 


0.08886 


0.07067 


0.05405 


0.04117 


R 1 

^2000 


0.02964 


0.09457 


0.08541 


0.07486 


0.05962 


0.04272 


R 1 

^lOOO 


0.1078 


0.2002 


0.1419 


0.1195 


0.09217 


0.06527 


R 1 

^OOIO 


0.05388 


0.1001 


0.07096 


0.05975 


0.04609 


0.03264 


R 1 

- D 2010 


0.05388 


0.1001 


0.07096 


0.05975 


0.04609 


0.03264 


R 1 

-DlOlO 


0.1729 


0.2853 


0.1743 


0.1455 


0.1137 


0.08389 


R 1 

^OlOl 


0.1325 


0.1546 


0.03437 


0.01420 


-0.001372 


-0.005843 


R 1 

-DllOl 


0.04635 


0.06415 


0.01984 


0.009433 


0.0005521 


-0.002624 


qSGA 
°0000 


-0.001319 


-0.0103680 


-0.009239 


-0.007183 


-0.002033 


-0.0003753 


qSGA 
°2000 


0.003957 


0.01615 


0.01473 


0.01151 


0.008378 


0.001691 


qSGA 
°1000 


0.01137 


0.02918 


0.02083 


0.01558 


0.01141 


0.002253 


QSGA 
°0010 


0.005683 


0.01459 


0.01042 


0.007792 


0.005703 


0.001126 


qSGA 
°2010 


0.005683 


0.01459 


0.01042 


0.007792 


0.005703 


0.001126 


QSGA 
°1010 


0.002638 


0.005783 


0.005487 


0.004327 


0.006345 


0.001315 


qSGA 
°0101 


-0.003445 


-0.02060 


-0.02827 


-0.02725 


-0.02214 


-0.011790 


QSGA 
°1101 


0.004170 


0.004308 


-0.006649 


-0.008751 


-0.008870 


-0.005381 


tdSGA 

-"-oooo 


-0.001994 


-0.005996 


-0.005351 


-0.004217 


-0.001217 


-0.0002282 


jdBGA 
-"-2000 


0.03369 


0.02395 


0.01183 


0.008305 


0.005490 


0.001041 


tdSGA 
-"-1000 


0.02578 


0.01984 


0.009774 


0.006855 


0.004741 


0.0009033 


TfSGA 
-"-0010 


0.02578 


0.01984 


0.009774 


0.006855 


0.004741 


0.0009033 


tdSGA 
-"-2010 


0.02578 


0.01984 


0.009774 


0.006855 


0.004741 


0.0009033 


tdSGA 
-"-1010 


0.003387 


0.002406 


0.001847 


0.001401 


0.001985 


0.0004025 


tdSGA 
-"-0101 


-0.005543 


-0.01367 


-0.02564 


-0.03005 


-0.03579 


-0.05010 


dSGA 

-n-noi 


0.02162 


0.008372 


-0.01426 


-0.02200 


-0.03161 


-0.04932 



Table 4.14: The values of B\ lmn1 Ski mn and i?|^ for e R e + scattering 
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Figure 4.18: The cross-sections in the helicity, beamline and off-diagonal bases at a/s = 
400 GeV. Here we use a "beamline basis" , in which the top quark axis is the positron direction 
in the top rest frame, for each e^e + and e] i e + scattering. For the soft gluon approximation 
(SGA) w max = 13.49 GeV. 
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Figure 4.19: The cross-sections in the helicity, beamline and off-diagonal bases at ^/s = 
500 GeV. Here we use a "beamline basis" , in which the top quark axis is the positron direction 
in the top rest frame, for each e^e + and e] i e + scattering. For the soft gluon approximation 
(SGA) u max = 38.60 GeV. 
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Figure 4.20: The cross-sections in the helicity, beamline and off-diagonal bases at y/s = 
800 GeV. Here we use a "beamline basis" , in which the top quark axis is the positron direction 
in the top rest frame, for each e^e + and e] i e + scattering. For the soft gluon approximation 
(SGA) u max = 106.7 GeV. 
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Figure 4.21: The cross-sections in the helicity, beamline and off-diagonal bases at y/s = 
1000 GeV. Here we use a "beamline basis" , in which the top quark axis is the positron direction 
in the top rest frame, for each ej j e + and e^e + scattering. For the soft gluon approximation 
(SGA) c; max = 149.7 GeV. 
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Figure 4.22: The cross-sections in the helicity, beamline and off-diagonal bases at ^/s = 
1500 GeV. Here we use a "beamline basis" , in which the top quark axis is the positron direction 
in the top rest frame, for each e~[e + and e^e + scattering. For the soft gluon approximation 

(SGA) uj max = 254.8 GeV. 
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Figure 4.23: The cross-sections in the helicity, beamline and off-diagonal bases at ^/s = 
4000 GeV. Here we use a "beamline basis" , in which the top quark axis is the positron direction 
in the top rest frame, for each e~[e + and e^e + scattering. For the soft gluon approximation 

(SGA) uj max = 775.9 GeV. 
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4.6 Fourth Generations Quarks 

The existence of three generations of quarks and leptons is well established in the Standard 
Model. Although there is no strong expectation of additional quarks and leptons, there is 
no reason that additional generations are ruled out. Several models with new generations or 
arguments flavoring new generations have been proposed ||3| . 



We have shown the efficiency and the usefulness of the study of spin correlations in the 
previous Sections and Chapters. Therefore, we consider the pair production of a fourth gen- 
eration quark (b 1 ) with charge —1/3. and investigate the polarized cross-sections at polarized 
e + e~ linear colliders. Here we assume that the b' mass is between 150 GeV and 300 GeV, and 
take the CM energy at 500 and 1000 GeV. In this analyses, we use three bases as shown in 
the previous Sections. One is the usual helicity basis which corresponds to cos£ = +1. The 
second is the beamline basis, in which the top quark spin is aligned with the positron mo- 
mentum in the top quark rest frame. The third is the off-diagonal basis which makes like-spin 
configuration vanish in the e~[e + scattering process. 

To present the differential cross-section, we investigate the b' quark speed, /3, dependence 
of couplings fjjS. These couplings are given by, 

0.97 <f LL < 1.02 , 0.208 < f LR < 0.218 , 
0.433 < f RR < 0.442 , -0.260 < f RL < -0.215 . 

The couplings /jj's have weak (3 dependence. Without Z boson couplings, the couplings 
become fn = Ilr = 1/3. So the contribution from the Z boson is constructive to the fn 
and f RR , and destructive to the f^ R and f R L- 



Next, we plot the (3 dependence of the ratios, Jlr/ III and Jrl/ Irr, in Fig. |4.24| . The ratio 
I'lr/ III has a weak (3 dependence, and it shows the similar behavior as the ratios of the top 
quark couplings, f LR / ' f LL and f RL / ' f RR . However, the ratio Jrl/Irr for b' productions gives 
the large contribution ~ 0.5 compared to the top quark case. As discussed previously, the 
parameters Jll and f RR control the behavior of the dominant configurations in the beamline 
and off-diagonal bases, and the parameters /lr/ III and JrlIJrr decide the behavior of the 
subdominant configuration. Therefore, in the e R e + -^ b'X(b' or b'g) process, the subdominant 
configuration gives larger contribution to the total cross-section than e~[e + — > b'X(b' or b'g) 
process. 
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Figure 4.24: b' quark speed, /3, dependence of ratios, Jlr/ III and /rl/Irr- 



The differential cross-sections in three bases are shown in Fig. [4.25| ~ Fig. |4.29| , and these figures 
contain the results of the tree level and the full QCD. For the e^e + — > b'X(b' or b'g) process, 
the dominant component is top quark spin-up configuration. In the beamline and off-diagonal 



bases, this configuration makes up more than 96% of total cross-section in Fig.[4.25| ~ Fig.[4.29 



While, in the helicity basis, the dominant configuration occupies 77% of the total cross-section 
with b' mass m = 200 GeV at ^Js = 500 GeV, and is less than 92% of the total cross-section 
in Fig. f4.25| ~ Fig. |4.29| . The off-diagonal and the beamline bases represent very strong spin 
correlations. Thus, b' quark is produced in an essentially unique spin configuration even after 
including the lowest order QCD corrections. This result is the same as the top quark case. 

For the e^e + — > b'X(b' or b'g) process, the dominant component is top quark spin-down 
configuration. In the helicity basis, the dominant configuration contributes to 80 ~ 87% of 
the total cross-section. In the beamline and off-diagonal bases, this configuration contributes 
to the 83 ~ 88% of total cross-section. This results show that the dominant configuration 
gives the same amount of contribution to the total cross-section in these three basis. 
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Figure 4.25: The cross-sections for the process e + e — > b'X(b' or b'g) in the helicity, beamline 
and off-diagonal bases at y/s = 500 GeV. The value of b' mass is 150GeV. 
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Figure 4.26: The cross-sections for the process e + e — > b'X(b' or b'g) in the helicity, beamline 
and off-diagonal bases at y^i = 500 GeV. The value of b' mass is 200 GeV. 
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Figure 4.27: The cross-sections for the process e + e — > b'X(b' or b'g) in the helicity, beamline 
and off-diagonal bases at a/s = f OOOGeV. The value of V mass is 200 GeV. 
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Figure 4.28: The cross-sections for the process e + e — > b'X(b' or b'g) in the helicity, beamline 
and off-diagonal bases at a/s = lOOOGeV. The value of V mass is 250 GeV. 
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Figure 4.29: The cross-sections for the process e + e — > b'X(b' or b'g) in the helicity, beamline 
and off-diagonal bases at a/s = lOOOGeV. The value of V mass is 300 GeV. 
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Chapter 5 

Summary and Outlook 



We have studied in detail the spin-spin correlations for top quark productions and spin-angular 
correlations for its decay at polarized e + e~ linear colliders as well as an brief introduction to 
the spinor helicity method which is a powerful calculational technique to obtain analytically 
compact expression for the complicated process. 

For the top quark productions, we presented the polarized cross-sections using a "generic" 
spin basis. We introduced three characteristic bases; "helicity" , "beamline" and "off-diagonal" 
bases, and discussed the differences among these bases. At threshold energy, the beamline basis 
is equivalent to the off-diagonal basis. In the high energy limit, the beamline and the helicity 
bases coincide. However, at moderate energies, there are remarkable differences among these 
bases. We first analyzed the polarized cross-sections in these bases at the center of mass energy 
far above the threshold at the leading order of perturbation theory. We have shown that the 
most optimal spin basis is the off-diagonal basis. In this basis, the like spin configurations, 
UU (if, if) and DD (t±, tj.), vanish, and UD (ij, t±) configuration dominates to total cross- 
section, for instance, occupies 99.88% at the 400 GeV e + e~ collider. This result show that 
the top (anti-top) quark spin is strongly correlated with the positron (electron) spin and top 
quark pairs are produced in a unique spin configuration in the off-diagonal basis. We have 
also discussed why the top (anti-top) quark spin is associated with the positron (electron) 
spin and the relation between the beamline and off-diagonal bases. We have clarified that the 
smallness of parameters fin ~ Jul ~ is a key point to answer these questions. 

For the top quark decay, we investigated the correlation between the direction of top quark 
spin and the direction of motion of the decay products. The consequence was that the anti- 
lepton and ci-type quark are maximally correlated with the top quark spin. The longitudinal 
and the transverse W boson contribution to the top quark decay was also explained. The 
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interference effect between the longitudinal and the transverse W boson produces significant 
contributions to the angular distribution of the decay products. Thus the topology of the top 
quark decay provides us important information to obtain a constraint to the Standard Model. 
We also suggested that the spin correlations can be measured by analyzing the decay products 
of the top and ant-top quark. 

We have extended the analyses of top quark productions at the leading order to the next 
to leading order of QCD, and studied the 0(a s ) QCD corrections to top quark productions 
in a generic spin basis. The QCD corrections introduce two effects not included in the tree 
level approximation. One is the modification of the coupling of the top and anti-top quarks 
to 7 and Z bosons due to the virtual corrections. The other is the presence of the real gluon 
emission process. First, we consider the QCD corrections in the soft gluon approximation to 
see the effects of the modified 7/Z — t — i vertex. Using this approximation, we found that 
the tree level off-diagonal basis continues to make the like spin components vanish and that 
the effects of the induced anomalous magnetic moment are small. Next we analyzed the full 
QCD corrections at the one-loop level. When we consider the three particle final state, the 
top and anti-top quarks are not necessarily produced back to back. So we have calculated 
the inclusive top (anti-top) quark production. We have given an exact analytic form for the 
differential cross-section with an arbitrary orientation of the top quark spin. Our numerical 
studies show that the 0(a s ) QCD corrections enhance the tree level results and modify only 
slightly the spin orientation of the produced top quark. In the kinematical region where the 
emitted gluon has small energy, it is natural to expect that the real gluon emission effects 
introduce only a multiplicative correction to the tree level result. Therefore only "hard" gluon 
emission could possibly modify the top quark spin orientation. What we have found, by explicit 
calculation, is that this effect is numerically very small. The size of the QCD corrections to 
the total cross-section is summarized in Table 4.2 of Chapter 4. At y/s = 400 GeV, the 
enhancement from the tree level result is ~ 30% whereas at higher energies, 800 ~ 4000 GeV, 
it is at the 2 ~ 6% level. Near the threshold, the QCD corrections have a singular behavior 
in j3, the speed of the produced quark. This factor enhances the value of the correction at 
smaller energies. The size of these corrections is reasonable for QCD. On the other hand, the 
change of the orientation of top quark spin is quite small. The deviation from the enhanced 
tree level result is less than a few percent. Actually, the fraction of the top quarks in the 
dominant (up) spin configuration for e^e + scattering in the off-diagonal basis is more than 
93% at all energies we have considered. We can, therefore, conclude that the results of the 
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tree level analysis are not altered even after including QCD radiative corrections except for 
a multiplicative enhancement. Since the the dominant enhancement can be described by the 
soft-gluon approximation, we have also analyzed the validity of the soft-gluon approximation 
and shown the efficiency of this approximation. 

Finally two additional discussions regarding the spin correlations have been presented. One 
is the analysis of the polarized cross-sections under the approximation Jlr/ III = Jrl/ Irr = 
at the QCD one-loop level. The other is the pair production of the "fourth" generation quark 
b' with the charge —1/3. We have estimated the polarized cross-sections for the b' quark 
productions with b' mass between 150 GeV and 300 GeV, at the center of mass energy 500 
and 1000 GeV. For the ej / e + scattering process, the differential cross-sections in three bases 
show the same behavior as the top quark productions except for magnitude of total cross- 
section. The dominant configuration is the spin-up b' quark, which is strongly correlated with 
the positron spin. However, for the e^e + scattering process, the dominant configuration in 
three bases almost give the same contribution to the total cross-section. This fact is caused 
by the ratio of the coupling, \Jrl/ J'rr\ ~ 0.5. This situation is the characteristic feature of 
b' quark productions. 

As has been discussed in this thesis and many articles, there are strong correlations between 
the orientation of the spin of produced top (anti-top) quark and the angular distribution of 
its decay products. And the radiative corrections to these correlation seem to be very small 
as we have explicitly shown for the production process. Therefore, measuring the top quark 
spin orientation will give us important information on the top quark sector of the Standard 
Model as well as possible physics beyond the Standard Model. 
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Appendix A 

Notation and Conventions 



A.l Natural Unit and Convention 



Natural Unit: 



1, (ft 



2tt' 



Transform rules of Unit: 



lkg 



1 m 



1 sec 



5.61 x 10 26 GeV : 



/GeV\ 



V c 2 



5.07 x 10 15 GeV" 1 : 



1.52 x 10 24 GeV" 1 : 



)■ 

' he \ 

K GeVJ 
ft \ 



GeV J ' 
1 TeV = 10 3 GeV = 10 6 MeV = 10 9 KeV = 10 12 eV, 

1 fermi = 1 f m = 1(T 13 cm = 1(T 15 m = 5.07 GeV" 1 , 



'1 fermi) 



■\2 



10 mb = 10 4 /ib = 10 7 nb = 10 10 pb = 10 13 fb, 



1 b 



10 2 fm 2 



-21 



28 2 



10"' 4 cm' = 10"' 8 m z , 



(1 GeV)" 



0.389 mb = 3.89 x 10 8 pb . 



A. 2 Metric Tensor 



Our Metric Tensor g^ in Minkowski Space {x M : ( \x = 0, 1, 2, 3 )} 



9i-iu 



[IV 



^+10 \ 
0-100 
0-10 

V o o o -l j 
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A. 3 Coordinates and Momenta 

The space-time coordinates (t, x, y, z) = (t, x) are denoted by the covariant four- vector: 

rpl^ ^z. ( T T T T 1 

<Aj 1 JU - JU - vC 1 - vC 1 I 

= (t,x,y,z) . 

The covariant four- vector x M is obtained by the sign of the space components: 

Xf, = (x ,x 1 ,x 2 ,x 3 ) 
= (t, -x, -y, -z) 
9^uX 

The inner product is given by 

ry . ry ^^Z rfH 1 rp 

= t 2 -x 2 . 

Momentum vectors are similarly defined by 

p» = (p°,p\p 2 ,p 3 ) 
= (E,p x ,p y ,p z ) 

= g^Pu- 

The inner product is given by 



Pi • Vi = PlP2 M 

= E X E 2 -p l p 2 

A. 4 Pauli Matrices 

• Pauli a Matrices: 



1 \ ( -i\ ( 1 

1 o ' ° 2 = [ i o ' as = -1 



The Relations of Pauli Matrices: 



{<Ji,(Jj} = 2Sij, 
[ai, (Tj] = 2ieij k a k 
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• The Properties of Pauli Matrices: 

a\ = di, det(o-i) = - 1, Tr[(Ji] = . 

• Useful Relation: 

( a ■ a){a ■ b) = (a ■ b) +i er ■ ( a x b) . 

A. 5 7 Matrices 

• Dirac 7 Matrices, 7^ = (7 , 7 ): 

{Y,Y} = tV + tY = 2iT, 

[ 7 ", 7 "] = f^ - 7 ^ = -21(7"", 

7 M = ^t" = (7°, - 7 ) • 



7 5 Matrix: 



75 = 7 5 = itVtV = - 7J^a ct 7 m 7"7 A 7' 7 , 

^ = — ^0123 = + 1 • 



The Relation of 7 5 Matrix: 



{75, Y} = 0. 

• Familiar Presentation of 7: 

7 ° = ( I -1 ) ' 7 = ( -°o- 0^ ) ' 75 = ( ! J ) ' 

where cr is given by 

cr = (cti, 0-2,(73). 

A. 6 Dirac Spinors 

Dirac spinors for fermion with momentum p, U\(p) and for anti-fermion with momentum p 
V\(p) satisfy following equations, (p: momentum , A: spin components , m: fermion mass.) 
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• Equations of Motion: 

(J) - m)u x (p) = 0, 
(^ + m)v\(p) = . 

• Normalization: 

u\(p)u y (p) = 2m 6\x>, 
v\(p)v\>(p) = -2m 8x\> ■ 

• Projection Operators (In the case of Spin Sum): 

^u x (p)ux(p) = f> + m, 
x 

^2v x {p)v x {p) = i> - m . 



A 



A. 7 Spin Vector 



spin vector s M for fermions is taken to be, 



m - (til - , &(?' *) 

S — I . 6 



m ' m(p° + m) J 

where 

(sf =1, s 2 = s% = -1, s» Plx = s-p = 

• The Relations between Spin Vector and Spinor: 

u(p, s)7 M 75«(p, s) = 2m s M . 

• Projection Operators: 

u(p,s)u(p,s) = (j> + m) — ^-, 

v (p, s)v(p, s) — (ft — m) . 



Appendix B 
Useful Formulae 



B.l Dirac 7 Matrices 



7 



t and 



7: 



7 



7 



7 



of 

r 

7' 



7V7 , 

7° , 7 fet 



7 5 crt 

/ 1 u pM 



7 



a^ 



^75 



= -7 fe (where fc = 1, 2, 3) , 

7 °rt 7 o , 



7 <T '^ 



a 



/«/ 



$$$$■■■$ 



7° («75) 7° = ^75 



Matrix Element: 



«(p', s') r u(p, s)\ = u(p, s) r u(p', s') , 



Trace Theorems: 

Tr [7 M1 7 W • • -7 Mn ] 

Tr 1 

Tr [W] 

Tr [^Yl X Y 
Tr [7 A 7 /i 7 iy 7 P 7 CT 7 r 



( for n — odd ) , 

4, 

4 (g^gX* 

g^Tr [7^7 P 7 CT 7 T ] - g Xu Ti [YYYY] + g Xp Tr [7^7 V7l 



g ^X g ua + g ^ g u\ 
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-g Aa Tr [YYYY] + 9 Tr [^YYl 



Tr [7 Ml 7^ 2 • • -7 



AM 



Tr [757 m 7 M2 •••7 Mn 

Tr [tsT'Vt V 
Tr [7 5 7 A 7 M 7 iy 7 P 7 ,T 7 r 



^l/ i 2'p r UM3^,M4 . . . ^ 



/'n 



^MlM3^ r U^2^,M4 . . .^ 



Mn 



_l_ . . . _|_ ^MiMn Tr [ 7 m 7 M2 . . . yw] ( f or n = even ) ? 
( for n < 3 ) 



Xpvp ar 



M{e^ vp g 



£ 



Xpva pr 



£ Xpvr pa 



_ wrap Xp 



£ 



gT + £ 



prap Xv Xrcrp pu 



9 



e "gr 



7 Identities: 



IpT 
1p1 X 1 u Y 

ipYrrr 
ipi x i u i a i p Y 



-27", 

^9xv , 

-27 7 7 , 

2(7V7Y + 7Y7V) 



Version: 



Other Relations: 



7 M ^7 m = 


-- -14, 


l*iW - 


= 4a -b , 


iptyw -- 


= -2^, 


7m 4 V i $Y = 


= 2 <y ^ y ^+ #p 4 $) 



pu^Xarp _ pX^uarp _|_ pa^uXrp _ pr ^vXap _|_ pp^uXar _ q 

(q py e XcrTp = q px e V(yTp + q pa e Xvrp + q p ' T e Xaup + q pp e XaTU ) 



B.2 Conventions for Dimensional Regularization 

• 7 Identities ( D-dimension ): 



7^7 



D 



7m7 V = (2-^)7", 



B.2. CONVENTIONS FOR DIMENSIONAL REGULARIZATION 
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Version: 






4^ + (£> - 4)7 V , 

_2 7 CT 7 ^7 A - (D - 4) 7 A 7 V , 



7 M 7 A 7 I/ 7 <T 7 P 7" = 2( 7 P 7 A 7 V + 7V7 V) + (D - 4)7 A 7 V7 P , 



7 M ^7 p 

7m i ^ 

1»4V ^7 P 

7m I ¥ 4 ii^ 



(2-D)t, 

4a -b + (D-A) £}}, 

-2 i y i - (d - 4) i y i , 

2 (a* ft </+ </ $ ft ft) + (D - 4) y y £ 



Expansion in e: 



a, 

r(e) 

r(e-l) 

T(e-n) 



2 3 

— . o. 

2 
2 " \ , //v,-2\ 



£ 2 o 

-- - (7^ - 1) - \{ll ~ 2lE + ^) + 0(e 2 ) 
e 2 b 



_ iy 



n\ 



1 /V- ^ 

£ fc=i K 



+ 0(e) 



T - Function: 



*r(*) 
r(*)r(i-*) 



r(2*) 
r(i) = r(2) = i,r(3) 



n])-^.r(§ 



7T 

sin7T2; 

7T 

COS7T2; 
2 2*-l 



T r W r, + - 



2, r(n + l) = n\, 



2 - r (§ 



4 



r n + 



-n 



(2n-l)!! ,- (2n)\ 

vtt = r 

2™ 4™n! 

f 4 ) n n! 



7T. 



(2n-l)!! 



7T 



(2n)! 



7T . 
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• ip - Function: 

n— 1 i 

^(1) = -7 , ^(2) = 1-7 ,^(n) = E - " 7 , 



i=i l 



*Q = -In 4 - 7 ,^(±n + ±) = 2^^-1114 - 7 



/3 - function: 



B(x, y ) = fr^i-er 1 (*,y > o), 

Jo 
B(x,y) = B(y,x) , 

B{x,y) = 2 I"' d9 cos 2 *' 1 9 sin 2 ^ 1 , 
Jo 

B(x,x) = 2 1 - 2x B{x, 1 
B(x,y) 



>, 2 j , 



T(x + y) 
Feynman Parametrization: 

1 r 6(1 — «i — «2 — ' ' ' — a r. 



(n — 1)! / da\da.2 ■ ■ ■ da n 



a\CL2 ■ ■ ■ CL n J (ftiOi + 0:202 + ' ' ' + Ot n a r , 



dx- 



2 ■ 



AB Jo [xA + (1 - x)B] 

1 Via + 13) f l , x a ~ 1 (l-x)' 3 - 1 



A«B0 T (a) T (P) Jo "" [xA + (1 - x)B] a+ ' 3 

dx dy- 

Jo Jo 



3 



ABC Jo Jo [xA + yB + (l-x- y)C] 

D-Dimension Integral: 

d°k 1 _ i , T(n-%j „ 

f (_) Tin) L ' 



(2vr) D [k 2 - L] n 


(4vr 


d D k k 2 


% 


(2tt) d [k 2 - L] n 


(4vr 


d D k k^k u 


2 



^(n-l-fJD 



Ls 



+l-n 



f v ; Tfn) 2 



rfn-1- #) „**" 



xn+] V 1)T r#+l-n 



(2tt)^ [P - Lf ( 47r )# V 7 r(n) 

rf D A; k»k u k x kP i , r(n-2-f 



(2tt) d [fc a - L] B (4tt)t 7 r(n) 

x I (>V P + # M V P + £ M V A ) L^ +2 ' n 



4 V '" 

17/ 

[2ti) d '" ■" J Y" J D J (2? 



ff^k^f(k 2 ) = 9 —if^k 2 f(k 
J (2ir) D ^ > D J 2n) D J V 
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B.3 SU(iV) Generators 

The SU(iV) Generators T a (a = 1,2,..., N 2 — 1) are hermitian and traceless matrices, which 
generate the closed SU(iV) algebra. 



rpa rpl 



• rabcrpc 



The fundamental representation is N-dimensional where T a satisfy an additional relation 

{T a ,T b } = ±5 ab + d abc T c , 



which is consistent with the normalization 



arpb 



TrlTT 



-S ab - 



Here d abc is totally symmetric in a, b and c. 
It's given by 



d: 



abc 



2Tr 



a rpb^rric 



{T\T b }T 



According to above equations we get following relation. 



T a T b = J_ 6 + \ d abc T c + ^rabcrpc 

2N 2 2 



The SU(iV) structure constants satisfy following relation. 



rade fbef scfd r 



abc 



and Jacobi Identities 



rabe rede , rcbe rdae , edbe race r\ 

rabejede < rcbejdae _i_ rdbejace __ r\ 

For The SU(3) Generators and structure constants are explicitly written down as follows. 



Ai 



/ 1 

1 

V o o o 

/ 1 

0-10 

V o o o 



\a 



I -i 

i 

V o o o 

1 



1 
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A, 



A 7 



where 




As 



, Ag 




Non-zero components of structure constants are 



./ 



12.3 



; 147 

/4 58 



1 ■ 

-/ 156 
/6 78 = 



= ; 246 

2 ' 



/ 



257 



./ 



:>> 15 



./ 



:-!(i7 



1 

2 ' 



Now we can get some formulae, of course these formulae are independent of the explicit 
representation of SU(iV) Generators and structure constants. 



(T c T c ) tJ 



Tr 



rpcrparpc 


Tr 


rparph 


1 
2 


T' 


irpbrpc 


1 

4 



C 2 (R)5ij , 

lc 2 (G) + C 2 (R) 



8° 



rpa 



abc\ 



Tr 



rparpbrpcrpd 

rabcrpbrpc 



(d abc + if abc ) , 



—S ab 5 cd + -{d abe + if abe )(d cde + if cde ) 

AN 8 

abc fbcdrpd 



_ rabc rbcdrp 

l -C 2 (G)T a . 



C 2 (R) 

racd ebcd 



N 2 



2N 
C 2 (G)5 ab 



N5 



ab 



Appendix C 
Feynman Rule 

C.l QED Lagrangian 

OED ^class t '-'gauge • 

Cdass = 4> (i P - m) i> - -^F^ , 

1 2 

'-gauge = ~~n~ \^ mJ 

• Covariant Derivative and Field Strength Tensor: 

D u = d^ - ieA u , 
r uv = (J^A U — u v A a , 

We split up the above Lagrangian into the free part C$ and the interaction part £/. 

£ QED - C o + C i > 



£ = $(iQ - m)^- 7 (d ll A v -d v A li )(d^A v -d v A' i ) 
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C.2 QED Feynman Rules 



Fermion Propagator: 



Photon Propagator: 



m + is 



v/VWWWV^ 

v k (j, 



— I 



k 2 + is 



g^u 



-(l- a) 



K,,K 



fl^V 



k 2 + is 



Photon - Fermion Vertex: 




win 
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C.3 QCD Lagrangian 

QCD class ~r L, gauge T" ^ ghost i 

Ccuss = ${%ip - m)* - \f; v f^ , 

1 / a\ 2 

L-gauge = — — (^ ^ M J , 

*a\ nafc„,b 



CgHost = W^Y. 

• %, X* represent ghost and anti-ghost field. 

• G M is given following, 

1. Coulomb Gauge: • • • G* = (0, V) 

2. Covariant Gauge: • • • G M = <9 M 

(a) Feynman Gauge: • • • ct = 1 

(b) Landau Gauge: • • • a = 

3. Axial Gauge: • • • G M = n^ ( Light-cone Vector ) 

4. Temporal Gauge: • • • G^ = (1, ) 

• Covariant Derivative and Field Strength Tensor: 

D M = d, - ig s A;T a , 

r p,v = O^Ay — u v A^ — %g s [A^, A v \ 

= (d,A a u - d v Al + g s f abc AlAl) T a . 

We split up the above Lagrangian into the free part £ an d the interaction part Cj 

£ QCD - C ° + £i > 



1 

T 



£ = ^ - m)iP--{d^A a v -d v Al){d' J 'A av -d v A^) 



^(9M«) 2 + (^ X * a )(^x a ) 



d = g s ^T a ^A a u - ^r^Al-d^A^A 



j s fabc/ Pi /[a <3 Aa\ AU-b avc 

2 



o 2 

^jabcjacte^c^c^e _ gjabc^^a^b^c _ 
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C.4 Notation and Conventions 

• Covariant Derivative: 



D^[A] = d,5 l0 - ig s Al[T\ 
Df[A] = dJ ab - g s f abc A c u . 



Field Strength: 








[D„D^ = 


-tg s F« u [T% , 




[D„D u f = 


-ig s F*[Tr b 



= -g s f abc F^, 
F% = d,A a v - d v Al + g s r hc A\Al 



Dual Tensor: 



1 — 2 L p<y i 



P^a = ,u P , fQ pA a _ Q^a + jja^c 



Bianchi Identity: 

[£>„,[£>„,£>„]] + [A,, [£>„£>„]] + [£>„,[£>„,£>„]] = 
Thus we get follow identity. 
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C.5 QCD Feynman Rules 

• Fermion Propagator: 
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Gluon Propagator: 



i5 ij 



— m + ie 



V ,b k (j, ,a 



Ghost Propagator: 



-iS, 



ab 



k 2 + ie 



9i*» - (1 - a) 



k 2 + ie 



—iS. 



ab 



q 2 + ie 
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APPENDIX C. FEYNMAN RULE 



Gluon - Fermion Vertex: 



Gluon 3-point Vertex: 




WslnF^ 



3 

3 



a, (J. 



P2 



<&? 



C 

b,v 



/ ^ 



^ 



c,p 



&/° C [^(Pl ~ V2)p + 9vp{P2 ~ Ps)p + 9pp(P3 ~ Pi) 



C.5. QCD FEYNMAN RULES 
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Gluon - Ghost Vertex: 






Gluon 4-point Vertex: 



g s r c P » 



a,u b ,v 



c 

cX 



^h^r 



/ "% 



d,p 



- *g 2 sf ace f bde (g,u9x P - g, P 9ux) 

ig 2 f a<le f Oe ^ _ g , 
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C.6 Electro- Weak Feynman Rule 



Z° Boson Propagator: 



v/VWWVW^ 

/ k \x 



k 2 -M 2 z + ie 



9iiv ~ (1 - a z ) 



Kfj,K v 



k 2 - a z M z + ie 



W Boson Propagator: 



v/VWVWW^ 

v k (j, 



— I 



k 2 - M^ + ie 



SV - (1 - a> w ] 



k^ky 



k 2 — a w M w + ie 



C. 6. ELECTRO- WEAK FEYNMAN RULE 
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Photon - Fermion Vertex: 




ieQ,(f)Y 



Z° Boson - Fermion Vertex: 




/( 



Qz(M + Qz(fh 



•>R,L, 



Coefficients Qj(f),Q% if): 

When "/" stands for lepton (e~,/i~ and r~), then coefficients Q 7 (/), Qz(f) an d Qz(f) 
are 

2 sin 2 9 W - 1 



While "/" stands for Quark, we get 

Q 7 (U) : 
Q 7 (D) : 



i , Qz(f) = ^r, Q L z(f) 



COS On, 



3 3 cos 9, m 



I ; gf( D ) = -®^L Q|(D) 

3' VZV ; 3cos#, ' ^ zv ; 



2 sin # w cos W 



3-4 sin 2 ^ 
6 sin 9 W cos 6^ 
-3 + 2 sin 2 9 V 



6 sin 9 W cos #„, 
where the suffices U, D mean up-type (u, c, t) and down-type (d, s, b) quarks. 
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APPENDIX C. FEYNMAN RULE 



W Boson - Fermion Vertex: 



H,W 




-Cull 



Coefficients Cu'- 

When "i, J" stand for leptons (e~, /j,~, r~ and (v e , u^, u T ), then coefficients Cu = 1. 



While "i, J" stand for Quark, we get Cu 
Maskawa mixing parameter. 



Vu, where Vu is Cabbibo-Kobayashi- 



Appendix D 

Spinor Helicity Method 

D.l Notation and Convention 

• Pseudo-scalar, Axial-vector: 



1. 1. 

1r = 7^^— , Tl = 7^^— • 



Massless Helicity States: 

u(p,+) = 1ru(p) = \p+) , 
u(p,-) = 1lu(p) = \p-} . 

Dirac Equation: 

p\p±) = (p±\p = 0, (p = 'fp lt ). 



Chirality Conditions: 



(1± 75 )|FF) = 0, 
(p±|(l± 75 ) = 0. 



Normalization and Completeness: 



(p±\r\p±) = W, 
\p+) (p+\ + \p~) (p-\ = P- 
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• Massless Spinor Relations: 

|P±)(P±I = — 2~^> 
(p+\q+) = (p-|g-) = 0, 

(p - \q+) = -(q- \p+) , 
(p - \p+) = ip + \p-) = o. 

• Fiertz rearrangement theorem: 

( J) ±\ 1 »\ q ±) 1 l±Jl = 2 |p=F> <g=F| , 

(p±|y|g±) 7 „I±2* = 2\q±)(p±\ , 
(A+\Y\B+)(C-\^\D-) = 2(A+\D-)(C-\B+) , 

{p-\q+)lR = \q+) (P-| - b+) (q-\ , 

(A - \B+) (C - \D+) = (A - \D+) (C - \B+) + (A- \C+) (B - \D+) . 

• Properties: 

(A T \B±) = -(B T \A±) , 

(A±\ 1 ,\B±) = (BTllMT) , 

(p±hW7/J?=F> = -(lihWT/JPT) (neven), 

(p± hW7/J?±) = (p±|7 /U1 ---7 /Un |p±) (nodd). 

• Polarization Vectors: 

1 (k±\^\p±) 



e%{k) 



'V2 (pT\k±) 



eU^h, = ± {p ^ k±) (\kT)(pT\ + \p±)(k±\ 
Properties of Polarization Vector: 

e ± (k)-k = 0, 

e%{k) = [eUW, 

e x (k)-e* y (k) = -5 xy (A,A'6(±,0)) , 

ex{k)-k = 0, 

k^k v 

£ 4(k)e?(k) = -g»» + L!L. 



D.l. NOTATION AND CONVENTION 
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Massive Spinor States: The massive fermions with momentum p and anti-fermion with 
momentum p 

(P2 Pi) 



«(p,t 

u(p,J, 
u(p,t 
u(p, I 
u(p,t 
v(p,i 

u(p,t 
u(p,J, 



bi+) + |p 2 -> 

bi-) + |P2+) 



m 

[P2 Pi] 



(Pl + I + 
(P1-I + 



[Pi P2] 

m 
(Pi P2) 



/// 



/;/ 



(P2~|, 

' (P2+| , 



bi-) - IP2+) 



[P2P1] 

m 

I- , \ I- \ fePi) 

bi+ - P2-) , 



(Pr 
(pV 



(Pi P2) 

m 
[P1P2] 

m 



(P2 + \ , 
(P2-|, 



with light like momenta pi, P2, Pi and P2, 

p + ms p 

7~2 

p + msp 



Pi 
Pi 



P2 = 
P2 = 



_ p — ms p 



2 
p — msp 



where s p is spin vector of the massive fermion, and Sp is spin vector of the massive 
anti-fermion. 

Polarization Vectors for Massive Boson: 

Three polarization vectors e^(k) are given in terms of two light-like vectors, 



k^ — ms^- 



A' 9 



2 



The polarization vectors for the incoming gauge boson: 

(fcii^lAki) 






(A; 1 +| 7 ^|A; 1 +)-(A; 2 +|7' t |A; 2 +) fcf - A# 



2m 



/;/ 
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Appendix E 

Phase space integrals over y and z 



The phase space integrals necessary to derive the cross section are summarized in this Ap- 
pendix. Although many of them have already appeared in the literatures fl20H , we will list 
all of them below for the convenience of the reader. After the integration over the angular 
variables, we are left with the following four types of integrals: 

dydz 



J; 



la 



dydzfi(y,z) 
dydz 



Ni 



■fi{y,z) 



K, 



1 -y) 2 -a 

dydz 



--fi(y,z) 



;fi{y,z) ■ 



The infrared divergences are regularized by the small gluon mass A and f3 



VT 



1 — Am 2 /s. For the type Li integrals, a shorthand notation u> = J(l — y / a)/(l + y/a) is 
used. The Ki integrals have a spurious singularity at the upper bound of the y integral, 
y + = 1 — y/a. Since this singularity turns out to be canceled out in the cross section, we 
regularize each integrals by deforming the integration region as y + — > 1 — y/a — e [p0|1 . Li 2 is 
the Spence function. 



Class J Integrals: 

J\ — / dydz 

J 2 = dydz — = dydz — 
J z J y 



H^r)-rH«hm 



J 3 = / dydz — = / dydz - = —(3 + 1 1 — -a ) In I — 



!+£' 

P; 
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Ji 
J 5 



dydz — 

z 2 



dydz — = -p - In 

y 2 a V 1 -/ 5 , 



dydz — = / dydz — 
J y 2 J z 2 



J fi 



y 

W f\ n ?L + 21na - 41n/3 - 41n2 + 2^ + 2 (l - -) In (i±| 

1 



yz 
In — - In a + 41n/3 + 21n2 j In ( ^-^ 



Li^ 



2/5 



1-/3, 



Li<; 



2/3 x 
.1 + 0, 



Class N Integrals: 



N, 

N 2 
N 3 



dydz 



i -y) 2 -a 



r- 1 1-J^i 

1 - V a - 7T a l n ^— 

2 V -v/a , 



dydz 



^(l-y) 2 -a y 
dydz y 

l-y) 2 ~a z 



1 , 

- lna + In (2 — \fa 



2-,/E 



H^- h (H*M^H-iW^ 



N 4 
N 5 



+ - In 2 - - v^ + 1 + Li 2 

4 4 



dydz 



1 + P\ t- A-/3 N 



+ Li 2 



-2Li 2 Y. 



l-y) 2 -a z2 a 
dydz 



1 -y) 2 -a 



!) 



1 

16 



N. 



-a 2 In a + 2 a 2 ln(2 - y/a) + 4(2 - ^Ja) 2 - 4 (2 - at) 
dydz 



1 
32 



1 -y) 2 -a 
12 - (2 + a) 2 - _ + ^ a 2 + 2(8 - a)a In 



2-v^ 



2-^/a 



iV 7 



dydz 



i -y) 2 -a 



133 



1 
1 + -a 
2 



Li2 (l±l) + Li2 (lz£ 



-lnfi±^U 



2 / 



2Li 2 I -V</ 



In 2 -o 



+ 3wln(- -J +-(18 + a) 



--(20 - a)v^ + f 3 - a + —a 2 ) In 
8 V Id 



1 + /3J 8 

/2-^ N 



iV s 



iVo 



ofa/Gb 1 „, /2 — Va N 

2 In - 



>/(i - y) 2 - a y 

dydz 1 

v/(l"!/) 2 -a J 
Li2 (i±^) + Li 2 (I^T 



2Li 2 



2-y/al 4 



+ - In 2 a 



+ ln 2 



iVi 



10 



dydz 



ln(l + /?) ln(l - /3) 



iV- 



ii 



l-y) 2 -a ^ 2 a 
dydz 1 



l-y) 2 -a r 
A 2 



N- 



12 



■In— + lna + 21n(l - y/a\ - 41n(2 - y/a) + 21n2 - 2 

s 

cfa/ob 1 



i -y) 2 -a 



N, 



1.3 



In— -lna + 21n(l-v^)-^-^lnf^-^ J +21n2 
s 13 V 1 -/ 3 / 

efa/dz 1 



l-y) 2 -a J/* 



^ ln (ttI. 



In— + ilna + 41n(2- Va)-41n(2/?)-21n(- — ^ 
s z y 1 + p 



1 / (I-/?) 2 \ 2 / v^(2-y^) \ / 2^(1-^/5) * 



2 
1 



Li, 



/ > /o(2-Vo) \ 

I (1+/5) 2 ; 



1 + /3 



Li 2 — - — + Li 2 



Li 2 



+Li2 ( <;-^ 



7T 
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Class L Integrals: 



Li 
L 2 



dydz 1 — a (\ + (3 s 

[i-y) 2 -a = 4^ n l v T^ / 

dydz z 



'l-yf -a y 

12 24 



-1 



4-o (4-o) 2 y 
cfa/ob y 



In 



" 1 + /A 2/3 
1-/3/ 4-o 



(1 -y) 2 -a z 
1 



2 1 1-0 



lna + ln(2 + y/a) -hi(l + y/a) -21n2 



+ 



+ 



1 - y/a 



Li 2 (a;) + Li 2 I p= a; | - (a; 



2-x/o 



-a; 



Li 2 



1+a; 



Li 2 (2 + y/a) 



1 +00 



Lio 



2v^ 



(2 + v^)(l + a;) 



— (a; — > —a; 



£4 

u 
u 

u 



dy dz y 2 

[1 — y) 2 — a z 2 



In 






2/3 



dydz 



y 



[l-y) 2 -a 

dydz 
[l-y) 2 -a 

dydz 1 

l-y) 2 -a z 

dydz y 
'1 — y) 2 — a z 2 



, , 1 6 

- 1 + -a 

1 2 4- 



3 a 



i l ,ii±4U/ ? 



4 
1 



In 



/9 N 



1-/3, 



+ 



20 



Li 2 (cj) + Li 2 



'2 + ^ 

2-y/E 



(a) \ — (LO 



-co 



a 1 1-/3, 



Class X Integrals: 



#i 



dydz 



{(l-y) 2 -o}3/2 



2(1 -v^) 2^(1-^/5) 



#9 



v^(2 - v^) 2 
dydz 



e(l + -v/o) (4 - a) 



4a ^{l_+ y/E){2 - y/E) 2 



2y/aa 



{(l-y) 2 -a}3/2 



.'/ 



2(1 -y/E) 2 2y/a~{l-y/a) 
In 



rr , (1 + y/a)(2-y/af 



y/E(2 - y/E) 2 e(l + y/a) (4 - a)' 



In 



2 v /aa) 
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K» 



K A 



K* 



K. 



K 7 



+ In 



_Va_ 

1 + Va) 

dydz 



{( l_ y) 2_ a} 3/2 



2(1-JE)\ 2Ja(l--Ja) a(a 2 + 20a - 32) , (1 + Ja){2 - JE) 2 
In ; =. 1 — — In — 



(2 - ^y 



2^E 



e(l + y/a) 2(4 -a) 3 

2o(l - ^/a) 



2y / aa 



-iln 

2 (1 + v^) (4-a)(2-v^) 2 



+ 



dydz 



{( l_ y) 2_ a} 3/2 



y2 



2(1 - v/a") 3 ^ 2yg(l - yg) a 2 (12 - 7a) ^ (1 + y/a)(2 - ^E) 2 



(2 - ^)3 



e(l + v 7 ^) 
2 



1^ 2y^ 

2 n (1 + y/a) ' (4 - a 



+ 



2(4 -a) 3 

;i - j- a f 



2^/aa 



(2 - V5) 5 



+ VO-1 



dydz 



y 



{(l-y) 2 -a} 3 / 2 ^ 



2(1 - v^) ln Vg - yjj) (4 - 3a) ^ (1 + v/a")(2 - ^) 2 



a(2 — yfa) e(l + ^/o) 
2v^ 



-iln 

a (1 + v^a) a(2 — ^/a) 



a(4 — a) 2y / aa 

1 + /5\ 



4(1 -V5) 2/3 

a \l-(3 



dydz 



1 



{(I-!/) 2 
2 



a 



j.3/2 



a(2 — y^a) 



b 2ya(l - y/g) | (4 + a) ^ (1 



y/a){2-y/af 



c(l 

+ 



4 



a (1 + -y/a) a(2 — y/a) 



dydz 



y 



a(4 — a) 

a/3 V 1 -/ 9 / 
2 2^(1-^/5) 



2^/aa 



{{l-yf-ay/ 2 z 2 a^E e(l + y/a) 
dy dz y 2 



{(l-t/) 2 -a} 3 / 2 z 2 
2(1 - yg / Vgl-Va) ' 
11 V e(H-v^) , 



aA/a 



a \1 + Va, 
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Appendix F 
Total cross section 



Let us examine the unpolarized total cross section including the top quark pairs using our 
formulae. It is given by 



a T (e~e + ^t + t + X) = - [cr T (ele% ^ t + t + X) + cr T (e R e+ -^t + t + X) 



Note that only k = 0, 2 and I, m, n — terms in Eq. (|4.43| ) contribute to the total cross section. 
Integrating over the angle 8, we get 



Mez4 -> t + t + X) 



ira 



P 



(III + Jlr) 2 (3 - (3 2 ) 
6 



x 1 + a s { Vi 



2 o2 



3-/3 



2 ! 8 



8 



/? J * /? 3 /?(3-/? 2 ) 



J 2 



+ 2(f LL -f LR yp 



x (i +4 .{v I+%+ ^-|j, + ?e^fflj i+ ^; 



The cross section for the process e R e~j^ is obtained by interchanging the coupling constant 
L «-» i? in the above expression. Parametrizing the total cross section as 



/Us) = 






a T (e~e+ -»• i + t + X) 



7T 



where cr p t = 47ra 2 /3s, we get the following numerical results at the CM energies y/s 



400 , 500 , 800 , 1000 , 1500 GeV. 
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v^ (GeV) 


R^(s) 


C2(R)R^(s) 


400 


1.0083 


8.9963 


500 


1.4190 


6.0267 


800 


1.6864 


3.3999 


1000 


1.7317 


2.8422 


1500 


1.7714 


2.2911 



These are consistent with the results in Ref. 44 



Appendix G 
Numerical Tables 



We have shown the numerical tables for e r L , R e + — »■ tX(t or £g) process in the previous Chapters 
and Sections. In this Appendix, we present the numerical tables, we used to analyze the spin 
correlations because of its usefulness. 

Now we show the numerical value of the coefficients Ckimn, ^kimn using following inputs. 



m = 


= 175 


(GeV) 


M z -- 


= 91.187 


(GeV) 


a = 


= 1/128 




a s {M 2 z ) - 


= 0.118 




sin 2 9 W = 


= 0.2315 





The multiplicative enhancement, kl/^l, the top quark speed f3, and QCD strong coupling 
constant a s , a s = C 2 (i?)/(47r)a s . 



v^ 


400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


4000 GeV 


k(L) 


0.2778 


0.1283 


0.05747 


0.04574 


0.03474 


0.02593 


k(R) 


0.2799 


0.1316 


0.06023 


0.04790 


0.03599 


0.02618 


P 


0.4841 


0.7141 


0.8992 


0.9368 


0.9723 


0.9962 


a s 


0.09804 


0.09564 


0.09096 


0.08890 


0.08539 


0.07794 


a s 


0.01040 


0.01015 


0.009652 


0.009433 


0.009060 


0.008270 
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G.l e~ L e + ^U + X 



Now we define coefficients D° klmn , D l klmn} S klmn: S l klmn and S] 



klrnn' 



' e L e R ~>tiX) = — J2 ( D klmn + 6t, S C k lmn) COS fc sin' 9 COS™ £ sin™ £ , 



dcos9 y L n ' ' As ,, 

klmn 



with 



D = D° -I- D 1 

J-sklmn — ^klmn ' ^klmrf 

^2000 = P VJLL + JLRIi 
-^1000 = ^p (j LI/ — J LR ), 

-^0010 = P KJLL ~ JLRJi 
-^2010 = P \JLL ~ J. 



Li?)' 



£>1010 = p-[(fLL + fLR) 2 +P 2 (fLL-fLR)% 

A)101 = —j= a (fLL + Ilr) , 
\/0 



1101 ~~ /- a \JLL JLRJ- 



a 



Dlooo = PlfL + fLR + lafLLfLRKasVj) 

-PWll + Ilr? - PUll ~ f L R) 2 ]& s V n , 

DloOO = /F(fL + fhd(P<sVl)+/P{fLL-fLR) 2 & a V II , 

DU = 2(3 2 (f 2 LL -f LR )(a s V I ), 
^ooio = WlL-flR)^), 
Dloio = ^(flL-fDi&sVj), 

Dl 0W = P[(fLL + fLR? + P 2 (fLL- fLRH&aVj) 

-2/3[(/ LL + / Li? ) 2 - /5 2 (/ L l - f L R) 2 ]a s V H , 

£>oioi = ^ F {fLL + fLR) 2 [a{a s V I )-{l + a)a s V II ] , 
\ a 



Dl 101 = ^=(f 2 LL - f 2 LR )[a(a s Vr) - (I - a)a s V n ] . 
\/a 
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i — ^( e z4 -»■ *t*) = — i — Z 5 ^™™ cosfe sin/ cos?n £ sin ™ £ > 

a cos 4s , , 



with 



W E (^ m n + ^ m n)cos^sin^cos m esin"e 



4s 

klmn 



-/5[2(/ll + /l/?) 2 - /9 2 (/ll - fLR?]Sn , 
5 2ooo = P (III + Ilk) Si + P (III ~ Ilr) S h , 
^1000 — %P {III ~ fLR)Si > 

A)010 = P \JLL ~ JLRJ^I , 
-°2010 = P (/ii — JLRJ^I ) 
#1010 = /3[(/LL + / M ) 2 + /3 2 (/LL-/Li ? ) 2 ]^ 

-mifLL + fLR? - /3 2 (/LL - fLR?]Sn , 

^oioi = -^(fLL + fLRflaSj-il + ^Sn], 
\/a 



B\ W1 = ^(fh-flRjlaSr-il-^Su] 



Uhlmn * ™s^klmn — ^klmn ' ^klmn 

= (l + KJ^L + ^mn ••• Ful1 

— (1 + K )^Hmn + ^Hmn • • S.G.A 

Therefore, S klmn , S\ lmn and Sl lmn are given by 

^klmn ^klmn ' ^ s^ klmn 

'-'klmn ~ ^klmn ' ^s^klmn ~ K ^klmn 

Q2 = nl _ K F)0 

klmn klmn klmn ' 



ol 

ryFull klmn 

It 



klmn (1 + k)D° 



It 



SGA '-'klmn 



klmn 

C2 

Or. 



kl ™ " (l + K)^^ 
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0(a s ) Analysis 



v^ 


400 GcV 


500 GcV 


800 GeV 


1000 GeV 


1500 GeV 


4000 GeV 


Coooo 


-11.66 


-30.13 


-58.68 


-71.01 


-94.84 


-167.3 


C2000 


-1.872 


-12.07 


-44.72 


-61.57 


-92.84 


-175.3 


C1000 


-6.010 


-25.83 


-74.58 


-98.212 


-142.3 


-263.3 


C0010 


-2.988 


-12.72 


-36.28 


-47.68 


-69.09 


-128.7 


C2010 


-3.023 


-13.14 


-38.56 


-50.97 


-74.00 


-136.2 


C1010 


-13.60 


-42.69 


-105.4 


-135.4 


-191.9 


-349.5 


C0101 


-11.28 


-26.85 


-39.76 


-40.73 


-38.70 


-27.28 


C1101 


-2.637 


-9.126 


-16.75 


-17.79 


-17.47 


-12.60 


-Doooo 


1.630 


2.025 


2.289 


2.374 


2.542 


3.048 


-D2000 


0.2675 


0.8100 


1.671 


1.949 


2.332 


3.012 


-D1000 


0.8650 


1.760 


2.858 


3.194 


3.674 


4.626 


Anno 


0.4325 


0.8799 


1.429 


1.597 


1.837 


2.313 


-D2010 


0.4325 


0.8799 


1.429 


1.597 


1.837 


2.313 


D1010 


1.890 


2.835 


3.960 


4.324 


4.874 


6.060 


-D0101 


1.563 


1.753 


1.439 


1.241 


0.9196 


0.4246 


Dnoi 


0.3770 


0.6100 


0.6139 


0.5475 


0.4185 


0.1969 




V~s 


400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


4000 GeV 


-^0000 


1.183 


1.528 


1.629 


1.625 


1.614 


1.603 


D° 


0.1882 


0.5960 


1.174 


1.323 


1.475 


1.583 


-^1000 


0.6111 


1.299 


2.010 


2.169 


2.325 


2.431 


-^0010 


0.3056 


0.6497 


1.005 


1.085 


1.163 


1.215 


-^2010 


0.3056 


0.6497 


1.005 


1.085 


1.163 


1.215 


-^1010 


1.371 


2.124 


2.802 


2.948 


3.089 


3.185 


-^0101 


1.137 


1.331 


1.040 


0.8652 


0.5978 


0.2294 


-^1101 


0.2674 


0.4548 


0.4397 


0.3797 


0.2712 


0.1063 


5*0000 


0.3260 


0.1910 


0.09417 


0.07918 


0.06780 


0.06193 


5*2000 


0.05993 


0.09144 


0.06577 


0.04599 


0.01608 


-0.02013 


5iooo 


0.1913 


0.1984 


0.1278 


0.09841 


0.05992 


0.01786 


5ooio 


0.09583 


0.1012 


0.07363 


0.06265 


0.04874 


0.03378 


^2010 


0.09548 


0.09688 


0.05161 


0.03170 


0.004244 


-0.02836 


5ioio 


0.3852 


0.2775 


0.1408 


0.09893 


0.04598 


-0.01484 


5oioi 


0.3092 


0.1488 


0.01564 


-0.008947 


-0.02875 


-0.03045 


5noi 


0.08220 


0.06263 


0.01257 


0.00001984 


-0.01105 


-0.01366 
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400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


4000GeV 


(1 + ^o°ooo 


1.511 


1.724 


1.722 


1.700 


1.671 


1.644 


(1 + ^2°000 


0.2404 


0.6724 


1.241 


1.383 


1.526 


1.624 


(1 + «Mooo 


0.7809 


1.466 


2.126 


2.269 


2.406 


2.494 


(1 + ^oVo 


0.3905 


0.7330 


1.063 


1.134 


1.203 


1.247 


(1 + ^2°010 


0.3905 


0.7330 


1.063 


1.134 


1.203 


1.247 


(1 + k)D\ 010 


1.751 


2.396 


2.963 


3.083 


3.197 


3.268 


(1 + k)D° 0W1 


1.452 


1.502 


1.100 


0.9048 


0.6186 


0.2353 


(1 + k)D» 101 


0.3416 


0.5131 


0.4650 


0.3970 


0.2807 


0.1091 


ol 

°0000 


-0.002552 


-0.005000 


0.0005645 


0.004836 


0.01172 


0.02039 


ol 

°2000 


0.007655 


0.01500 


-0.001682 


-0.01450 


-0.03516 


-0.06116 


ol 

°1000 


0.02154 


0.03170 


0.01224 


-0.0008102 


-0.02085 


-0.04516 


ol 

°0010 


0.01094 


0.01785 


0.01586 


0.01304 


0.008357 


0.002274 


°2010 


0.01059 


0.01356 


-0.006158 


-0.01791 


-0.03614 


-0.05987 


ol 
°1010 


0.004433 


0.005069 


-0.02030 


-0.03590 


-0.06134 


-0.09742 


ol 
°0101 


-0.006564 


-0.02193 


-0.04413 


-0.04852 


-0.04952 


-0.03639 


ol 
°1101 


0.007920 


0.004308 


-0.01270 


-0.01734 


-0.02047 


-0.01642 





400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


4000 GeV 


jdFuU 

-"-oooo 


-0.001689 


-0.002900 


0.0003278 


0.002845 


0.007014 


0.01240 


tdFuU 
^OOO 


0.03184 


0.02231 


-0.001356 


-0.01048 


-0.02304 


-0.03767 


tdFuU 

-"-iooo 


0.02758 


0.02163 


0.005757 


-0.0003571 


-0.008667 


-0.01811 


tdFuU 

-"-ooio 


0.02803 


0.02435 


0.01492 


0.01150 


0.006948 


0.001823 


pFull 
-"-2010 


0.02711 


0.01849 


-0.005793 


-0.01579 


-0.03004 


-0.04802 


tdFuU 

-"-ioio 


0.0025310 


0.002115 


-0.006849 


-0.01165 


-0.01919 


-0.02981 


jdFuU 

-"-oioi 


-0.004520 


-0.01460 


-0.04013 


-0.05362 


-0.08005 


-0.15465 


tdFuU 
-"-1101 


0.02318 


0.008395 


-0.02731 


-0.04369 


-0.07294 


-0.15050 



144 



APPENDIX G. NUMERICAL TABLES 



Soft Gluon App. 

We define the u max to satisfy the following relation. 



T ' L/Rf e7, D e + ^t s X(iortg)) 



dcos6 V L i R 



s=T4 





400 GcV 


500 GeV 


800 GcV 


1000 GeV 


1500 GeV 


4000 GeV 


^max 


13.49 


38.60 


106.7 


149.7 


254.8 


775.9 


Bqooo 


1.508 


1.718 


1.716 


1.695 


1.668 


1.644 


-B2000 


0.2482 


0.6903 


1.258 


1.397 


1.534 


1.625 


-B1000 


0.8023 


1.499 


2.151 


2.288 


2.416 


2.496 


Bqow 


0.4011 


0.7495 


1.075 


1.144 


1.208 


1.248 


-B2OIO 


0.4011 


0.7495 


1.075 


1.144 


1.208 


1.248 


-B1010 


1.757 


2.408 


2.975 


3.092 


3.202 


3.269 


-B0101 


1.446 


1.486 


1.0749 


0.8790 


0.5962 


0.2235 


-Bnoi 


0.3496 


0.5187 


0.4593 


0.3889 


0.2717 


0.1037 


R 1 
^0000 


0.3259 


0.1900 


0.08785 


0.06979 


0.05347 


0.04103 


R 1 
^2000 


0.06009 


0.09431 


0.08469 


0.07414 


0.05908 


0.04258 


R 1 
^1000 


0.1911 


0.1996 


0.1407 


0.1183 


0.09133 


0.06506 


R 1 
^OOIO 


0.09556 


0.09979 


0.07034 


0.05916 


0.04567 


0.03253 


R 1 
- D 2010 


0.09556 


0.09979 


0.07034 


0.05916 


0.04567 


0.03253 


R 1 
-DlOlO 


0.3860 


0.2843 


0.1725 


0.1439 


0.1126 


0.08361 


R 1 
-DOIOI 


0.3099 


0.1540 


0.03373 


0.01373 


-0.001588 


-0.005863 


R 1 


0.08219 


0.06395 


0.01957 


0.009227 


0.0004544 


-0.002634 
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400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


4000 GeV 


o2 
°0000 


-0.002607 


-0.005957 


-0.005747 


-0.004551 


-0.002615 


-0.0005169 


q2 
°2000 


0.007821 


0.01787 


0.01724 


0.01365 


0.007846 


0.001551 


Q2 
°1000 


0.02134 


0.03293 


0.02514 


0.01910 


0.01057 


0.002038 


C2 
°0010 


0.01067 


0.01646 


0.01257 


0.009549 


0.005284 


0.001019 


C2 
°2010 


0.01067 


0.01646 


0.01257 


0.009549 


0.005284 


0.001019 


o2 

°1010 


0.005214 


0.01191 


0.01149 


0.009101 


0.005231 


0.001034 


q2 
°0101 


-0.005870 


-0.01676 


-0.02604 


-0.02585 


-0.02235 


-0.01181 


q2 
°1101 


0.007915 


0.005621 


-0.005706 


-0.008136 


-0.008968 


-0.005391 


tdSCA 

-"-oooo 


-0.001725 


-0.003455 


-0.003337 


-0.002677 


-0.001566 


-0.0003144 


jdSGA 
-"-2000 


0.03253 


0.02658 


0.01389 


0.009870 


0.005141 


0.0009552 


tdSGA 
-"-1000 


0.02733 


0.02246 


0.01183 


0.008418 


0.004393 


0.0008172 


tdSGA 

-"-ooio 


0.02733 


0.02246 


0.01183 


0.008418 


0.004393 


0.0008172 


tdSGA 
-"-2010 


0.02733 


0.02246 


0.01183 


0.008418 


0.004393 


0.0008172 


r>SGA 
-"-1010 


0.002977 


0.004972 


0.003879 


0.002952 


0.001636 


0.0003164 


tdSGA 
-"-0101 


-0.004042 


-0.01115 


-0.02368 


-0.02857 


-0.03614 


-0.05018 


jdSGA 
-"-1101 


0.02317 


0.01096 


-0.01227 


-0.02049 


-0.03195 


-0.04941 



SGAI 



V~s 


400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


4000 GeV 


^max 


10 GeV 


30 GeV 


90 GeV 


130 GeV 


230 GeV 


730 GeV 


-Boooo 


1.498 


1.690 


1.675 


1.653 


1.629 


1.610 


-B20OO 


0.2466 


0.6794 


1.228 


1.363 


1.498 


1.592 


-S1000 


0.7970 


1.475 


2.099 


2.232 


2.360 


2.444 


Bqoio 


0.3984 


0.7376 


1.050 


1.116 


1.180 


1.222 


-B2010 


0.3984 


0.7376 


1.050 


1.116 


1.180 


1.222 


-B1010 


1.745 


2.370 


2.903 


3.017 


3.127 


3.201 


-B0101 


1.436 


1.461 


1.047 


0.8566 


0.5817 


0.2187 


-Bnoi 


0.3472 


0.5104 


0.4480 


0.3790 


0.2651 


0.1015 


Rl 

^0000 


0.3157 


0.1621 


0.04600 


0.02769 


0.01430 


0.007045 


Rl 

^2000 


0.05846 


0.08341 


0.05453 


0.03989 


0.02330 


0.009018 


Rl 

^lOOO 


0.1858 


0.1758 


0.08901 


0.06213 


0.03493 


0.01351 


Rl 

^OOIO 


0.09292 


0.08791 


0.04451 


0.03107 


0.01746 


0.006754 


Rl 

- D 2010 


0.09292 


0.08791 


0.04451 


0.03107 


0.01746 


0.006754 


Rl 
-°1010 


0.3742 


0.2455 


0.1005 


0.06758 


0.03760 


0.01606 


Rl 
°0101 


0.3000 


0.1297 


0.007001 


-0.008684 


-0.01609 


-0.01073 


Rl 

-DllOl 


0.07988 


0.05563 


0.008266 


-0.0006056 


-0.006126 


-0.004889 
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G.2 e~ R e + 



t\ +x 



0(a s ) Analysis 



V~s 


400 GeV 


500 GcV 


800 GcV 


1000 GcV 


1500 GeV 


4000 GeV 


Coooo 


-5.088 


-13.73 


-28.17 


-34.56 


-46.82 


-83.45 


c*2000 


-0.9132 


-5.943 


-22.22 


-30.66 


-46.31 


-87.57 


C1000 


-3.387 


-14.68 


-42.75 


-56.41 


-81.91 


-151.7 


C0010 


-1.684 


-7.227 


-20.80 


-27.39 


-39.76 


-74.14 


C2010 


-1.704 


-7.467 


-22.11 


-29.27 


-42.59 


-78.47 


C1010 


-6.033 


-19.91 


-51.38 


-66.60 


-95.22 


-174.4 


C0101 


-4.816 


-11.64 


-17.53 


-18.03 


-17.21 


-12.20 


C1101 


-1.486 


-5.187 


-9.600 


-10.22 


-10.05 


-7.263 


-Doooo 


0.7129 


0.9271 


1.104 


1.159 


1.257 


1.521 


-D2000 


0.1308 


0.3992 


0.8305 


0.9708 


1.164 


1.505 


-D1000 


0.4875 


1.000 


1.638 


1.835 


2.115 


2.666 


Axuo 


0.2438 


0.5001 


0.8192 


0.9174 


1.057 


1.333 


-D2010 


0.2438 


0.5001 


0.8192 


0.9174 


1.057 


1.333 


-D1010 


0.8437 


1.3263 


1.9343 


2.130 


2.420 


3.025 


Ahoi 


0.6679 


0.7616 


0.6364 


0.5506 


0.4097 


0.1897 


-D1101 


0.2125 


0.3467 


0.3520 


0.3145 


0.2408 


0.1135 



v^ 


400 GeV 


500 GcV 


800 GeV 


1000 GeV 


1500 GeV 


4000 GeV 


-^0000 


0.5167 


0.6982 


0.7840 


0.7927 


0.7978 


0.7995 


-^2000 


0.09176 


0.2932 


0.5829 


0.6583 


0.7357 


0.7905 


D° 
-^1000 


0.3444 


0.7385 


1.152 


1.246 


1.338 


1.401 


D° 

-^0010 


0.1722 


0.3692 


0.5762 


0.6231 


0.6690 


0.7004 


-^2010 


0.1722 


0.3692 


0.5762 


0.6231 


0.6690 


0.7004 


D° 
-^1010 


0.6085 


0.9915 


1.367 


1.451 


1.534 


1.590 


D° 
-^0101 


0.4856 


0.5786 


0.4598 


0.3840 


0.2663 


0.1025 


-^1101 


0.1507 


0.2585 


0.2521 


0.2181 


0.1561 


0.06128 


5*0000 


0.1433 


0.08952 


0.04786 


0.04071 


0.03479 


0.03117 


5*2000 


0.02954 


0.04566 


0.03323 


0.02331 


0.008236 


-0.01002 


5iooo 


0.1078 


0.1127 


0.07325 


0.05653 


0.03448 


0.01029 


5ooio 


0.05401 


0.05750 


0.04221 


0.03599 


0.02805 


0.01947 


52010 


0.05381 


0.05506 


0.02959 


0.01821 


0.002443 


-0.01634 


5*1010 


0.1725 


0.1328 


0.07156 


0.05096 


0.02413 


-0.007143 


5oioi 


0.1322 


0.06490 


0.007412 


-0.003494 


-0.01250 


-0.01371 


5iioi 


0.04633 


0.03560 


0.007208 


0.000011398 


-0.006358 


-0.007873 
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v^ 


400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


4000 GeV 


(1 + ^o°ooo 


0.6613 


0.7901 


0.8313 


0.8307 


0.8266 


0.8204 


(1 + ^2°000 


0.1174 


0.3318 


0.6180 


0.6899 


0.7622 


0.8112 


(1 + «Mooo 


0.4408 


0.8357 


1.222 


1.306 


1.386 


1.437 


(1 + ^oVo 


0.2204 


0.4178 


0.6109 


0.6529 


0.6931 


0.7187 


(1 + ^2°010 


0.2204 


0.4178 


0.6109 


0.6529 


0.6931 


0.7187 


(1 + k)D\ 010 


0.7787 


1.122 


1.449 


1.521 


1.589 


1.632 


(1 + k)D° 0W1 


0.6216 


0.6547 


0.4875 


0.4024 


0.2759 


0.1052 


(1 + k)D» 101 


0.1929 


0.2925 


0.2673 


0.2285 


0.1617 


0.06289 


ol 
°0000 


-0.001287 


-0.002358 


0.0006297 


0.002742 


0.006079 


0.01024 


ol 
°2000 


0.003862 


0.007074 


-0.001884 


-0.008222 


-0.01824 


-0.03072 


ol 
°1000 


0.01143 


0.01556 


0.003835 


-0.003159 


-0.01367 


-0.02639 


ol 
°0010 


0.005814 


0.008914 


0.007500 


0.006143 


0.003972 


0.001128 


°2010 


0.005612 


0.006475 


-0.005121 


-0.01164 


-0.02164 


-0.03468 


ol 
°1010 


0.002248 


0.002290 


-0.01078 


-0.01854 


-0.03107 


-0.04878 


ol 
°0101 


-0.003764 


-0.01122 


-0.02028 


-0.02189 


-0.02208 


-0.01640 


ol 
°1101 


0.004154 


0.001587 


-0.007977 


-0.01043 


-0.01198 


-0.009478 



v^ 


400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


4000 GeV 


■QFull 

-"-oooo 


-0.001947 


-0.002985 


0.0007576 


0.003301 


0.007354 


0.01248 


T)Full 
-"-2000 


0.03289 


0.02132 


-0.003048 


-0.01192 


-0.02393 


-0.03787 


jdFuU 

-"-iooo 


0.02593 


0.01862 


0.003139 


-0.002419 


-0.009864 


-0.01836 


tdFuU 

-"-ooio 


0.02638 


0.02133 


0.01228 


0.009409 


0.005731 


0.001570 


jdFuU 
-^2010 


0.02546 


0.01550 


-0.008382 


-0.01782 


-0.03122 


-0.04826 


tdFuU 

-"-ioio 


0.002886 


0.002041 


-0.007438 


-0.01219 


-0.01955 


-0.02989 


tdFuU 

-"-oioi 


-0.006056 


-0.01714 


-0.04161 


-0.05439 


-0.08003 


-0.1559 


tdFuU 
-"-1101 


0.02154 


0.005427 


-0.02984 


-0.04566 


-0.07406 


-0.1507 
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Soft Gluon App. 

We define the u max to satisfy the following relation. 



dcosO J2 , T ' R n ( e L/R e+ -»• t 8 X(t or tg)) 
1 t i a cosy ' 



/ dcosO y^ 



dascA 



R, 



dcose y L ' R 



t s t) 





400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


4000 GeV 


uJmax 


13.55 


38.84 


107.2 


150.2 


255.2 


776.1 


Bqooo 


0.6600 


1.719 


1.717 


1.696 


1.668 


1.644 


-B2000 


0.1214 


0.6906 


1.259 


1.397 


1.535 


1.625 


-B1000 


0.4522 


1.500 


2.152 


2.289 


2.417 


2.496 


-60010 


0.2261 


0.7498 


1.076 


1.144 


1.209 


1.248 


-B20IO 


0.2261 


0.7498 


1.076 


1.144 


1.209 


1.248 


-B1010 


0.7814 


2.409 


2.976 


3.094 


3.203 


3.269 


-B0101 


0.6181 


1.486 


1.074 


0.879 


0.5964 


0.2235 


-B1101 


0.1970 


0.5189 


0.4596 


0.3891 


0.2718 


0.1037 


R 1 

^oooo 


0.1433 


0.1907 


0.08886 


0.07067 


0.05405 


0.04117 


R 1 
- D 2000 


0.02964 


0.09457 


0.08541 


0.07486 


0.05962 


0.04272 


R 1 
^lOOO 


0.1078 


0.2002 


0.1419 


0.1195 


0.09217 


0.06527 


R 1 
^0010 


0.05388 


0.1001 


0.07096 


0.05975 


0.04609 


0.03264 


R 1 
-02010 


0.05388 


0.1001 


0.07096 


0.05975 


0.04609 


0.03264 


R 1 
-OlOlO 


0.1729 


0.2853 


0.1743 


0.1455 


0.1137 


0.08389 


R 1 

-Ooioi 


0.1325 


0.1546 


0.03437 


0.01420 


-0.001372 


-0.005843 


R 1 
-OllOl 


0.04635 


0.06415 


0.01984 


0.009433 


0.0005521 


-0.002624 
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400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


4000 GeV 


C2 
°0000 


-0.001319 


-0.0103680 


-0.009239 


-0.007183 


-0.002033 


-0.0003753 


q2 
°2000 


0.003957 


0.01615 


0.01473 


0.01151 


0.008378 


0.001691 


Q2 
°1000 


0.01137 


0.02918 


0.02083 


0.01558 


0.01141 


0.002253 


C2 
°0010 


0.005683 


0.01459 


0.01042 


0.007792 


0.005703 


0.001126 


q2 
°2010 


0.005683 


0.01459 


0.01042 


0.007792 


0.005703 


0.001126 


o2 
°1010 


0.002638 


0.005783 


0.005487 


0.004327 


0.006345 


0.001315 


q2 
°0101 


-0.003445 


-0.02060 


-0.02827 


-0.02725 


-0.02214 


-0.011790 


q2 
°1101 


0.004170 


0.004308 


-0.006649 


-0.008751 


-0.008870 


-0.005381 


LfSGA 

-"-oooo 


-0.001994 


-0.005996 


-0.005351 


-0.004217 


-0.001217 


-0.0002282 


dSGA 

-"-2000 


0.03369 


0.02395 


0.01183 


0.008305 


0.005490 


0.001041 


TfSGA 

-"-1000 


0.02578 


0.01984 


0.009774 


0.006855 


0.004741 


0.0009033 


TfSGA 

-"-0010 


0.02578 


0.01984 


0.009774 


0.006855 


0.004741 


0.0009033 


TfSGA 
-n^oio 


0.02578 


0.01984 


0.009774 


0.006855 


0.004741 


0.0009033 


TfSGA 
-n-lOlO 


0.003387 


0.002406 


0.001847 


0.001401 


0.001985 


0.0004025 


TfSGA 

-"-0101 


-0.005543 


-0.01367 


-0.02564 


-0.03005 


-0.03579 


-0.05010 


TfSGA 
-"-1101 


0.02162 


0.008372 


-0.01426 


-0.02200 


-0.03161 


-0.04932 



SGAI 



V~s 


400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


4000 GeV 


^max 


10 GeV 


30 GeV 


90 GeV 


130 GeV 


230 GeV 


730 


-Boooo 


0.6554 


0.7741 


0.8080 


0.8077 


0.8057 


0.8031 


-B2000 


0.1206 


0.3350 


0.6106 


0.6787 


0.7476 


0.7950 


-S1000 


0.4492 


0.8384 


1.2035 


1.282 


1.358 


1.4086 


Bqoio 


0.2246 


0.4192 


0.6017 


0.6409 


0.6791 


0.7043 


-B2010 


0.2246 


0.4192 


0.6017 


0.6409 


0.6791 


0.7043 


-B1010 


0.7760 


1.109 


1.419 


1.486 


1.553 


1.598 


-B0101 


0.6138 


0.6349 


0.4629 


0.3802 


0.2592 


0.09768 


-Bnoi 


0.1957 


0.2901 


0.2568 


0.2177 


0.1526 


0.05847 


Rl 

^0000 


0.1388 


0.07592 


0.02394 


0.01493 


0.007892 


0.003678 


Rl 

^2000 


0.02883 


0.04174 


0.02772 


0.02035 


0.01190 


0.004559 


Rl 

^lOOO 


0.1047 


0.09992 


0.05103 


0.03569 


0.02010 


0.007784 


Rl 

^OOIO 


0.05237 


0.04996 


0.02552 


0.01784 


0.01005 


0.003892 


Rl 

- D 2010 


0.05237 


0.04996 


0.02552 


0.01784 


0.01005 


0.003892 


Rl 
-°1010 


0.1676 


0.1177 


0.05166 


0.03528 


0.01979 


0.008237 


Rl 
°0101 


0.1282 


0.05634 


0.003095 


-0.003854 


-0.007169 


-0.004792 


Rl 

-DllOl 


0.04502 


0.03162 


0.004739 


-0.0003479 


-0.003526 


-0.002817 
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G.3 Off-Diagonal Basis 



a 



e L e^ 



ij + X) /ctt,l 



v^ 


400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


Tree 


0.99876 


0.99255 


0.97350 


0.96514 


0.95345 


0(a s ) 


0.99850 


0.99093 


0.96815 


0.95833 


0.94480 


SGA 


0.99871 


0.99225 


0.97276 


0.96434 


0.95275 



a \e T e 



t i + X)/a T)L 



Vi 


400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


Tree 


0.00124 


0.00745 


0.02650 


0.03486 


0.04655 


0(a s ) 


0.00150 


0.00907 


0.03185 


0.04167 


0.05520 


SGA 


0.00129 


0.00775 


0.02724 


0.03566 


0.04725 



„(. 



e R e+ 



t i +X)/a T . 



\R 



V~s 


400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


Tree 


0.99833 


0.99371 


0.98141 


0.97619 


0.96897 


0{a s ) 


0.99786 


0.99184 


0.97615 


0.96958 


0.96057 


SGA 


0.99809 


0.99318 


0.98065 


0.97545 


0.96837 



a\e R e' 



t T + X) ja T 



\R 



V~s 


400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


Tree 


0.00167 


0.00629 


0.01859 


0.02381 


0.03103 


0(a s ) 


0.00214 


0.00816 


0.02385 


0.03042 


0.03943 


SGA 


0.00191 


0.00682 


0.01935 


0.02455 


0.03163 
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G.4 Helicity Basis 



t L + X) /g t ,l 



v^ 


400 GcV 


500 GcV 


800 GcV 


1000 GeV 


1500 GeV 


Tree 


0.66359 


0.75084 


0.83177 


0.84997 


0.86797 


0(a s ) 


0.66814 


0.75658 


0.83500 


0.85160 


0.86712 


SGA 


0.66809 


0.75657 


0.83579 


0.85299 


0.86962 



a(e L e + ^t R + X) /a T ^ L 



Vi 


400 GcV 


500 GcV 


800 GeV 


1000 GeV 


1500 GeV 


Tree 


0.33641 


0.24916 


0.16823 


0.15003 


0.13203 


0(a s ) 


0.33187 


0.24342 


0.16500 


0.14840 


0.13288 


SGA 


0.33191 


0.24343 


0.16421 


0.14701 


0.13038 



a 



e R e 



t R + X) far, 



\R 



V~s 


400 GeV 


500 GcV 


800 GcV 


1000 GeV 


1500 GeV 


Tree 


0.70979 


0.80926 


0.89267 


0.91039 


0.92758 


0{a s ) 


0.71528 


0.81541 


0.89546 


0.91145 


0.92608 


SGA 


0.71524 


0.81550 


0.89660 


0.91327 


0.92913 



a 



e R e 



t L + X) jo T , 



\R 



V~s 


400 GeV 


500 GcV 


800 GeV 


1000 GeV 


1500 GeV 


Tree 


0.29021 


0.19074 


0.10733 


0.08961 


0.072419 


0(a s ) 


0.28472 


0.18459 


0.10454 


0.08855 


0.073918 


SGA 


0.28476 


0.18450 


0.10340 


0.08673 


0.07087 
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a[e L e 



t T + X) /ct,l 



^ 


400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


Tree 


0.98811 


0.96769 


0.93097 


0.91822 


0.90218 


0(a s ) 


0.98851 


0.96784 


0.92922 


0.91565 


0.89853 


SGA 


0.98869 


0.96889 


0.93210 


0.91913 


0.90272 



a e r e 



t x + X) /cr T ,L 



V~s 


400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


Tree 


0.01189 


0.03231 


0.06903 


0.08178 


0.09782 


0(a s ) 


0.01149 


0.03216 


0.07078 


0.08435 


0.10147 


SGA 


0.01131 


0.03111 


0.06790 


0.08087 


0.09728 



a[e R e 



t l + X) /a T . 



^ 


400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


Tree 


0.99436 


0.98432 


0.96581 


0.95934 


0.95119 


0(a s ) 


0.99470 


0.98437 


0.96400 


0.95666 


0.94727 


SGA 


0.99490 


0.98545 


0.96693 


0.96025 


0.95174 



-„+ 



a\e R e 



t [ + X) /a T . 



v^ 


400 GeV 


500 GeV 


800 GeV 


1000 GeV 


1500 GeV 


Tree 


0.00564 


0.01568 


0.03419 


0.04066 


0.04881 


0(a s ) 


0.00530 


0.01563 


0.03600 


0.04334 


0.05273 


SGA 


0.00510 


0.01455 


0.03307 


0.03975 


0.04826 
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